基于来自随机化的斜率的比例,大于或小于观察到的数据的斜率,我想计算获得观察到的斜率的预期概率。观察到的斜率为-0.2717。
非常感谢任何帮助,我是新手。
histdata<- numeric(10000)
for (i in 1:10000) {histdata[i]<-(summary.lm(lm(sample(tcons)~tleave))
[[4]][[2]])}
hist(histdata)
abline(v=-0.2717, lwd=3, lty=2)
box()
data3<- -0.2717>histdata
这^^给了我9954不大于原来的46和46更大。
答案 0 :(得分:0)
如果您有rand_vals
中的随机化过程的结果和obs_val
中的观察值,那么单尾 p值(量化对null的支持)假设与观察值大于空值的替代假设是
mean(rand_vals>=obs)
mean(c(rand_vals,obs)>=obs)
;如果您的随机化集很大,这将不会有太大的影响。mean(abs(rand_vals)>=abs(obs))