sasrules Tracker
https://communities.sas.com/kntur85557/tracker
sasrules TrackerFri, 19 Apr 2024 16:22:33 GMT2024-04-19T16:22:33ZRe: PROC GLIMMIX: Different p-values for different reference groups?
https://communities.sas.com/t5/Statistical-Procedures/PROC-GLIMMIX-Different-p-values-for-different-reference-groups/m-p/264432#M13949
<P>Thank you!</P>Sun, 17 Apr 2016 17:40:20 GMThttps://communities.sas.com/t5/Statistical-Procedures/PROC-GLIMMIX-Different-p-values-for-different-reference-groups/m-p/264432#M13949sasrules2016-04-17T17:40:20ZPROC GLIMMIX: Different p-values for different reference groups?
https://communities.sas.com/t5/Statistical-Procedures/PROC-GLIMMIX-Different-p-values-for-different-reference-groups/m-p/263092#M13886
<P>Hello All,</P>
<P> </P>
<P>I'm running a multinomial logistic model with a random effect using PROC GLIMMIX. I've noticed that I get different p-values for my independent variable depending on which reference group I use. This seems strange to me and doesn't happen if I remove the random effect. I am wondering if I might be doing something incorrectly with the random effect. In particular, I don't understand why I have to add the "group" option for the random effect, but SAS seems to require that I do so! Please see my code below. Thank you!</P>
<PRE><CODE class=" language-sas">proc glimmix data=lib.submixed;
class author placement;
model placement(ref='Middle')=year/dist=mult link=glogit; /*p=0.0026 for year*/
random int / subject=author group=placement;
run;
proc glimmix data=lib.submixed;
class author placement;
model placement(ref='End')=year/dist=mult link=glogit; /*p=0.0089 for year*/
random int / subject=author group=placement;
run;<BR /></CODE></PRE>Tue, 12 Apr 2016 04:14:46 GMThttps://communities.sas.com/t5/Statistical-Procedures/PROC-GLIMMIX-Different-p-values-for-different-reference-groups/m-p/263092#M13886sasrules2016-04-12T04:14:46ZRe: Predicted probabilities for mixed multinomial model
https://communities.sas.com/t5/Statistical-Procedures/Predicted-probabilities-for-mixed-multinomial-model/m-p/263085#M13885
<P>Thank you so much for your help, Steve! I had pretty much given up on being able to do this with SAS, but I'll take a look at the example you suggested.</P>Tue, 12 Apr 2016 03:32:06 GMThttps://communities.sas.com/t5/Statistical-Procedures/Predicted-probabilities-for-mixed-multinomial-model/m-p/263085#M13885sasrules2016-04-12T03:32:06ZPredicted probabilities for mixed multinomial model
https://communities.sas.com/t5/Statistical-Procedures/Predicted-probabilities-for-mixed-multinomial-model/m-p/254780#M13448
<P>I'm trying to create a mixed model with a multinomial outcome and a random effect. I want to model the probability of a speech tag in a sentence coming at the beginning, middle, or end of a sentence, and I have a large sample of sentences. The outcome can have three levels (beginning, middle, end), the only fixed effect is decade of publication, and I want to account for author as a random effect, since there are multiple sentences in the sample written by the same author.</P>
<P> </P>
<P>I want to graph the predicted probability of each type of speech tag (beginning, middle, or end) by decade of publication. I figured I could use LSMEANS to get the predicted probability of each outcome (beginning, middle, end). But, I can't seem to figure out how to do this using SAS.</P>
<P> </P>
<P>I've tried to use PROC LOGISTIC, PROC GENMOD, and PROC GLIMMIX. PROC LOGISTIC won't take the random effect into account, but will give me the correct LSMEANS if I use link=glogit. PROC GENMOD will account for the random effect but will not allow link=glogit, and I am not sure how to interpret the LSMEANS when link=cumlogit because the probabilities of the three outcomes do not add up to one. PROC GLIMMIX won't allow LSMEANS for the multinomial distribution. Attached are the codes I've tried. Does anyone have any ideas?</P>
<P> </P>
<P> </P><BR /><IMG src="https://communities.sas.com/t5/image/serverpage/image-id/12354i2A9A5F5599AC8C94/image-size/large?v=1.0&px=600" border="0" alt="Screen Shot 2016-03-05 at 4.12.45 PM.png" title="Screen Shot 2016-03-05 at 4.12.45 PM.png" />Sun, 06 Mar 2016 02:43:03 GMThttps://communities.sas.com/t5/Statistical-Procedures/Predicted-probabilities-for-mixed-multinomial-model/m-p/254780#M13448sasrules2016-03-06T02:43:03ZPrint multiple survival plots on the same page with PROC LIFETEST
https://communities.sas.com/t5/SAS-Procedures/Print-multiple-survival-plots-on-the-same-page-with-PROC/m-p/230613#M54505
<P>Hi All,</P><P> </P><P>I am trying to print 14 survival plots on the same page using PROC LIFETEST. I could use another procedure if anyone has any suggestions. Please see the example code below. I would like a plot for each level of "group" (group has 14 levels). But, I want them all on the same page rather than one per page. There seem to be many examples of how to do this using PROC GPLOT, but not for a procedure like PROC LIFETEST.</P><P> </P><P>Thanks in advance for any help!</P><P>Anna</P><P> </P><P>proc lifetest data=testdata plots=survival;<BR />by group;<BR />time surtim * surind(0);<BR />strata ethnicity;<BR />run;</P>Mon, 19 Oct 2015 17:25:48 GMThttps://communities.sas.com/t5/SAS-Procedures/Print-multiple-survival-plots-on-the-same-page-with-PROC/m-p/230613#M54505sasrules2015-10-19T17:25:48Z