楼主: lnlhckao123
3882 6

[原创博文] 原因变量为有序变量的R*2列联表资料检验方法的选择 [推广有奖]

  • 0关注
  • 3粉丝

副教授

62%

还不是VIP/贵宾

-

威望
0
论坛币
14810 个
通用积分
10.6125
学术水平
7 点
热心指数
4 点
信用等级
4 点
经验
12566 点
帖子
754
精华
0
在线时间
405 小时
注册时间
2010-8-21
最后登录
2025-6-7

楼主
lnlhckao123 发表于 2011-6-27 00:51:56 |AI写论文
200论坛币
我在胡良平编的书上看到一段话“对于R*2列联表,如果原因变量为多值有序变量要看分析的目的,如果要分析原因变量间的频数分布是否相同,则可用一般卡方检验或Fisher精确检验;如果要分析第一列(或第二列)上的频率变化是否呈线性趋势,则要选用线性趋势检验”。我想请教,怎样理解频数分布是否相同与“频数变化是否呈线性趋势”,谢谢高手指点!

最佳答案

stata18 查看完整内容

data prg15_1; do dose=1 to 4; do response=0 to 1; input f@@; output; end; end; cards; 120 30 100 50 90 60 50 100 ; run; proc freq;weight f;tables dose*response/trend nocol nopct;run; proc freq;weight f;tables dose*response/chisq nocol nopct;run;
关键词:检验方法 因变量 列联表 Fisher Fish 因变量 资料

沙发
stata18 发表于 2011-6-27 00:51:57
data prg15_1;                           
do dose=1 to 4;
  do response=0 to 1;
   input f@@;
   output;
  end;
end;
cards;
120   30
100   50
90   60
50  100
;
run;
proc freq;weight f;tables dose*response/trend nocol nopct;run;
proc freq;weight f;tables dose*response/chisq nocol nopct;run;

藤椅
lnlhckao123 发表于 2011-7-5 01:21:30
为了高手理解我的意思,我具体列出程序,请高手指点!程序如下:
为了研究某种化合物的致癌性,将60只大鼠随机分为4组,分别接受100mg、200mg、300mg和400mg剂量的该化合物 饲料喂养,随访观察1年以后观察其肿瘤发生情况,实验数据见下表,欲考察肿瘤的发生率是否随着剂量水平的变化而呈线性趋势。
                                             60只大鼠肿瘤发生情况
化合物剂量                                                   例数
                               肿瘤发生情况             未发生              发生                       合计
100mg                                                         12                   3                           15                 
  200mg                                                       10                   5                            15
  300mg                                                        9                    6                            15  
400mg                                                          5                  10                            15
合计                                                           36                  24                            60


检验上述数据显性趋势的程序如下
data prg15_1;                           
do dose=1 to 4;
  do response=0 to 1;
   input f@@;
   output;
  end;
end;
cards;
12  3      
10  5
  9  6
  5  10
;
run;

proc freq data=prg15_1;
tables dose*response/trend norow nocol  nopercent
scores = table;
weight f;
exact trend;
ods html;
run;  

运行上述结果如下,显示结果P(0.0132)<0.05,

                    Statistics for Table of dose by response

                                              Cochran-Armitage Trend Test
                                              ----------------------------
                                              Statistic (Z)        -2.5927

                                              Asymptotic Test
                                              One-sided Pr <  Z     0.0048
                                              Two-sided Pr > |Z|    0.0095

                                              Exact Test
                                              One-sided Pr <=  Z    0.0066
                                              Two-sided Pr >= |Z|   0.0132

而当用普通的卡方检验,如下
data ex7_1;                              
do dose=1 to 4;
  do response=1 to 2;
   input f@@;
   output;
  end;
end;
cards;
12  3      
10  5
  9  6
  5  10
;
run;
proc freq;
  weight f;
  tables dose*response
  /chisq
  expected;
run;
结果如下,P(0.06514216730615)>0.05,这是怎么回事,请高手赐教!!
                   Statistics for Table of dose by response

                            Statistic                                     DF        Value                Prob
                            ----------------------------------------------------------------
                            Chi-Square                                 3        7.2222       0.06514216730615
                            Likelihood Ratio Chi-Square       3        7.3681       0.06104542043560
                            Mantel-Haenszel Chi-Square      1         6.6102      0.01013971166655
                            Phi Coefficient                                       0.3469
                            Contingency Coefficient                        0.3278
                            Cramer's V                                            0.3469
即使在人大经济论坛这个网络世界,我仍以真诚为基础与我的好友进行交往!

板凳
stata18 发表于 2011-7-5 13:10:39
谈点个人看法:
  这个应该使用趋势检验;  趋势检验结果p<0.05,可解释为“肿瘤发病率有随着剂量增加而增加的趋势”;如果用普通卡方检验,不考虑自变量的有序性,那么,
12   3
10   5
9   6
5  10
==============
9   6
5  10
12   3
10   5
没什么不同,卡方检验结果相同;但是趋势卡方检验结果不同。
只是普通卡方检验结果>0.05确实不太好解释,也许是样本量不够?

报纸
stata18 发表于 2011-7-5 13:11:57
趋势卡方在线性趋势存在的情况下更敏感?

地板
stata18 发表于 2011-7-5 13:13:24
似乎确实与样本量有一定关系,此样本数*10,卡方检验结果就<0.05了

7
stata18 发表于 2011-7-5 13:18:41
”频数分布是否相同”应该是各剂量组的发病率有无差异;无效假设:rate1=rate2=rate3=rate4
“频数变化是否呈线性趋势”应该是rate1<rate2<rate3<rate4或者rate1>rate2>rate3>rate4,大概是这样理解吧

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2026-1-4 23:24