べき分布フィッティングの自動化

べき分布フィッティングの自動化。

べき指数の推定のKS testと最尤法を使う方法はかなり使える。

コードは、最後のフィッティングの補正の箇所。

h1d<-h1b[h1b>=1 & b1b==code]

plot(sort(h1d),length(h1d):1/length(h1d),log="xy",type="l")
alpha<-alpha123b[which(Db==min(Db))]

xmin<-xminbb[which(Db==min(Db))]

h1e<-h1b[h1b>=xmin & b1b==code]
ph1<-length(h1e)/length(h1d)

curve(ph1*xmin^alpha/((x)^alpha),add=T,col=2)

#titleをプロット
title(paste(sprintf("%.3f",alpha),sprintf("%d",as.integer(xmin))))

(本質的な部分の
　コードは、結構簡単で、実装済み..。)。