2

我希望在 R 中做一些基本的模拟来检查 p 值的性质。我的目标是查看大样本量是否趋向于小 p 值。我的想法是生成 1,000,000 个数据点的随机向量，将它们相互回归，然后绘制 p 值的分布并寻找偏斜。

到目前为止，这是我的想法：

x1 = runif(1000000, 0, 1000) 
x2 = runif(1000000, 0, 1000) 
model1 = lm(x2~x1)

使用取自另一个线程的代码：

lmp <- function (modelobject) {
    if (class(modelobject) != "lm") stop("Not an object of class 'lm' ")
    f <- summary(modelobject)$fstatistic
    p <- pf(f[1],f[2],f[3],lower.tail=F)
    attributes(p) <- NULL
    return(p)
  }
lmp(model1) 
0.3874139

关于如何为 1000 个模型甚至更多模型执行此操作的任何建议？谢谢！

4

1 回答 1

1

看到?replicate...但是您正在计算的 p 值假设高斯误差不是统一的（这在 n=10^6 时并不重要）

具体来说，是这样的：

nrep <- 1000
ndat <- 1000000
results <- replicate(nrep, {
     x1=runif(ndat, 0, 1000);
     x2=runif(ndat, 0, 1000);
     model1=lm(x1 ~ x2);
     lmp(model1)
     })

应该可以工作，但是运行需要很长时间。

我建议将 nrep 和 ndat 做得更小以尝试一下。

于 2016-04-24T05:47:03.063 回答