采用具有块随机效应的线性模型 ( ) 和等。
如果,则模型简化为非常简单的线性模型。要检验原假设与替代,可以使用似然比检验。但是,null 下的极限分布不是,而是,即 50:50 的混合零点质量和一个自由度的卡方分布。
我们如何确定该混合物的临界值(或使用的 p 值)?
采用具有块随机效应的线性模型 ( ) 和等。
如果,则模型简化为非常简单的线性模型。要检验原假设与替代,可以使用似然比检验。但是,null 下的极限分布不是,而是,即 50:50 的混合零点质量和一个自由度的卡方分布。
我们如何确定该混合物的临界值(或使用的 p 值)?
0 以上任何点右侧的区域是的一半。因此,要获得水平测试,请查找点。
....只要。
当然,p 值的工作方式类似。将值向上查找,就好像它是并将生成的 p 值减半。`
如果有兴趣,一些 R 代码:
使用 TcGSA 包:
ss <- 0.3
sample_mixt <- TcGSA:::rchisqmix(n=1e5, s=0, q=1)
TcGSA:::pval_simu(s=ss, sample_mixt)
使用基地:
1/2*(1-pchisq(ss,df=1))