不知道有没有表述清楚,如下为整段过程,目的是练习求卡方值和P值.

data _null_;
array x {1:2,1:3} x11 x12 x13 x21 x22 x23(934 113 20 1070 92 8);
array exp{1:2,1:3} e11 e12 e13 e21 e22 e23;
array s{1:2} s1_ s2_;
array s_{1:3} s_1 s_2 s_3;
sum=sum(of x{*});
do i=1 to 2;
s=sum(of x{i, });end;
do i=1 to 3;
s_=sum(of x{ ,i});
end;
do i=1 to 2 ;
do j=1 to 3;
exp[i,j]=s*s_[j]/sum;
chisq +(x[i,j]-exp[i,j])**2/exp[i,j];
end;
end;
put chisq=;
p=cdf('chisq',chisq,3);
put p=;
put _all_;
run;