@sas_user4 wrote: The following code worked just fine:   data try; set new; if cmiss(of _all_) then delete; run;   Thanks! The above code will delete a row if any of the values are blank. I thought you wanted to delete a row of all of the values are blank. ... View more

PROC GLMSELECT can't be applied to logistic regression(binomial distribution), Try PROC HPGENSELECT.   The following is 10 fold cross validation for 20 times for proc logistic. You can start with it.   %macro k_fold_cv_rep(r=1,k=10); ods select none; %do r=1 %to &r; proc surveyselect data=sashelp.heart group=&k out=have; run; %do i=1 %to &k ; data training; set have(where=(groupid ne &i)) ; run; data test; set have(where=(groupid eq &i)); run; ods output Association=native(keep=label2 nvalue2 rename=(nvalue2=native) where=(label2='c')) ScoreFitStat=true(keep=dataset freq auc rename=(auc=true)); proc logistic data=training outest=est(keep=_status_ _name_) ; class sex; model status(event='Alive')=sex height weight; score data=test fitstat; run; data score_r&r._&i; merge true native est; retain rep &r id &i; optimism=native-true; run; %end; %end; data k_fold_cv_rep; set score_r:; run; ods select all; %mend; %k_fold_cv_rep(r=20,k=10); /********************/ data all; set k_fold_cv k_fold_cv_rep indsname=indsn; length indsname $ 32; indsname=indsn; run; proc summary data=all nway; class indsname; var optimism; output out=want mean=mean lclm=lclm uclm=uclm; run; ... View more

Are you not getting "my data" for the label of the horizontal axis? The axis1 statement coupled with haxis=axis1 should show that.   Or do you mean the text of the tickmarks? "Title" is generally associated with a whole graph to describe the overall purpose, Axis Lables discribe what the individual axis measures or shows. Tickmark values would want a FORMAT associated with the variable(s) to control text.   You need to provide a bit more concrete example of what you want and it might help to show the current result you are getting.   You'll also get suggestions to move to SGPLOT where a STEP plot with x=time_months y=survival / group=sex. ... View more

To add, I got this error message (log says  P_CVSCORE in TEST statement test_proportionality is not a parameter in the model) when I used this code:   proc phreg data = ABC; class p_cvscore race1c education salary ; model asthmatt*asthma(0) = p_cvscore  tp_cvscore /risklimit; tp_cvscore=p_cvscore*log( asthmatt); test_proportionality: test p_cvscore, tp_cvscore; run;     ... View more

How many variables do you have in your model? A categorical variable with N levels counts as N-1 variables. I'm thinking you have more variables than data... ... View more

Probably use also an "ELSE" statement. DATA NEW3; SET NEW; IF SEXUALORIENT1 in (1, 3) THEN SEXORIEN1 = 10; ELSE IF SEXUALORIENT1 = 2 THEN SEXORIEN1 = 20; ELSE IF SEXUALORIENT1 in ( -7, -8, -9, 4, 5) THEN SEXORIEN1 = 8; IF SEXUALORIENT2 in (1, 2) THEN SEXORIEN2 = 10; ELSE IF SEXUALORIENT2 in (-1, -9, 3, 4, 5, 6, 7) THEN SEXORIEN2 = 8; If missing(SEXORIEN) then SEXORIEN = max(SEXORIEN1, SEXORIEN2); RUN; ... View more