Hello, I am currently trying to learn how to run a logistic regression with three main effects and two interaction terms in R studio for a study while comparing my results to my mentor’s results in SPSS. However, we seem to be getting different results for some the main effects, but the same results for the interaction terms. I think that I am writing the glm() code wrong, but I am brand new to R and not sure how to fix it. Please know I have spent HOURS trying to figure this out… and could really use your help!
For this study, we are researching how gender (moderator: Gend) moderates the effect of people’s endorsement of sexual scripts (first predictor: sexualscriptendorsementsmean) and attitudes towards casual sex (second predictor: casualsexscore) on their desire to engage in consensual non-monogamy (outcome: desire.yes_no).
We also want to run a second model with the same moderator and predictors, but with actual engagement as the outcome (second outcome: engagement.yes_no).
I attached my R studio syntax and output.
We have the exact same Beta, std error, and wald z-stat for the main effect of gender and the two interaction terms, but our main effects of casual sex attitudes and sexual script endorsement are different. The same thing happens for our second model with the outcome of engagement, to a greater extent as my output shows results for the main effects of casual sex and sexual scripts that are not significant but they are significant for the SPSS output.