求助，关乎行列均为顺序变量的线性趋势检验程序问题？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

   题目是这样的：
为了研究晶状体混浊程度是否与年龄相关，将资料整理如下，试编写sas程序，分析年龄与晶状体混浊程度的相关关系！
                   20~                            30~                                   40~
+                215                            131                                    148
++               67                             101                                    128
+++               44                            63                                      132


data ex13_10;
do a=-1,0,1;
do b=-1,0,1;
input f @@;
output;
end; end;
cards;
215 131 148
67 101 128
44 63 132
;
proc freq;
tables a*b/chisq;
weight f;
run;
proc reg;
freq f;
model a=b;
run;
data b;
r=3; c=3;kf=74.402;
b=0.241002;sb=0.028982;
dff=(r-1)*(c-1);dfr=1;
dfb=dff-dfr;vb=sb*sb;
kr=round(b*b/vb,0.001);
kb=kf-kr;pr=1-probchi(kr,dfr);
pr=round(pr,0.0001);
if pr<0.0001 then pr=0.0001;
pb=1-probchi(kb,dfb);
pb=round(pb,0.0001);
if pb<0.0001 then pb=0.0001;
proc print;run;
我的问题：请问下面这段程序怎么来的：
r=3; c=3;kf=74.402;
b=0.241002;sb=0.028982;
dff=(r-1)*(c-1);dfr=1;
dfb=dff-dfr;vb=sb*sb;
kr=round(b*b/vb,0.001);
kb=kf-kr;pr=1-probchi(kr,dfr);
pr=round(pr,0.0001);
if pr<0.0001 then pr=0.0001;
pb=1-probchi(kb,dfb);
pb=round(pb,0.0001);
if pb<0.0001 then pb=0.0001;
proc print;run;
我搞不懂，书本上面直接给出也没有解释，求高手
帮我解释下，谢谢

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Weight tables output ROUND Input tables