読者です 読者をやめる 読者になる 読者になる

Rを通じて統計学を学ぶ備忘録ブログ

SPSSからRに移行したい私のような人向けのR解説ブログ兼学習用備忘録。

第12回補足 1変数の要約

第12回補足は1変数の要約について書きます。

ヒストグラムと確率密度の視覚化を行っていきます。

 

par(mfrow=c(1,2)) #画面を1行2列に分割する
hist(demo$mathtest) #縦軸を「頻度」でヒストグラムを返す
hist(demo$mathtest,freq=FALSE) #縦軸を「密度」でヒストグラムを返す

f:id:teruaki-sugiura:20150712021455p:plain

sum(hist(demo$mathtest,plot=FALSE)$density*2) #確率密度であるため1が返る
[1] 1
densitymath <- density(demo$mathtest) #カーネル法による確率密度関数の推定
xlim <- range(densitymath$x) #xの範囲を指定
ylim <- range(densitymath$y) #yの範囲を指定
hist(demo$mathtest,probability=T,col="green",xlim=xlim,ylim=ylim)  
#確率密度を書く
lines(densitymath,lwd=3,col="black") #確率密度曲線を追加する

f:id:teruaki-sugiura:20150712022925p:plain

第12回補足は以上でお終い。

広告を非表示にする