SAS中文论坛

标题: 请教3个有关SAS的小case！ [打印本页]

作者: shiyiming 时间: 2006-4-26 14:09
标题: 请教3个有关SAS的小case！
大家好！
很高兴加入大家庭！我有几个问题想请求帮忙：
1.关于样本量的问题。在我的实验研究中要用到线性相关分析的。在收集的34份标本中，每份标本检测出有5种物质（A，B，C，D，E）的含量。现在我需要观察数据A分别与数据B，C，D，E的线性关系。即：A对B，A对C，A对D
，A对E。那么，我现在的34份标本足够吗？如果不够的话，应该是多少份才合适呢？为什么？
2.以下是我用SAS编程的两个程序，请检查是否正确。若有问题请指点！谢谢！

A 直线相关spearman编程
data ratio_CHOL;
input y x @@;
cards;
7.65       7.40
7.33       5.80
5.99       6.80
5.71       7.50
4.87       6.20
4.65       5.80
4.51       7.30
4.50       4.80
4.31       5.90
4.25       5.20
3.96       4.80
3.95       6.00
3.76       4.40
3.69       5.50
3.56       5.60
3.32       6.30
3.27       4.50
3.08       5.50
3.04       3.70
3.04       3.30
3.02       5.90
2.85       4.60
2.74       4.70
2.70       4.90
2.68       6.00
2.60       6.20
2.59       6.00
2.23       4.40
2.14       4.10
2.09       5.20
2.09       4.00
2.04       4.10
1.71       6.20
1.66       4.40
1.37       5.80
;
run;
proc gplot data=ratio_CHOL;
plot y*x;
run;
proc corr spearman data=ratio_CHOL;
var x y;
run;
proc reg data=ratio_CHOL;
model y=x;
run;

B 当我的数据需要用到Wilcoxon符号秩和检验的时候，SAS的编程该如何呢？我的样本属于是单一样本与总体中位数的比较。请看我的编程：
date ratio_1;
input x1 x2@@;
d=x1-x2;
cards;
7.65  1.00
7.33  1.00
5.99  1.00
5.71  1.00
4.87  1.00
4.65  1.00
4.51  1.00
4.50  1.00
4.31  1.00
4.25  1.00
3.96  1.00
3.95  1.00
3.76  1.00
3.69  1.00
3.56  1.00
3.32  1.00
3.27  1.00
3.08  1.00
3.04  1.00
3.04  1.00
3.02  1.00
2.85  1.00
2.74  1.00
2.70  1.00
2.68  1.00
2.60  1.00
2.59  1.00
2.23  1.00
2.14  1.00
2.09  1.00
2.09  1.00
2.04  1.00
1.71  1.00
1.66  1.00
1.37  1.00
;
proc print;
proc univariate normal;
var d;
run;

请大家帮帮忙！万分感激！在线等！
,

作者: shiyiming 时间: 2006-4-26 16:30
标题: re
呵呵呵，欧也菜鸟一个。弱弱的回答：
据我的理解，问题b，lz编的程序是对的亚，结果看S的p值。
不过既然加了normal选项，应该先看是否正态分布，如果不是，就看S对应的p值，否则可以用配对t检验的说。

作者: shiyiming 时间: 2006-4-27 09:19
标题: to greenhand2008
第二题，你要看数据是否正态，然后不正态用nap1way里的秩和检验

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