Hello,   I would like to ask whether the below graphs are all different for some reasons and whether I can fix them out or I can have supporting idea to explain them.  Three of them come from exactly the same data but with different proc procedures.  Major concern is that proc plm gave me a probability with CI, which I wanted to have, but that is different from logistic with spline effect also different from what I have from excel graph. I have the blue graph to estimate the probability value based on the equation, which looks similar to adaptivereg and logistic spline. Whilst the plm graph are linear, why the rest of them have u shape? Is this because of the logic behind each procedures are different of I have done in wrong way round?  Another question is, I guess logistic spline and the excel graph are supposed to be  the same as the both have u shape. But then logistic spline has its lowest probability at age 26-27 whilst the blue graph, which I have made via excel has 30-32 at its lowest.. Does it look like excel seems to be wrong? But then I have checked many times but havent found what was wrong.  But anyway, please let me know whether it seems the blue one needs to be changes since logistic with spline effect are right ... View more

Hello,  I am carrying on a project that is about to have predicted probability and would like to ask you several questions. 1. To see whether there is relationship between age and health outcome, I did logistic with spline effect as followed then have sgplot. PROC LOGISTIC DATA=HI PLOTS=NONE;   EFFECT SPL= SPLINE(AGE95/DETAILS NATURALCUBIC BASIS=TPF(NOINT) KNOTMETHOD=PERCENTILES(5) );   MODEL HEALTH (EVENT=‘1’)= SPL;   OUTPUT OUT=SPLINE_HEALTH PREDICTED =PREDPROB; QUIT; PROC SGPLOT DATA= OUT=SPLINE_HEALTH NOAUTOLEGEND;   SCATTER SCATTER X=AGE95 Y=HEALTH(EVENT=‘1’) / JITTER   SERIES SCATTER X=AGE95 Y=PREDPROB;    YAXIS MIN=0 MAX=1 GRID LABEL=“PREDICTED PROBABILITY”; RUN; Then I also want to have CI on the graph using band statement but it does not give me what I wanted for (e. g. attach 1 from SAS/STAT® 14.1 User’s Guide The GAMPL Procedure p. 2918). Can you recommend what I need do for CI on a graph?    2. The other method I am doing is non-parametic with proc adaptivereg  PROC ADAPTIVEREG DATA=HI_2; MODEL HEALTH (EVENT=‘1’)=AGE95; OUTPUT OUT METABOLIC P(ILINK); RUN; PROC SGPLOT DATA=METABOLIC NOAUTOLEGEND; SCATTER X=AGE95 Y=HEALTH / JITTER; SERIES X=AGE95 Y=PRED; YAXIS MIN=0 MAX=1 GRID LABEL=“PREDICTED PROBABILITY”; RUN; Here, again I would like to have CI on the graph but it has been challenging.  Is it because Ci cannot be calculated or is it me not to have CI?  Please let me know how to deal with this issue.    Many thanks in advance! ... View more

Hello,  I would like to ask questions related to my project that I am currently doing.  Have a data which is nonparametric, which fairly recently relaised..  So before I notice I need to apply nonparametric method, I did logistic regression (univariate, univariable).  Then as I said data has characteristics and now studying nonparametric method but there are many things that I am no sure.  So my questions are: 1. Can I use "proc adaptivereg" when there is one variable and one variate?     This is because I have looked at below and kindly explained that the proc adaptivereg is for multivariate.. http://support.sas.com/documentation/cdl/en/statug/68162/HTML/default/viewer.htm#statug_adaptivereg_examples03.htm 1-1. when it is multivariate, does that mean I can have only one variable or it doesnt matter in multivariate setting?    2. If question 1 is no, then can you recommend syntax for nonparametric data for simple logistic regression?  (I used proc logistic before)     It would be great help for me.    Many thanks and best wishes,   Danielle.  ... View more