按楼主的题设,50个样本,u=170, sd = 100, 500, 1000 , 分别模拟一百次,每次模拟生成50个随机数字并以mu=165做t检验计算P值,三个条件下P值结果汇总在图里
结论是P区别没有 😉
library(ggplot2)
test1= sapply(1:100,function(i){
x=rnorm(n = 50,mean = 170,sd=100)
p=t.test(x,mu = 165)$p.value
return(p)
})
test2= sapply(1:100,function(i){
x=rnorm(n = 50,mean = 170,sd=500)
p=t.test(x,mu = 165)$p.value
return(p)
})
test3= sapply(1:100,function(i){
x=rnorm(n = 50,mean = 170,sd=1000)
p=t.test(x,mu = 165)$p.value
return(p)
})
df=data.frame(
group=c(rep("sd100",100),rep("sd500",100),rep("sd1000",100)),
p=c(test1,test2,test3)
)
df %>% ggplot(aes(x=group,y=p))+
geom_boxplot()