Greetings all. I fit an logistic regression model which the variables contains treatment(binary) and some explanatory variable(continuous). Now I want to edit the treatment estimate from the same model to calculate the scores for the purpose of obtaining the ROC analysis. I saved model information to a data set with the OUTMODEL= option, and I'm trying to build new data set("model2") which replace the estimate and the standard error from previous model information ("mode11"), then use it to recalculate the scores. But the following error is issued and the step terminates: " ERROR: Computations are terminated because the INMODEL=WORK.model2 data set has lost information. " My code is: /*logistic model*/ proc logistic data=data1 descending outmodel=model1; class y t(ref="0")/param = ref; model y(event='1')=t c1 c2 u; run; /*output predict probability*/ proc logistic inmodel=model1; score data=data1 out=probpred1; format P_1 10.7; run; /*new model (have error?)*/ data model2(type=logismod); set model1; if (_TYPE_='E') AND (_NAME_='EFFECT') AND (_CATEGORY_='E') AND (_NAMEIDX_=0) AND (_CATIDX_=0) then _MISC_=0.7802293905; /*replace treatment parameter estimate*/ if (_TYPE_='E') AND (_NAME_='EFFECT') AND (_CATEGORY_='V') AND (_NAMEIDX_=.) AND (_CATIDX_=0) then _MISC_=0.123159491; /*replace standard error*/ run; /*recalculate score (have error!!)*/ proc logistic inmodel=mode12(type=logismod); score data = data1 out=probpred2; run; Any help would be appreciated! Thanks!
... View more