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

統計学備忘録(R言語のメモ)

since2016 ときどきTEXのメモ

2標本のグラフ、記述統計

2標本のまとめを簡単に

サンプル
dat1 <- c(-0.65, -0.58, 1.98, 0.98, 1.55, -0.05, 0.10, -1.35, 1.24, -0.28)
dat2 <- c(1.31, 2.39, 2.87, 2.25, 2.33, 5.32, 3.13, 1.21)

まずはデータセットを作成

data <- c(dat1, dat2)
group <- c(rep("A", length(dat1)), rep("B", length(dat2)))
df <- data.frame(group, data)

ノンパラの検定とグラフ
グラフの変数名A、Bを1回目、2回目に変更

wilcox.test(dat1, dat2, paired=F, exact =F, alternative="t")

boxplot(data ~ group, data=df, xaxt="n", xlab="")
axis(1,at=c(1,2),formatC(c("1回目","2回目")))
stripchart(data ~ group, data=df, vert=T, xaxt="n", xlab="", 
           method="jitter", jit=0.05, pch=16, cex=0.5, xlim = c(0.5,2.5))
axis(1,at=c(1,2),formatC(c("1回目","2回目")))

f:id:yoshida931:20200709182117p:plain:w500

記述統計
四分位、中央値、平均値、標準偏差、標準誤差

by(df$data, df$group, summary)

df$group: A
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
 -1.350  -0.505   0.025   0.294   1.175   1.980 
--------------------------------------------------------- 
df$group: B
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  1.210   2.015   2.360   2.601   2.935   5.320 

#install.packages("psych")
library(psych)
#関数describeを使用
by(df$data, df$group, describe, skew=F)

df$group: A
   vars  n mean   sd   min  max range   se
X1    1 10 0.29 1.09 -1.35 1.98  3.33 0.34
--------------------------------------------------------- 
df$group: B
   vars n mean   sd  min  max range   se
X1    1 8  2.6 1.29 1.21 5.32  4.11 0.45
‘‘‘