从λ= 0.10的指数分布生成三个大小为10的样本。进行F检验和Kruskal-Wallis检验,看三组中是否存在治疗差异。重复此操作1000次,记录两次测试的p值。通过比较由这些p值制成的两个直方图来比较模拟结果。
这是我到目前为止所拥有的:
x<-rexp(10, rate = 0.10)
y<-rexp(10, rate = 0.10)
z<-rexp(10, rate = 0.10)
x
y
z
A=c(x,y,z)
这将生成随机数。 我知道我应该有一些形式
results <-aov(A~groups, data=A)
summary(results)
但这不起作用。