About Rajesh3

Rajesh3 · ‎11-27-2019

I'm not sure about the test for the proportional odds assumption in GLIMMIX. But the error that you are getting is because your data is sparse and is not sufficient for the model you are specifying. Possible solutions would be to reduce the number of covariates or try to make your outcome binary. You can run a simple proc freq on your outcome stratified by cluster to get a general idea on how sparse your data is. I hope that helps! Thanks, Rajesh.

Rajesh3 · ‎04-16-2019

Yeah. You would need to use two different models for your hypotheses. One would be MY=DMINT (Hypothesis: DMINT has an effect on MY). I think you already have the second one. Thank you, Rajesh.

Rajesh3 · ‎04-15-2019

So you mean, you got the same error when using the code I posted? I removed event='1' from your code and added 'desc' option in the first statement. Thanks, Rajesh.

Rajesh3 · ‎04-15-2019

You can also try this code instead if your outcome is coded as a numeric variable with 0 and 1. proc logistic data = traina desc; class Gender Married N_Dependents Self_Employed Property Education (ref='1')/ param = ref; model L_Status = Credit_History N_Dependents Self_Employed L_Amount_Term Gender Married Property L_Amount Education Applicant Co_Applicant; run; Thanks, Rajesh.

Rajesh3 · ‎04-15-2019

If I am understanding you correctly, you are trying to add 'MY' next to 'DMINT' in the model statement. If this is the case, the syntax is wrong. You need to use different models for different outcomes: feed intake(DMINT), water intake(WI), milk yield(MY) and body weight (Live_wgt). The outcomes need to be continuous. In other words, you cannot simultaneously examine the effect of lunar phase / other interactions on all outcomes at once. Thanks, Rajesh.

Rajesh3 · ‎04-15-2019

The error, I suspect is for your 'Education' variable in the class statement. Thanks, Rajesh.

Rajesh3 · ‎04-08-2019

You must specify all levels of a variable. Please use variable level information from the output to make sure you are correctly understanding the reference category. Hope the below code helps! PROC GENMOD DATA=TEST; CLASS ID THERAPY (REF="0" PARAM=ref) AGE_GROUP (REF="0" PARAM=ref) SEX (REF="F PARAM=ref); MODEL COMPLETION=THERAPY AGE_GROUP SEX/ DIST=poisson LINK=log; REPEATED SUBJECT=ID/TYPE=EXCH COVB CORRW; ESTIMATE 'THERAPY 0 v. 1' THERAPY 0 1/EXP; ESTIMATE 'AGE_GROUP 0 v. 1' AGE_GROUP 0 1 0 0/EXP; ESTIMATE 'AGE_GROUP 0 v. 2' AGE_GROUP 0 0 1 0/EXP; ESTIMATE 'AGE_GROUP 0 v. 3' AGE_GROUP 0 0 0 1/EXP; ESTIMATE 'SEX female v. male' SEX 0 1/EXP; RUN; Thanks, Rajesh.

Rajesh3 · ‎03-28-2019

If you included 'dl' in the class statement, the below SAS code will apply. estimate "dl" dl 1 -1/ exp; It means for a change in day from 1 to 3/1 to 5 (this will change based on the reference chosen but the concept is this), the change (increase/decrease) in the odds of having reactivity=1. If you did not include 'dl' in the class statement, the estimate statement should be: estimate "dl" dl 1/ exp; It means, on average, for a unit increase in day, the change (increase/decrease) in the odds of having reactivity=1. Thank you, Rajesh.

Rajesh3 · ‎03-28-2019

It will not give estimates for each time point separately. Now the interpretation becomes "change in reactivity for a day increase in time". If you want to interpret it for 2 days increase multiply the estimate for 'dl' by 2 and exponentiate it to get the appropriate change. Else, replace the relevant estimate statement with the below code. estimate 'dl' dl 2/exp; Thanks, Rajesh.

Rajesh3 · ‎03-28-2019

proc genmod data=C descending order=data; class group cow po; model react=group dl po/ d=bin link=logit; repeated subject=cow/ type=cs; estimate "group" group1 -1 / exp; estimate "dl" dl/ exp; run; The 'dl' variable need not be mentioned in the class statement. If you do not mention it in the class statement, it will be treated as a continuous variable. Also, if you do not need to control for 'po' variable, please remove it from the model. Please check if the model converges with the above code. Based on the error that you are getting, you have specified a more complex model than your data can support. Hope this helps! Thank you, Rajesh.

Rajesh3 · ‎03-27-2019

proc genmod data=fake descending order=data; class treatment animal; model reactivity=treatment day/ d=bin link=logit; repeated subject=animal / type=cs; estimate "treatment" treatment 1 -1 / exp; estimate "day" day / exp; run; The above code fits a GEE model with the binomial distribution. Thank you, Rajesh.

Rajesh3 · ‎03-27-2019

Hi, If you have an add-in for SAS in excel, I think you just need minimal VBA knowledge to accomplish your task. You can execute SAS code via VBA. Once you have the developer tab activated in the excel toolbar, you can create customize a sheet as you like and add a button. However, if you want SAS to export to the same excel file at the end, you may need to close your file first so that SAS can use it. Thank you, Rajesh.

Rajesh3 · ‎03-27-2019

Hello, What you have there are crude odds ratios. Since 'insurancetype' is a categorical variable with four levels, one level should be used as a reference. Otherwise, the computation of odds would be meaningless. Hope it helps! Thank you, Rajesh.

Rajesh3 · ‎03-26-2019

Hello, If you are sure that you are not generating cartesian product, maybe you can try an alternative like using arrays or hash objects. http://support.sas.com/resources/papers/proceedings09/071-2009.pdf Thanks, Rajesh.

Rajesh3 · ‎03-26-2019

Hello, You need to be cautious when increasing the random parameters in the model. Based on the error that you got, the variance-covariance matrix might be showing negative values for variance or 0 for covariance. This means a simpler model (without a random intercept/linear trend) is a better alternative in your case. Also, I don't understand why you used 'rep' (i am assuming this to be the id variable) in the model as a covariate. Thanks, Rajesh.

Online Status	Offline
Date Last Visited	‎05-06-2020 06:58 PM

Re: Multinomial models for clustered data using GLIMMIX

Re: Proc mixed unbalanced data in animal science with multiple variabl...

Re: Class event parameter "1" not known error when running Logistic mo...

Re: Class event parameter "1" not known error when running Logistic mo...

Re: Proc mixed unbalanced data in animal science with multiple variabl...

Re: Class event parameter "1" not known error when running Logistic mo...

Re: Seeking help with the ESTIMATE function in PROC GENMOD

Re: Binomial analysis

Re: Binomial analysis

Re: Binomial analysis

Re: Proc mixed unbalanced data in animal science with multiple variabl...

Re: Proc Glimmix vs Proc Mixed

Re: How to count missing records in each variable

Re: Multinomial models for clustered data using GLIMMIX

Re: Proc mixed unbalanced data in animal science with multiple variabl...

Re: Class event parameter "1" not known error when running Logistic mo...

Re: Class event parameter "1" not known error when running Logistic mo...

Re: Proc mixed unbalanced data in animal science with multiple variabl...

Re: Class event parameter "1" not known error when running Logistic mo...

Re: Seeking help with the ESTIMATE function in PROC GENMOD

Re: Binomial analysis

Re: Binomial analysis

Re: Binomial analysis

Re: Binomial analysis

Re: I have a question about the sas add in for excel

Re: How to get crude odds ratio for binary outcome variable and catego...

Re: Proc SQL Performance

Re: Glimmix repeated measures