全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 SAS专版
5276 7
2010-02-08
问题是这样的
有一个列联表
       男    女
挨个āi   46  341
挨个ái   25   2
现在想知道不同发音在同一个性别上是否存在显著差异,比如性别为男时,发一音和发二音在置信水平0.05时有没有差异。我想用单因素方差分析,是这样理解的A因素为发音,一共有两个水平,分别为发一音和发二音,B因素只有一个水平,性别为男 来进行检验A因素中两个水平的差异。我用SAS的anova语句检验,输出的好多统计量为".",包括p值,是有问题的有人告诉我可以考虑一下加权重,用weight语句,但是我翻了好多书,没有见到在方差分析过程中加weight语句的资料,我想问这样的问题怎么解决呢,谢谢
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2010-2-8 11:32:33
单因素方差分析属于定量统计分析;列联表属于定性统计分析。不能混在一起用。

下面代码中,sex-1表示男;sex-2表示女,pronounce(发音)1表示1音;2表示2音。

DATA ab;
DO sex=1 TO 2;
DO pronounce=1 TO 2;
INPUT count @@;
OUTPUT;
END;
END;
CARDS;
46 25
341 2
;
PROC FREQ;
WEIGHT count;
TABLES sex*pronounce / CHISQ;
RUN;
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2010-2-8 11:57:13
你好,谢谢你能解答我的问题,不过你写的代码是列联表的独立性检验的代码,相当于检验性别与发音是否独立。我想知道的问题是比如性别为男时,发一音和发二音在置信水平0.05时有没有差异,性别为女时,发一音和发二音在置信水平0.05时有没有差异,分开来考虑,有没有比较好的办法?谢谢你
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2010-2-8 14:08:42
根据你的设计,只能采用相对危险度(Relative Risk, RR)去检验,
sex-1表示男;sex-2表示女,pronounce(发音)1表示1音;2表示2音
DATA ab;
DO pronounce=1 TO 2;
DO sex=1 TO 2;
INPUT count @@;
OUTPUT;
END;
END;
CARDS;
46 341
25 2
;
PROC FREQ;
WEIGHT count;
TABLES pronounce*sex/ all;
RUN;

                                                                pronounce*sex的汇总统计量

                                           Cochran-Mantel-Haenszel 统计量(基于表得分)

              统计量    对立假设                  自由度          值      概率
                                 ----------------------------------------------------------------
                                   1       非零相关                       1    115.4216    <.0001
                                   2       行均值得分差值            1    115.4216    <.0001
                                   3       一般关联                       1    115.4216    <.0001

                                                               普通相对风险的估计值(行 1/行 2)

                              研究类型             方法                             值         95% 置信限
                           ---------------------------------------------------------------------------
                           案例对照             Mantel-Haenszel        0.0108       0.0025      0.0471
                             (优比)           Logit                           0.0108       0.0025      0.0471
                           Cohort               Mantel-Haenszel          0.1284       0.0959      0.1718
                             (第 1 列风险)    Logit                      0.1284       0.0959      0.1718
                           Cohort               Mantel-Haenszel        11.8953       3.1332     45.1611
                             (第 2 列风险)    Logit                    11.8953       3.1332     45.1611



                                                         总样本大小 = 414

你的抽样样本部分偏少,最好重新调整。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2010-4-13 08:05:05
从这个结果总,如果书面结论该怎么说呢??

一直对这个迷糊,解释不清
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2010-4-13 19:16:43
1# nankaimy
可不可以用卡方检验呢,就是检验两个水平的人数是不是各占50%。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群