s**f
发帖数: 365 | 1 请大家帮个忙，谢谢了！
我现在想要做个图（不是家庭作业！），同时算一下simple linear regression的参数
，并test一下斜率是否为0.数据是下面的格式：
Cancer Race Gender AgeGrp Year_dx
0 White F 1 2002
0 White M 3 1998
1 White M 2 2003
0 Black M 1 1997
. . .
. . .
. . .
cancer 0是没有，1是有
agegrp 1是20-44，2是45-64，3是65+
year_dx是1998，1999，。。。，2010
请问如何用SAS很简洁做下面的图？我能想到的方法可能太繁琐了。
X轴是年份（1980-2010）
Y轴是white male的cancer的rate（数量，就是把cancer中的1都加起来）
在上面这个图中画三条曲线（3个年龄组），每... 阅读全帖 |
|
l**********9
发帖数: 148 | 2 高手的代码果然简练，佩服佩服...
不过我记得sgplot中的reg不能输出estimate parameter，想要检验斜率只能用肉眼了.
...如果想输出的话，还是得用proc reg吧...
proc sort data=one;
by Agegrp;
run;
proc reg data=one outest=temp;
by Agegrp;
model cancer = year_dx ;
run;
data slope(drop=_MODEL_ _TYPE_ _DEPVAR_ _RMSE_ Intercept cancer);
set temp(rename=(year_dx=slope));
run;
这样能显示斜率...但是代码好长不美观>_< |
|
J********J
发帖数: 571 | 3 Thanks for the reply. I was asking how to imput "≤", "≥ " as values, not
an expression.
For example; if age >= 50 then agegrp = "Age ≥ 50 ";
Anybody knows? |
|
C*******d
发帖数: 80 | 4 You can code as: if age >= 50 then agegrp = 1; etc.
Then use proc format to display 1 as "Age ≥ 50 "
Hope it helps. |
|
