I used the following SAS code to read in all .csv file within a work directory. If I need to read in line 3 to line 153 of each .csv file, how should I modify the code? Thanks much for the help!   %macro drive(dir,ext); %local cnt filrf rc did memcnt name; %let cnt=0; %let filrf=mydir; %let rc=%sysfunc(filename(filrf,&dir)); %let did=%sysfunc(dopen(&filrf)); %if &did ne 0 %then %do; %let memcnt=%sysfunc(dnum(&did)); %do i=1 %to &memcnt; %let name=%qscan(%qsysfunc(dread(&did,&i)),-1,.); %if %qupcase(%qsysfunc(dread(&did,&i))) ne %qupcase(&name) %then %do; %if %superq(ext) = %superq(name) %then %do; %let cnt=%eval(&cnt+1); %put %qsysfunc(dread(&did,&i)); proc import datafile="&dir\%qsysfunc(dread(&did,&i))" out=dsn&cnt dbms=csv replace; run; %end; %end; %end; %end; %else %put &dir cannot be opened.; %let rc=%sysfunc(dclose(&did)); %mend drive; %drive(c:\temp,csv)     ... View more

I am using PROC GLIMMIX to do an analyses, My code is below   PROC GLIMMIX DATA=mfirst METHOD=LAPLACE NOCLPRINT; CLASS SCHOOLNUM; MODEL PASS (EVENT="1")=peat timeaftergrad/CL DIST=BINARY LINK=LOGIT SOLUTION; RANDOM INTERCEPT peat timeaftergrad/ SUBJECT=SCHOOLNUM S CL TYPE=VC; /*parms (0.2) (1) / hold=2;*/ COVTEST /WALD; ods output solutionR=mfirst_rest_mod2; run; I got the following error message   "Obtaining minimum variance quadratic unbiased estimates as starting values for the covariance parameters failed"   Please help. Thanks!   ... View more

I am using logistic regression for a project, I want to examine whether prediction residuals differ between groups. From SAS output, it provides both Pearson residual and Deviance residual. If I want to examine whether prediction residual differ between groups, which one should I use, Pearson residual or Deviance residual? As SAS code below, reschi=pr  is Pearson residual,  resdev=dr  is Deviance residual. I spent some time searching, but I haven't got clear idea yet, so I come here to ask the experts. thanks!    proc logistic data=Have; where first=1; class aaa ppp; model pass (event="1")=aaa ppp timeaftergrad schoolpass1; output out=out_first p=prob xbeta=logit resdev=dr h=pii reschi=pr difchisq=difchi; run;   ... View more

I am using logistic regression for a project, I want to examine whether prediction residuals differ between groups. From SAS output, it provides both Pearson residual and Deviance residual. If I want to examine whether prediction residual differ between groups, which one should I use, Pearson residual or Deviance residual? As SAS code below, reschi=pr  is Pearson residual, resdev=dr  is Deviance residual. I spent some time searching, but I haven't got clear idea yet, so I come here to ask the experts. thanks!    proc logistic data=Have; where first=1; class aaa ppp; model pass (event="1")=aaa ppp timeaftergrad schoolpass1; output out=out_first p=prob xbeta=logit resdev=dr h=pii reschi=pr difchisq=difchi; run; ... View more