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, 1.50, 0.97)

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

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

par(mfrow=c(1, 2))#1行2列の枠を作成

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回目"))
)

par(mfrow=c(1, 1))#1行1列に戻しとく

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

library(psych)
describe(df$data, skew=F)
by(df$data, df$group, describe, skew=F)
#skew:skew and kurtosis=歪度と尖度