独立样本t检验
library(MASS)
t.test(Prob~So, data = UScrime,var.equal=TRUE)
t.test(Prob~So, data = UScrime,var.equal=FALSE)
![在这里插入图片描述](https://i-blog.csdnimg.cn/blog_migrate/b0341da2dd656d54a4235d42ba8a01ce.png)
配对样本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)
![在这里插入图片描述](https://i-blog.csdnimg.cn/blog_migrate/29f5e341bd7b527f799fb6ac1d049c48.png)
非参检验
wilcox.test(Prob ~So,data = UScrime)
with(UScrime,wilcox.test(U1,U2,paired = TRUE))
wilcox.test(UScrime$U1,UScrime$U2,paired = TRUE)
![在这里插入图片描述](https://i-blog.csdnimg.cn/blog_migrate/ba5cc60617ea252bbd0c7c106f9c1a34.png)
多余两组的比较
方差分析
state <- data.frame(state.region,state.x77)
library(multcomp)
model <- aov(Income ~state.region,state)
summary(model)
![在这里插入图片描述](https://i-blog.csdnimg.cn/blog_migrate/bbc2e0b5125ed625851eaf44e3819719.png)
非参检验
kruskal.test(Income ~state.region,state)
![在这里插入图片描述](https://i-blog.csdnimg.cn/blog_migrate/5b7ac8e630cfb2c34adec7534e00ad2a.png)