请高手帮忙解决一下

0关注
1粉丝

大专生

78%

还不是VIP/贵宾

-

0%

威望: 0 级
论坛币: 58 个
通用积分: 0
学术水平: 0 点
热心指数: 0 点
信用等级: 0 点
经验: 2106 点
帖子: 23
精华: 0
在线时间: 91 小时
注册时间: 2009-11-10
最后登录: 2016-6-4

楼主

yan_0103 发表于 2011-6-15 11:48:48 |AI写论文

是否 +2 论坛币

k人参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

proc iml;
use _cox_;
read all into coxm;
n=nrow(coxm);
col=ncol(coxm);
p=col-3;
group=coxm[,1];
t=coxm[,2];
cens=coxm[,3];
z=t(coxm[ ,4:col]);
free coxm;
gvalue=unique(group)`;
gb=nrow(gvalue);
gm=j(gb,1,0);
gd=j(gb,1,0);
gindex=j(n,1,0);
do i=1 to gb;
  do j=1 to n;
if group[j]=gvalue[i] then
   do;
      gindex[j]=i;
      gm[i]=gm[i]+1;
      if cens[j]=1 then gd[i]=gd[i]+1;
   end;
  end;
end;

print '~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~';
print 'Data  Profile';
print '~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~';
reset nocenter;
if %index(a,.)=0 then  %let dsname=work.a;
else %let dsname=a;
print " Data set: %upcase(&dsname)";
print " Grouping Variable: group";
print " Dependent Variable: t";
print " Censoring Variable: cens";
print " Censoring Value: 0";
print " Ties Handling: BRESLOW";
reset center;
label='';
mattrib n       colname='Total'    label=label
      nevent colname='Event'    label=label
      ncens colname='Censored'  label=label
      percentc colname='%Censored' label=label format=6.2;
nevent=sum(cens);
ncens=n-nevent;
percentc=100*ncens/n;
print ,;
print ' Summary of the Number of ';
print ' Event and Censored Values';
print n nevent ncens percentc;
print ,;
call symput('mac_gb',trim(left(char(gb,12,0))));
reset nocenter;
print " Total number of groups: &mac_gb";
reset center;

if %upcase(greport)=Y then %str(
print /;
print '~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~';
print 'Summary of Group Size and the Number of Events Per Group';
print '~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~';
gseq=(1:gb)`;
mattrib gseq colname='Sequence#' label=label
         gvalue colname='Group'    label=label
         gm    colname='Size'    label=label
         gd    colname='Events' label=label;
Print gseq gvalue gm gd;
);
quit;
在日志中出现了(execution) Matrix has not been set to a value.
不知道哪块出问题了？还请大家帮忙看一下，不胜感激