rp100 <- function(){
set.seed(1234)
n=50
x1=runif(n);x1
x2=runif(n);x2
x3=runif(n);x3
x4=runif(n);x4
y=10*x1+10*x2+10*x3+10*x4+rnorm(n,0,1);y
x=cbind(x1,x2,x3,x4)
M1=lm(y~x);summary(M1)
e=resid(M1);e
yy=y-e;yy
Fm=lm(e^2~yy+yy^2)
summary(Fm)
R2=summary(Fm)$r.sq
m=ncol(Fm$model)
W=n*R2;W
P=1-pchisq(W,m-1);P
return(P)
}
# 设置一个新的变量
x <- vector()
# 重复10000最后的x就是你需要的结果
for(i in 1:1000)
x[i] <-rp100()
x
x[x<0.05]