Censored pymc for power-law distribution of column density

Questions v5
, ,
3

I tried testing this using simulation. But this method is failing for highly scatter dataset. Also if the number of input data are less then also it is giving a very bad fit. Any comment?