理学療法士がまとめたノート

統計学備忘録 since2016

Rを使って統計学を勉強するブログです

記述統計

ダミー変数の作成 makedummies

下のようなカテゴリカルデータをダミー変数に変更します treat outcome NO NO NO YES YES NO NO YES YES YES YES YES NO YES YES NO NO NO NO YES NO NO NO NO YES YES YES YES YES NO YES NO NO NO NO NO YES YES NO NO YES NO YES NO YES YES NO YES NO YE…

平均と分散

平均と分散 投稿日2016.11.9 更新日2017.8.25 未だにピンときてない自分のために 更新します ( 1 , 2 , 3 , 4 , 5 , 6 ) がデータxの場合 x<-c( 1 , 2 , 3 , 4 , 5 , 6 ) データxの平均 sum( x ) / length( x ) = mean( x ) = 3.5 データxの二乗平均 sum(…

並行箱ひげ図

エクセルでデータのみコピーt<-read.table("clipboard")x<-t$V1y<-t$V2boxplot(x,y,xaxt="n",main="タイトル") #x軸のラベルを消去name<-c("ラベル1","ラベル2") #x軸のラベルを定義axis(side=1,at=c(1,2),labels=name) #x軸のラベルを更新 このような小技が…

データの要約

irisx<-iris$Sepal.Length[1:50]y<-iris$Sepal.Length[51:100]xy<-data.frame(x,y)summary(xy)boxplot(x,y)t<-c(1,2)ave<-c(mean(x),mean(y)) #平均値points(t,ave,pch=16) #平均値(・)を挿入 x y Min. :4.300 Min. :4.900 1st Qu.:4.800 1st Qu.:5.600 Me…

共分散と相関係数

(修正日2017.6.5)何度やっても忘れてしまう(悩)ので、大切なことは、真面目にコツコツ書くことにしましたたとえばエクセルからコピーしたデータをベクトルにする方法などx<-read.table("clipboard")x<-x$V1これすらも覚えられないので重症(´;ω;`)共分散デ…

確率とパーセント点

正規分布の累積分布(確率) pnorm(1.959964,lower=F) #Z≧1.959964となる上側確率=0.025 pnorm(-1.959964) #Z≦-1.959964となる上側確率=0.025 正規分布の確率点(パーセント点) qnorm(0.025,lower=F) #上側確率が0.025となるZ値=1.959964 qnorm(0.025) #…

幹葉図、箱ひげ

# 箱ひげ、パーセンタイル説明のために偏ったデータを作成 x<-c( 5, 8, 11, 23, 33, 35, 41, 44, 44, 44, 44, 46, 47, 49, 50, 51, 52, 52, 52, 53, 55, 55, 55, 56, 56, 56, 58, 58, 58, 59, 59, 59, 60, 60, 60, 61, 61, 62, 63, 63, 64, 65, 65, 65, 65, …