java计算R值,R：使用模拟计算p值

我编写了这段代码来对两个随机分布的观察x和y运行一个测试统计

mean.test

alternative=c("two.sided","less","greater"))

{

p.value

alternative

s

t

p.value

TRUE, log.p = FALSE)) #try to calculate p value

data.name

names(t)

zero

names(zero)

return(structure(list(statistic = t, p.value = p.value,

method = "mean test", data.name = data.name,

observed = c(x,y), alternative = alternative,

null.value = zero),

class = "htest"))

}

代码使用蒙特卡罗模拟来生成检验统计量均值(x) - 均值(y)的分布函数，然后计算p值，但显然我错过了定义这个p值，因为：

> set.seed(0)

> mean.test(rnorm(1000,3,2),rnorm(2000,4,3))

输出应该如下：

mean test

data: c(rnorm(1000, 3, 2), rnorm(2000, 4, 3))

difference in means = -1.0967, p-value < 2.2e-16

alternative hypothesis: true difference in means is not equal to 0

但我得到了这个：

mean test

data: c(rnorm(1000, 3, 2), rnorm(2000, 4, 3))

difference in means = -1.0967, p-value = 0.8087

alternative hypothesis: true difference in means is not equal to 0

有人可以向我解释这个错误吗？