在R中查找直方图数据的p值

Question

根据随机化的斜率比例，大于或小于观察数据的斜率，我想计算获得观察到的斜率的预期概率。观察到的斜率为 -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 个大于原来的值。

Answer 1

如果您在 rand_vals 中有随机化程序的结果，在 obs_val 中有观察值，那么单尾 p 值（对零假设的量化支持与观测值大于零值的备择假设）是

mean(rand_vals>=obs)

请注意，这不是 ☢☣（找不到骷髅和交叉骨表情符号）"probability of getting the observed slope"。它是*如果原假设为真，则观察值 大于或等于观察到的斜率 的概率。
在某些情况下，将观测值也包含在 "randomization" 集中可能是合适的，即 mean(c(rand_vals,obs)>=obs)；如果您的随机化集很大，这不会有太大区别。
双尾 p 值类似于 mean(abs(rand_vals)>=abs(obs))

在R中查找直方图数据的p值

Finding p value of histogram data in R

for-loop

r

histogram

p-value