Take Your Time

仕事や研究、コンピューターとの付き合い方

hist, lines, density, fitdistr()

マイクロデータの粗粗の分析方法。
1.hist()でヒストグラムを書く。
2.plot(density())でEDFを書いてみる。
3.fitdistr()でパラメータ推定してみる *libraryからMASS必須
5.2.で書いたグラフにline()でパラメタライズした分布関数を重ねて当てはまりを見る。このとき、y軸だけじゃなくてx軸も入力する。でないと変なことになる。

例:
1
> library(MASS)
> x<-rnorm(100)
> plot(density(x))
> x[5:10]<-NA
> plot(density(na.omit(x)))
> lines(dnorm(c(-3:3),mean=result$est[1],sd=result$est[2]),col="yellow")
> lines(c(-3:3),dnorm(c(-3:3),mean=result$est[1],sd=result$est[2]),col="red")