Re: GENLINMIXED EMMEANS don't agree with descriptives

Posted by Bruce Weaver on Feb 27, 2021; 1:44am
URL: http://spssx-discussion.165.s1.nabble.com/GENLINMIXED-EMMEANS-don-t-agree-with-descriptives-tp5740214p5740225.html

Thanks for providing that example.  I am not surprised that the EMMEANS do
not match the MEANS, given that you are using a binomial error distribution
and a logit link.  If you change those to normal error distribution and
identity link (as shown below), the EMMEANS do match the MEANS.  (If this
"risk difference" kind of model makes sense, you'll probably want to use the
ROBUST covariance estimator too, as I have done below.)  


VARIABLE LEVEL CHOICE (SCALE).
GENLINMIXED /DATA_STRUCTURE SUBJECTS=SUB
 /FIELDS TARGET=CHOICE TRIALS=NONE OFFSET=NONE
 /TARGET_OPTIONS DISTRIBUTION=NORMAL LINK=IDENTITY REFERENCE=0
 /FIXED EFFECTS= AC WT AC*WT USE_INTERCEPT=TRUE
 /RANDOM USE_INTERCEPT=TRUE SUBJECTS=SUB
 /EMMEANS TABLES=AC*WT
 /BUILD_OPTIONS TARGET_CATEGORY_ORDER=aSCENDING
INPUTS_CATEGORY_ORDER=aSCENDING COVB = ROBUST.

MEANS CHOICE BY AC BY WT.



Beckstead, Jason-2 wrote

> Thanks Bruce. Below is small data and syntax.
>
>
> * Encoding: UTF-8.
> *-9 within-subjects conditions (3 x 3; AC x WT) & binary choice response.
> *-below is hypothetical data from N=6 subjects.
> *-At issue: estimated means do not agree with observed means.
> DATA LIST FREE /SUB AC WT CHOICE.
> BEGIN DATA.
> 1 1 1 0
> 1 1 2 0
> 1 1 3 0
> 1 2 1 0
> 1 2 2 1
> 1 2 3 1
> 1 3 1 1
> 1 3 2 1
> 1 3 3 1
> 2 1 1 0
> 2 1 2 0
> 2 1 3 1
> 2 2 1 0
> 2 2 2 0
> 2 2 3 1
> 2 3 1 0
> 2 3 2 1
> 2 3 3 1
> 3 1 1 0
> 3 1 2 1
> 3 1 3 0
> 3 2 1 1
> 3 2 2 1
> 3 2 3 0
> 3 3 1 1
> 3 3 2 1
> 3 3 3 0
> 4 1 1 0
> 4 1 2 0
> 4 1 3 0
> 4 2 1 1
> 4 2 2 1
> 4 2 3 1
> 4 3 1 1
> 4 3 2 1
> 4 3 3 1
> 5 1 1 0
> 5 1 2 0
> 5 1 3 0
> 5 2 1 0
> 5 2 2 0
> 5 2 3 0
> 5 3 1 0
> 5 3 2 0
> 5 3 3 0
> 6 1 1 1
> 6 1 2 1
> 6 1 3 1
> 6 2 1 1
> 6 2 2 1
> 6 2 3 1
> 6 3 1 1
> 6 3 2 1
> 6 3 3 1
> END DATA.
>
> VARIABLE LEVEL SUB AC WT (NOMINAL).
> GENLINMIXED /DATA_STRUCTURE SUBJECTS=SUB
>  /FIELDS TARGET=CHOICE TRIALS=NONE OFFSET=NONE
>  /TARGET_OPTIONS DISTRIBUTION=BINOMIAL LINK=LOGIT REFERENCE=0
>  /FIXED EFFECTS= AC WT AC*WT USE_INTERCEPT=TRUE
>  /RANDOM USE_INTERCEPT=TRUE SUBJECTS=SUB
>  /EMMEANS TABLES=AC*WT
>  /BUILD_OPTIONS TARGET_CATEGORY_ORDER=aSCENDING
> INPUTS_CATEGORY_ORDER=aSCENDING.
>
> MEANS CHOICE BY AC BY WT.

-----
--
Bruce Weaver
[hidden email]
http://sites.google.com/a/lakeheadu.ca/bweaver/

"When all else fails, RTFM."

NOTE: My Hotmail account is not monitored regularly.
To send me an e-mail, please use the address shown above.

--
Sent from: http://spssx-discussion.1045642.n5.nabble.com/

=====================
To manage your subscription to SPSSX-L, send a message to
[hidden email] (not to SPSSX-L), with no body text except the
command. To leave the list, send the command
SIGNOFF SPSSX-L
For a list of commands to manage subscriptions, send the command
INFO REFCARD