独立样本t检验
library(MASS)
t.test(Prob~So, data = UScrime,var.equal=TRUE)
t.test(Prob~So, data = UScrime,var.equal=FALSE)
配对样本t检验
library(MASS)
sapply(UScrime[c("U1","U2")],function(x)(c(mean=mean(x),sd=sd(x))))
with(UScrime,t.test(U1,U2,paired = TRUE))
t.test(UScrime$U1,UScrime$U2,paired = TRUE)
非参检验
wilcox.test(Prob ~So,data = UScrime)
with(UScrime,wilcox.test(U1,U2,paired = TRUE))
wilcox.test(UScrime$U1,UScrime$U2,paired = TRUE)
多余两组的比较
方差分析
state <- data.frame(state.region,state.x77)
library(multcomp)
model <- aov(Income ~state.region,state)
summary(model)
非参检验
kruskal.test(Income ~state.region,state)