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

統計学備忘録 since2016

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

記述統計

平均と分散

平均と分散 投稿日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, …