SAS中文论坛

标题: 三组卡方检验两两比较无意义但可信区间不含0的理解 [打印本页]

作者: chenli9045 时间: 2014-6-5 11:54
标题: 三组卡方检验两两比较无意义但可信区间不含0的理解
各位大神，请教一下。
三组资料A、B、C有效率比较，三组比较P=0.0164，有统计学意义；
两两比较AB组P=0.0079，有统计学差异，BC组P=0.0321和AC组P=0.6919无统计学差异.
疑问是，做有效率组间率差的95%可信区间时，BC组的95%可信区间在0的一侧，是为什么？
如果没有统计学差异95%CI不是应该包含0吗？

作者: mono 时间: 2014-6-6 08:58
这种情况是可以出现的，因为你的数据不一定是正态分布的的，而且你的数据分布很不均衡。

作者: chenli9045 时间: 2014-6-6 09:04

mono 发表于 2014-6-6 08:58
这种情况是可以出现的，因为你的数据不一定是正态分布的的，而且你的数据分布很不均衡。 ...

我的变量取值是有效和无效。所以用的是卡方检验。卡方检验还要求正态分布吗？

作者: mono 时间: 2014-6-6 12:27
本帖最后由 mono 于 2014-6-6 12:33 编辑

chenli9045 发表于 2014-6-6 09:04
我的变量取值是有效和无效。所以用的是卡方检验。卡方检验还要求正态分布吗？ ...

95%可信区间是用正态分布来估计的，你问问你自己，为什么你的95%的可信区间总是对称的，你可以用profile 95%CI 来坐估计，这样你的95%的可信区间就不是对称的，卡放分布本身就不是对称的。
读读下面的东西，可能对你有益。
In general, the confidence interval based on the standard error strongly depends on the assumption of normality for the estimator. The "profile likelihood confidence interval" provides an alternative.
I am pretty sure you can find documentation for this. For instance, here and references therein.
Here is a brief overview.
Let us say the data depend upon two (vectors of) parameters, θ  and δ , where θ  is of interest and δ is a nuisance parameter.
The profile likelihood of θ  is defined by
L p (θ)=max δ L(θ,δ)
where L(θ,δ) is the 'complete likelihood'. L p (θ)  does no longer depend on δ since it has been profiled out.
Let a null hypothesis be H 0 :θ=θ 0 and the likelihood ratio statistic be
LR=2(logL p (θ ^ )−logL p (θ 0 ))
where θ ^ is the value of θ  that maximises the profile likelihood L p (θ)  .
A "profile likelihood confidence interval" for θ  consists of those values θ 0 for which the test is not significant.

作者: chenli9045 时间: 2014-6-6 14:35

mono 发表于 2014-6-6 12:27
95%可信区间是用正态分布来估计的，你问问你自己，为什么你的95%的可信区间总是对称的，你可以用profile ...

你好，谢谢回答。这样解释的时候会有问题吗？算是前后方法不一致得到了错误的结论吗？

作者: gogotiger 时间: 2014-6-7 22:28
BC组p=0.0321, BC组的95%可信区间在0的一侧，
这两者有矛盾吗？

作者: chenli9045 时间: 2014-6-9 08:58

gogotiger 发表于 2014-6-7 22:28
BC组p=0.0321, BC组的95%可信区间在0的一侧，
这两者有矛盾吗？

因为P值的结果不是应该和可信区间的方向是一致的吗？这里BC检验的时候不是调整了一类错误吗，所以P值结果是没有意义的，但是可信区间是有意义的。
不过我已经知道我的问题了。就是想知道可能会出现这样的结果吗？

作者: gogotiger 时间: 2014-6-23 00:50
你是单边检验的P值？

作者: chenli9045 时间: 2014-6-23 21:27

gogotiger 发表于 2014-6-23 00:50
你是单边检验的P值？

不是，是双侧检验的P值。

作者: gogotiger 时间: 2014-6-24 07:43
双边检验，p-value=0.03<0.05, 表明差异性。
CI 在同侧不含0，也表明差异性。
这没矛盾呀，还是我想的过于简单了？

作者: chenli9045 时间: 2014-6-24 10:31

gogotiger 发表于 2014-6-24 07:43
双边检验，p-value=0.03

因为我当时想的是0.03相当于是两两比较的结果，涉及多重比较就要调整P值，所以调整后变成没差异了，但是可信区间仍然显示有差异。
所以可能是因为这里有不对应的问题在，因为两两比较的时候相当于把P值的检验水准调成0.0167了，但是可信区间仍然计算的还是95%CI，所以前后是自相矛盾的了。

欢迎光临 SAS中文论坛 (https://mysas.net/forum/)