Hi there,   I am trying to fit a cox regression with time dependent covariate. I am pretty new to this method. I looked related paper and the SAS PHREG guide, however none of them fit my case. the measurement X was repeatedly taken and it is time dependent. I want to estimate the effect of X on the death. However, the follow-up for this study is very long, it can be up to 10 years. The unit of the time to death is days. I cannot use any array function when preparing the dataset  for the PHREG. The array is always out of range. Does any body have any recommendation? Thank you very much in advance.  Data example.  ID  X X_date death time to death 1 23 1/1/2003 1 1880 1 15 2/1/2007 1 1880 1 27 … 1 1880 1 28 … 1 1880 1 20 … 1 1880 2 … … 0 2990 2 … … 0 2990 2 … … 0 2990 2 … … 0 2990 2 … … 0 2990 2 … … 0 2990 2 … … 0 2990 2 … … 0 2990 … … … … … ... View more

Hi there,   I was trying to fit an ordinal mixed effects model using PROC GLIMMIX procedure. The code is as below. The outcome Y has 3 ordered value (1,2,3). Y_baseline is the outcome measured at baseline time. The model includes random intercept and random-effect of time, and the baseline outcome as covariate. The data includes 53 subjects with 181 observations. However the model did not converge, even if I combined the three values into 2 values and run regular logistic mixed effects mode, it still did not converge. Is it due to that the sample size is too small? Is there any solution? Thank you in advance.    PROC GLIMMIX data=final1; CLASS pt Y_baseline(ref=first); MODEL Y(desc order=internal)= Y_baseline days/DIST=multinomial LINK=cumlogit SOLUTION CL; RANDOM int days/subject=pt; RUN;   Kane ... View more

Hi there,   I have a statistical question. I am trying to run mixed effects model on repeatedly measured ordinal outcome, but the data looks not meed the proportional odds assumption. Is there any procedure can take care of this issue? Thank you very much in advance.      ... View more

Hi there, I want to display uneven value in the axis using proc sgplot, but it looks like the VALUE statement does not work. My code is attached below. But the output plot only shows 0, 12, 36, and 60 in the axis.    proc sgplot data=cif1; step x=months_to_scenario1 y=cif/group=study_id name="curve"; scatter x=months_to_scenario1 y=cenp / markerattrs=(symbol=plus) GROUP=study_id ; text x=t1 y=s1 text=texc1/position=center textattrs=(size=10) outline; xaxis offsetmin=0.015 offsetmax=0.015 max=60 values=(0 1 12 24 36 48 60); yaxis max=0.25; keylegend "curve"; run;   Thank you.  ... View more

Hi there,   I already download and install proc traj procedure but when I run it, I got the error message below...Does any one have idea how to solve it? Thanks.  ERROR: An exception has been encountered. Please contact technical support and provide them with the following traceback information: The SAS task name is [TRAJ] ERROR: Read Access Violation TRAJ Exception occurred at (2E80279C) Task Traceback Address Frame (DBGHELP API Version 4.0 rev 5) 000000002E80279C 000000002EA8D400 traj:mcn_main+0x1179C 000000002E817637 000000002EA8D5B0 traj:mcn_main+0x26637 000000002E810330 000000002EA8D650 traj:mcn_main+0x1F330 000000002E7F1C6B 000000002EA8FBF0 traj:mcn_main+0xC6B 00000000033A89DB 000000002EA8FBF8 sashost:Main+0x10EBB 00000000033AE61D 000000002EA8FF50 sashost:Main+0x16AFD 00007FFCE63B3034 000000002EA8FF58 KERNEL32:BaseThreadInitThunk+0x14 00007FFCE7FF1471 000000002EA8FF88 ntdll:RtlUserThreadStart+0x21 ... View more