[R-sig-ME] logistic regression on posttest (0, 1) with pretest(0, 1)*Group(Treatment, Ctrl) interaction

Sorkin, John j@ork|n @end|ng |rom @om@um@ry|@nd@edu
Sun Apr 21 20:01:59 CEST 2019


It is both difficult and dangerous to add a comment to a thread that one has not followed, and in doing so possibly making an inappropriate suggestion. Please forgive what may be an not fully informed thought.

The model you suggest, posttest ~ pretest*Group  (ignoring random effects) is unusual. In a model that contains an interaction,  I would expect to see, in addition to the interaction, all main effects included in the interaction, i.e.
posttest ~ group+pretest+pretest*Group


John David Sorkin M.D., Ph.D.
Professor of Medicine
Chief, Biostatistics and Informatics
University of Maryland School of Medicine Division of Gerontology and Geriatric Medicine
Baltimore VA Medical Center
10 North Greene Street
Baltimore, MD 21201-1524
(Phone) 410-605-7119
(Fax) 410-605-7913 (Please call phone number above prior to faxing)

From: R-sig-mixed-models <r-sig-mixed-models-bounces using r-project.org> on behalf of Souheyla GHEBGHOUB <souheyla.ghebghoub using gmail.com>
Sent: Sunday, April 21, 2019 4:57 AM
To: Ren�; r-sig-mixed-models using r-project.org
Subject: [R-sig-ME] logistic regression on posttest (0, 1) with pretest(0, 1)*Group(Treatment, Ctrl) interaction

Dear Rene, and any member of the list who is willing to read this : )

I have decided to use the interesting model that you have structured for me,
To refresh your memory here is what you said :

*Anyway, there is a second possibility to define your model, depending on
> how you want to interpret it. In the previous model you can say something
> about the type-of-change likelihoods depending on the treatment group.
> But you could implement the model as binomial as well (i.e. logistic
> regression) mod2 <- brm(posttest ~ pretest*Group + (1|Subject) +
> (1+Group|Word),...)  And what you would expect here would be an interaction
> between pre-test and Group. For instance; if pretest=0 & treatment 1 then
> posttest larger than with pretest=0 & treatment 2; but not when pretest=1
> (because this is a plausible no gain situation). And so on...*

But I found the interpretation of the Pretest*Group interaction results to
be tough. I still can't grasp how to claim there is an effect of Group on
posttest outcome when considering pretest. The pretest slope does not come
with 0 or 1 in the output, the Group does come with one category, but its
confusing what the intercept and slope estimates refer to in this case?

Thank you very much,

        [[alternative HTML version deleted]]

R-sig-mixed-models using r-project.org mailing list

	[[alternative HTML version deleted]]

More information about the R-sig-mixed-models mailing list