arupaka-_-arupaka https://blog.hatena.ne.jp/arupaka-_-arupaka/ arupaka-_-arupakaの日記 https://arupaka-arupaka.hatenadiary.org/ べき分布 R データ解析 べき分布フィッティングの自動化。べき指数の推定のKS testと最尤法を使う方法はかなり使える。 コードは、最後のフィッティングの補正の箇所。 h1d=1 & b1b==code]plot(sort(h1d),length(h1d):1/length(h1d),log="xy",type="l") alphaxminh1e=xmin & b1b==code] ph1curve(ph1*xmin^alpha/((x)^alpha),add=T,col=2)#titleをプロット title(paste(sprintf("%.3f",alpha),sprintf("%d",as.integer(x… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Farupaka-arupaka.hatenadiary.org%2Fentry%2F20090819%2F1250670034" title="べき分布フィッティングの自動化 - arupaka-_-arupakaの日記" class="embed-card embed-blogcard" scrolling="no" frameborder="0" style="display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;"></iframe> Hatena Blog https://hatena.blog 2009-08-19 17:20:34 rich https://arupaka-arupaka.hatenadiary.org/entry/20090819/1250670034 1.0 100%