零点质量的临界值和一个自由度的卡方分布

机器算法验证 混合模式 混合分布
2022-04-14 12:26:50

采用具有块随机效应的线性模型 ( ) 等。bi

yij=μ+bi+βxij+eij
biN(0,σb2)

如果,则模型简化为非常简单的线性模型。要检验原假设与替代,可以使用似然比检验。但是,null 下的极限分布不是,而是,即 50:50 的混合零点质量和一个自由度的卡方分布。σb2=0H0:σb2=0H1:σb2>0χ1212χ02+12χ12

我们如何确定该混合物的临界值(或使用的 p 值)?

2个回答

0 以上任何点右侧的区域是的一半。因此,要获得水平测试,请查找χ12α2αχ12

....只要α<0.5

当然,p 值的工作方式类似。将值向上查找,就好像它是并将生成的 p 值减半。`χ12

如果有兴趣,一些 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))