BookmarkSubscribeRSS Feed
Joanna7726
Calcite | Level 5

Hello, 

 

I am using the below syntax to retreivew estiamte values but as you can see in the output, the syntax is only estimating the last three commands and ignoring the rest of the syntax. Could anyone let me know what might be wrong with the syntax or why I am not getting outputs for the other values? Thank you very much for your help!

 

*0 = no stressor ; *1 = stressor;

proc mixed noclprint covtest;
 class m2id stressyn racenew gender chronic workyn marital;
 model NA =  gender C_edu racenew C_age Z_N_hhinc chronic C_N_neuro workyn marital C_PA m_stress C_lifedis|stressyn C_daydis|stressyn
/solution CL ddfm=bw;
   random intercept stressyn /subject=m2id;

estimate 'Estimate of outcome variable for non-stressor days in -1SD Discrimination' intercept 1 stressyn 0 C_daydis*stressyn -4.7007771;
estimate 'Estimate of outcome variable for non-stressor days in Mean Discrimination' intercept 1 stressyn 0 C_daydis*stressyn 1;
estimate 'Estimate of outcome variable for non-stressor days in +1SD Discrimination' intercept 1 stressyn 0 C_daydis*stressyn 4.7007771;
estimate 'Estimate of outcome variable for stressor days in -1SD Discrimination'     intercept 1 stressyn 1 C_daydis*stressyn -4.7007771;
estimate 'Estimate of outcome variable for stressor days in Mean Discrimination'     intercept 1 stressyn 1 C_daydis*stressyn 1;
estimate 'Estimate of outcome variable for stressor days in +1SD Discrimination'     intercept 1 stressyn 1 C_daydis*stressyn 4.7007771;

estimate '-1SD Discrimination Effect for Non-Stressor Days' C_daydis -4.7007771 C_daydis*stressyn 0/cl;
estimate 'Mean Discrimination Effect for Non-Stressor Days' C_daydis 1 C_daydis*stressyn 0/cl;
estimate '+1SD Discrimination Effect for Non-Stressor Day'  C_daydis 4.7007771 C_daydis*stressyn 0/cl;
estimate '-1SD Discrimination Effect for Stressor Day'      C_daydis -4.7007771 C_daydis*stressyn -4.7007771/cl;
estimate 'Mean Discrimination Effect for Stressor Day'      C_daydis 1 C_daydis*stressyn 1/cl;
estimate '+1SD Discrimination Effect for Stressor Day'      C_daydis 4.7007771  C_daydis*stressyn 4.7007771/cl;
run;

Joanna7726_0-1652758408090.png

1 REPLY 1
PaigeMiller
Diamond | Level 26

Non-est means you can't estimate this from the model fit. Sometimes this happens because you have interactions in the model, but there are empty cells i.e. combination of variables in an interaction such as C_lifedis*stressyn (or any other interaction in the model) has at least one empty cell.

 

Also, to get predicted values and confidence intervals at a particular set of x-values is via these methods. The ESTIMATE statement is not the way to do it.

https://blogs.sas.com/content/iml/2014/02/17/the-missing-value-trick-for-scoring-a-regression-model....

https://blogs.sas.com/content/iml/2014/02/19/scoring-a-regression-model-in-sas.html

 

--
Paige Miller

SAS Innovate 2025: Call for Content

Are you ready for the spotlight? We're accepting content ideas for SAS Innovate 2025 to be held May 6-9 in Orlando, FL. The call is open until September 25. Read more here about why you should contribute and what is in it for you!

Submit your idea!

What is ANOVA?

ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.

Find more tutorials on the SAS Users YouTube channel.

Discussion stats
  • 1 reply
  • 299 views
  • 0 likes
  • 2 in conversation