SAS中如何实现多样本趋势检验？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

[tr][/tr]

现观察三组人，他们都是男性。第一组人的年龄在20到30岁之间，第二组人的年龄在30到40岁之间，第三组人的年龄在40到50岁之间。他们的β脂蛋白的测量值如下表，问这三组人的测量值是否符合人们的经验：年龄越大的人的β脂蛋白的含量越大？ 第  1 组 260 200 240 170 270 205 190 200 250 200 第  2 组 310 310 190 225 170 210 280 210 280 240 第  3 组 320 260 360 310 270 380 240 295 260 250 参考了一些资料，发现是要写程序的，网上有相关代码，但是也不是很清楚具体要怎么写。 哪些数据是要提前计算的，哪些数据是要带入程序替换的？

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：如何实现 测量值 脂蛋白 网上有 如何 脂蛋白 样本 测量 程序

对于分类变量的趋势检验，我们可以用卡方检验，即proc freq，对于连续性变量，可以用线性回归模型检验其trend，sas code as follows:

1. data test;
   
2. input group result @@;
   
3. cards;
   
4. 1 260 1 200 1 240 1 170 1 270 1 205 1 190 1 200 1 250 1 200
   
5. 2 310 2 310 2 190 2 225 2 170 2 210 2 280 2 210 2 280 2 240
   
6. 3 320 3 260 3 360 3 310 3 270 3 380 3 240 3 295 3 260 3 250  
   
7. ;
   
8. run;
   
9. proc reg data=test;
   
10. model result = group;
    
11. run;

复制代码

heperwong 发表于 2014-1-7 17:57
对于分类变量的趋势检验，我们可以用卡方检验，即proc freq，对于连续性变量，可以用线性回归模型检验其tre ...

太感谢啦  热心同学   我先去操作下试试  [em24][em24]

稍微补充以下:楼上的建议之外,也可以考虑检测线性趋势. 从数据结果来看,趋势存在并且明显.京剧
see more about R code-
http://blog.sina.com.cn/s/blog_a3a926360101hv59.html

1. data have;
   
2.   infile datalines missover;
   
3.   input group @;
   
4.   do i =1 to 10;
   
5.     input y @;
   
6.     output;
   
7.     end;
   
8. datalines;
   
9. 1  260 200 240 170 270 205 190 200 250 200
   
10. 2  310 310 190 225 170 210 280 210 280 240
    
11. 3  320 260 360 310 270 380 240 295 260 250
    
12. ;
    
13. proc glm data =have;
    
14.   class group;
    
15.   model y =group;
    
16.   contrast 'linear trend' group -1 0 1;
    
17.   run; quit;

复制代码

jingju11 发表于 2014-1-8 00:59
稍微补充以下:楼上的建议之外,也可以考虑检测线性趋势. 从数据结果来看,趋势存在并且明显.京剧
see more a ...

你好！这句（contrast 'linear trend' group -1 0 1;）是什么意思啊！麻烦解释一下好吗？

jingju11 发表于 2014-1-8 00:59
稍微补充以下:楼上的建议之外,也可以考虑检测线性趋势. 从数据结果来看,趋势存在并且明显.京剧
see more a ...

你好 谢谢你的热心解答  ，还有个问题  我的数据标准差都很大    比较离散 ， 需要把各组内的数据转化成正太分布吗？赶脚转化成正太貌似不是很合适吖   ~

heperwong 发表于 2014-1-7 17:57
对于分类变量的趋势检验，我们可以用卡方检验，即proc freq，对于连续性变量，可以用线性回归模型检验其tre ...

你好 能不能再麻烦你解释一下分类变量的趋势检验如何做   谢谢你

lmq707 发表于 2014-1-8 09:03
你好！这句（contrast 'linear trend' group -1 0 1;）是什么意思啊！麻烦解释一下好吗？

是做一个线性趋势检验，“”里的内容是标签，group是用于分组的变量，-1 0 1依次是各组系数，有三个数字应该是group分了三组，趋势检验注意系数和要为零，系数之间的差值要相等。个人的一点经验，有错的地方还请各位指正。

jdc546576585 发表于 2015-1-24 14:46
是做一个线性趋势检验，“”里的内容是标签，group是用于分组的变量，-1 0 1依次是各组系数，有三个数字应 ...

不好意思，我想问下，contrast 后面的系数-1 0 1 是什么意思，如果-1 0 1 代表线性趋势，那么如果我要第一组和第三组比较，后面的系数该怎么写？

Emily225928 发表于 2014-1-8 18:48
你好 能不能再麻烦你解释一下分类变量的趋势检验如何做   谢谢你

分类变量的趋势检验，对二元变量（表现为2×k的列联表），可以使用Cochran-Armitage test

1. /* SAS: proc freq + trend option*/
   
2. proc freq data = test;
   
3.     tables X * Y / Trend;
   
4.     weight count;
   
5. run;

复制代码

对多元变量（m*k）,可以使用Jonckheere-Terpstra非参数检验，SAS中proc freq + JT option。补充，对二元变量，也可以用proc logistic中的score test实现。