Hi I am working on a regression model with repeated measures, and my response variable is continuous. Each pair of physician-patient-drug is unique, but the physician may have several patients who are prescribed different medications. My dataset looks like this (My dataset contains approximately 15000 rows, approximately 15,000 patients are treated by around 150 physicians.) data drugs; input physician patient drug $ outcome age gender$; datalines; 1 1 a1 83.1 55 F 1 1 b5 76.2 55 F 1 1 b7 99.9 55 F 1 4 f2 2.7 82 M 2 4 a1 40.2 82 M 2 5 a2 88.1 65 M 3 6 c8 26.1 49 F 4 2 a1 23.3 89 F 5 3 b1 92.1 67 M 5 3 b2 92.4 67 M 6 7 e4 12.1 28 M 6 8 a2 5.1 73 M 6 8 b2 3.1 73 M 7 9 a2 98.9 79 M ;   proc mixed data=drugs noclprint; class patient physician drug gender(ref="F") ; model outcome= gender age /solution outp=preddata ; repeated DRUG /subject=physician(patient) type=cs ; RUN; I am working on evaluating whether certain drugs are prescribed correctly, which will be my outcome measure. Do you believe the model is accurate? This is my first time working with this type of data. Thank you for your response. Teresa ... View more

Hi I am working on an alternative logistic regression model with repeated measures The drug is the unit of analysis, and there are TWO levels of clustering—within patient, and patient within MD My dataset looks like this (My dataset contains approximately 15000 rows) data drugs; input patient MD Gender Age drug $ indication outcome; datalines; 1 1 2 54 a 1 1 1 1 2 54 b 1 1 1 1 2 54 c 0 0 2 4 1 41 a 1 0 2 4 1 41 c 1 0 2 4 1 41 e 1 1 3 1 1 24 h 0 0 4 5 2 29 c 1 0 5 1 2 72 a 1 1 5 1 2 72 a 0 1 6 2 2 72 i 0 1 6 2 2 72 b 0 0 7 1 1 36 a 0 0 8 3 1 25 a 0 1 ; PROC GENMOD DATA= drugs; CLASS MD PATIENT gender(ref="1") indication(ref="1"); MODEL outcome (event='1')= gender indication age / DIST=BIN ; REPEATED SUBJECT= PATIENT / TYPE=EXCH ; RUN ;   What is wrong with my model? Thank you very much for your help ... View more

I need to create a ROC curve. My data looks like : data Data1; input id score diagnosis; datalines; 1 0.1 0 1 0.2 0 1 1.0 1 1 0.8 0 1 0.9 0 2 1.8 1 2 1.2 1 3 1.7 1 3 1.1 0 4 0.1 0 4 1.9 0 4 1.8 1 4 1.2 1 4 2.0 0 5 0.7 0 5 1.2 1 ; I have several tests per ID. My sas code following ods graphics on; proc logistic data=data1 plots(only)=(roc(id=obs) effect) descending; model diagnosis(event='1')=score / ctable outroc=RocStats ; ROC; ROCCONTRAST; run; ods graphics off; My question: I don't know how to integrate the repeated ID into my model. I use only diagnosis and score. O am confused about this issue Thanks: ... View more