Re: Interactions binary logistic regression

Posted by Rich Ulrich on Jan 23, 2017; 5:56pm
URL: http://spssx-discussion.165.s1.nabble.com/Interactions-binary-logistic-regression-tp5733602p5733734.html

I think that I would have to modify this procedure, a couple of ways.

    "The 3-way interaction was
    significant, thus I did not interpret main effects or 2-way interaction and
    I focused only on the significant 3-way interaction. In order to explain
    this significant 3-way interaction, I ran pairwise comparison corrected with
    the bonferonni procedure and I described comparisons that answered my
    hypotheses."

1) How big are the various effects?  If the N is big enough, it is easy for trivial effects to

be magnified to "significance" at the arbitrary level.  "Be cautious" about interpreting

effects is a more cautious way of giving the conventional warning.  You may hear

"do not interpret" from people who are looking at your analysis if they see that the high-

order interaction is about as big as the others, but they are skipping the statement of the

assumption, that these are all the same size. (Professors who have not done much data

analysis can also be guilty of this oversight.)

2) Before worrying about followup tests, look at the means.  Is there an obvious story

there?  Do you have "basement" or "ceiling" effects which would introduce interactions?

 - those are the sort that are easiest to ignore, as potential artifacts of scaling.  (Dichotomies

created from continuous scales introduce additional potential for artifact, but, like analyzing

ranks, these are apt to be too small to show up unless the N is very large.)

Paired tests are the right sort to run for follow-up, but I don't see that they should deliver

a good "answer" for your hypotheses, because interactions disrupt the narrative of the main

effects.  I think it is awkward if your design puts the hypothesis into the interaction.  What an

unexpected and strong interaction sometimes shows you is that your whole model should be

constructed differently.  -- I learned that from a Usenet newsgroup poster, 15 or 20 years ago.

--

Rich Ulrich

---

From: SPSSX(r) Discussion <[hidden email]> on behalf of PhD student <[hidden email]>
Sent: Monday, January 23, 2017 4:47:05 AM
To: [hidden email]
Subject: Re: Interactions binary logistic regression

 

After the reading of this very helpful book (Heck, Thomas & Tabata, 2010), it
would appear that GEE are particularly appropriate for my data. However I
still have a little question in how to explain my significant 3 way
interaction.

Just as a reminder, my participants (Anxiety Disorders vs. Control) carried
out 8 trials resulting from the combination of a 2 (Formulation : Gain vs
Loss) and a 4 (Level of risk : 20% vs 40% vs 60% vs 80%) within subject
factors. In each trial, they had to choose between a safe and a risky
option.

As I know that it is easy to make mistakes with these “recent” analyses (I
did not find any well-written example in my field and any of my
supervisors/colleagues know much about this analysis), I will briefly
explain how I proceed.

First, I chose a logit link and a binomial error distribution. Then, in
order to select the most appropriate working correlation structure, I ran 3
different models using either the independent, autoregressive or
exchangeable correlation structure. I found that QIC for the independent,
autoregressive and exchangeable were in the same range (i.e., vary less than
2 points between these 3 structures). As complete independence between
trials within the same individual is not credible and as any of my
participants report changing his strategy across trials, I selected an
exchangeable working correlation structure.

Second, I ran a full factorial model (including 3 main effects: Group,
Formulation and Level of Risk; three 2-way interaction: Group X Formulation,
Group X Level of Risk and Formulation X Level of Risk; one 3 way
interaction: Group X Formulation X Level of Risk). The 3-way interaction was
significant, thus I did not interpret main effects or 2-way interaction and
I focused only on the significant 3-way interaction. In order to explain
this significant 3-way interaction, I ran pairwise comparison corrected with
the bonferonni procedure and I described comparisons that answered my
hypotheses.
What do you think about this analysis ?


Here is my syntax :
* Generalized Estimating Equations.
GENLIN choice (REFERENCE=LAST) BY Group Formulation Risk (ORDER=ASCENDING)

/MODEL GROUP FORMULATION RISK GROUP*FORMULATION GROUP*RISK FORMULATION*RISK      
GROUP*FORMULATION*RISK INTERCEPT=YES
DISTRIBUTION=BINOMIAL LINK=LOGIT
  /CRITERIA METHOD=FISHER(1) SCALE=1 MAXITERATIONS=100 MAXSTEPHALVING=5
PCONVERGE=1E-006(ABSOLUTE)
SINGULAR=1E-012 ANALYSISTYPE=3(WALD) CILEVEL=95 LIKELIHOOD=FULL
  /EMMEANS TABLES=GROUP*FORMULATION*RISK SCALE=ORIGINAL
COMPARE=GROUP*FORMULATION*RISK CONTRAST=PAIRWISE PADJUST=BONFERRONI
  /REPEATED SUBJECT=ID WITHINSUBJECT=Trial SORT=YES CORRTYPE=EXCHANGEABLE
ADJUSTCORR=YES COVB=MODEL MAXITERATIONS=100 PCONVERGE=1e-006(ABSOLUTE)
UPDATECORR=1
  /MISSING CLASSMISSING=EXCLUDE


Thank you again for your amazing support!




--
View this message in context: http://spssx-discussion.1045642.n5.nabble.com/Interactions-binary-logistic-regression-tp5733602p5733733.html
Sent from the SPSSX Discussion mailing list archive at 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

===================== 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