R
!!!fitdistrplus
分布のあてはまりの確認

 install.packages("fitdistrplus")
 library(fitdistrplus)

!!descdist()
分布の描写
* これで、どんな分布か、あたりをつける
 descdist(SbyS.dat.wide$MDD)
{{ref_image descdist.png}}
! boot オプション
* ブートストラップしたサンプルも追加でプロット
 descdist(SbyS.dat.wide$MDD, boot=500)
{{ref_image descdistboot.png}}
!!fitdist()
* 分布を当てはめてみる
 fitdist(SbyS.dat.wide$MDD, "gamma")

* 結果をプロットする
{{pre
MDD.gamma.fit <- fitdist(SbyS.dat.wide$MDD, "gamma")
plot(MDD.gamma.fit)
}}





!!Reference
https://www.jstatsoft.org/article/view/v064i04