Thank you ksharp. I tried repeated year/ subject=id type=chol; and I don't think chol is an option for the variance covariance structure. I got the following error:   ERROR 22-322: Syntax error, expecting one of the following: ANTE, AR, ARH, ARMA, CS, CSH, FA, FA0, FA1, HF, LIN, LINEAR, SIMPLE, SP, TOEP, TOEPH, UN, UNAR, UNCS, UNR, VC.     any additional ideas? Is it possible for me to try proc glm with REML?   Thanks. ... View more

Hi,   I am trying to run a proc mixed analysis for the following data:   ID            Score        year         group        othervariables 1                80             0              1 1                90             1              1 1                77             2              1 2                66             0              0 2                90             3              0 2                82             6              0 .                  . .                  . .                  .   So the years with scores are from 0 to 7 but some have scores years 0, 1,2, 3, 6 and some 0, 3, 6 etc. It varies. The dataset is fairly large and includes 17155 observations and 223 variables...I need to adjust for some variables in the model (around 7 etc).   When running the unadjusted model using the following code: proc mixed data=cup.cup; class year group; model F393MSE = group year group*year; repeated year/ subject=id type=ar(1); run;   I get the following error: WARNING: ODS graphics with more than 5000 points have been suppressed. Use the PLOTS(MAXPOINTS= ) option in the PROC MIXED statement to change or override the cutoff. WARNING: Stopped because of infinite likelihood.   I added the  PLOTS(MAXPOINTS= NONE) in the statement, but it didn't work. Surprisingly, when I added covariates in the model (around seven) the model runs. However, I do need to run the unadjusted.   Also, one other question. How does one run a GLM with REML method and specify the random effects...   Thank you for your help! ... View more

No, I haven't done that part as I assumed it is already built in in SAS. Excuse my ignorance but what is the best way of downloading the template?   Also, I have been reading the following instrucyions: https://support.sas.com/documentation/cdl/en/statug/63962/HTML/default/viewer.htm#statug_templt_a0000000061.htm    How does one use the layout overlay statement with the proc lifetest..   I really think this should be easier, as in other programs like STATA it is much easier...   any ideas?     ... View more

This is the error from the log:   12018 %CompileSurvivalMacros; - 180 WARNING: Apparent invocation of macro COMPILESURVIVALMACROS not resolved. ERROR 180-322: Statement is not valid or it is used out of proper order.   Many thanks!! ... View more

tried that. It didn't work...   %ProvideSurvivalMacros; %let xOptions = label="Time to Event (months)" linearopts=(viewmin=0 viewmax=12) tickvaluelist=(0 2 4 6 8 10 12)); %CompileSurvivalMacros; proc lifetest data=Work.Work atrisk plots=survival(atrisk cb) outs=out; strata type; time visit*status(0); run;   am I missing something? ... View more

Thanks! no worries. you just answered my question. I was wondering if there was another way of getting the median etc and variance directly....     ... View more

thanks, yes I know proc means gives you the summary statistics and you can get the mean, median and variance. but I want to generate a whole new dataset showcasing the mean, median, and variance for each feature by ID...   ... View more

thank you. Sorry, I thought I updated the group. It worked out in the end.I agree that there were some unnecessary steps in my code.       ... View more

Thanks.. sure. No errors in the log...The results is just a status variable with all 0 inputs..   So this how I am trying to make the data look for the status variable:   ID         TXdate                     Time           Status 1                                              0                  0 1                                              1                  0 1                                              2                  0 1                                              3                  0 1                                              4                  0 2                                              0                  0 2                                              1                  0 2             MAR122010             2                  1 2                                              3                  1 2                                              4                  1 3                                              0                  0 3             MAY092012              1                  1 3                                              2                  1  3                                              3                  1   Thanks   ... View more

Thank you. I tried the following, but it doesn't seem to work. any ideas?    data want (drop=get) set have; by ID; retain get;   if first.id then get=0; status=0; if get eq 0 and TXdate^= ' ' then do; status=1; get=1; end; run; ... View more

Thank you...However that didn't work   I am trying the following:   data want (drop=get) set have; by ID; retain get;   if first.id then get=0; status=0; if get eq 0 and TXdate^= ' ' then do; status=1; get=1; end; run;   This also isn't working. Any ideas, of what I could be doing wrong   ... View more

Thank you!   Quick question...Can I assign a 1 to the status variable as soon as it hits a Txdate?   This time instead of an ID being 1 if there is TXdate. The status variable only turns to 1 once it hits a Txdate at a certain time point...   Can I make the data look like this:   ID          TXdate              Time         status 1                                        0               0 1                                        1               0 1           MAR102016         2               1 1                                        3               1 1                                        4               1     ... View more