1
$\begingroup$

I am trying to use GLMMs models to analysis the morbidity status of child (yes or no) with mother’s demographic and environmental factors like

  1. Wealth with factors ("Lower quartiles”,"Middle quartiles") and “Higher quartiles"
  2. Residence (Urban" and Rural")
  3. Education ("Never attended”, “Primary education", “Secondary education", "Above secondary education")
  4. Marital status (“Married or with partner", Widowed or divorced”, ” Never married"))
  5. Age (Age between 15-24”, “Age between 25-34", "Age above 34")
  6. Parity (“Zero Parity", “Parity between 1 and 2", “Parity between 3 and 4”, Parity above 5")

The following is r code formula that I try to fix determine the model

Fit1<-glmer (Morbidity~Wealth+Residencey+education+Marital+Age+Parity +(1|region)+(1|participant_id),family=binomial,data=M)

Dear respected professionals, I try to use GLMMs model, Then I face a bit more serious problems on the following issues like

Q1. How many interaction is possible between the given independent variables and how can specify it? If I try to test the interaction between certain variables, it is honestly insignificance.

Q2. The appropriate R code to plot the interaction effect between variables two variable and the responsice variable, morbidity. If somebody can help me on it, I can’t thank you enough your quick response.

$\endgroup$
3
  • $\begingroup$ I edited your post for clarity and English usage. I also formatted the code. $\endgroup$
    – Peter Flom
    Commented Oct 20, 2023 at 19:36
  • $\begingroup$ I'm confused. Region is where the participants are from no? Wouldn't that mean that participants are nested within region? If so, that would likely be a better modeling of the random effects than what you currently have. $\endgroup$ Commented Oct 22, 2023 at 10:49
  • $\begingroup$ This problem is not set up correctly. Age and parity were inappropriately binned. Go back to the raw data. Also, morbidity is not a binary variable. There are degrees of morbidity. Using the variables incorrectly will damage the analysis more any any proper analysis can fix. And interactions need to be pre-specified using subject matter considerations. Binning variables will also create false interactions that try to make up for the loss of information from the original variables. $\endgroup$ Commented Oct 22, 2023 at 13:04

1 Answer 1

0
$\begingroup$

Welcome to CV. You ask

Q1. How many interaction is possible between the given independent variables and how can specify it? If I try test the interaction between certain variables, it is honestly insignificant.

I don't know what "honestly insignificant" means, but the number of two way interactions among $p$ variables is $p(p-1)/2$. If you want to look at higher-way interactions, there are a lot more.

Q2. The appropriate R code to plot the interaction effect between variables two variable and the response variable, morbidity.

There is no such thing as a interaction between the independent variables and the response variable. Interactions are among independent variables.

You didn't ask about this, but categorizing quantitative variables (wealth and age) is usually a bad idea, unless you have some strong, substantive reason for doing so.

I'm also not clear why you need a GLMM (that is a mixed model). There doesn't seem to be any clustering, either in time or space or otherwise, and the errors would seem to be independent (but maybe there is more going on that you have not told us).

$\endgroup$
1
  • $\begingroup$ Dear Peter Flow, Honestly, I can't thank you enough for your help. The word I used, "honestly insignificant," is simply to say insignificant. Based on my thoughts, I try to use the GLMM to see the presence of regional and cluster heterogeneity in child morbidity, given that there are child morbidity records based on different regions of the country (6 regions). The one thing that I want to say is that, say, I have five predictors. Then how many interactions are possible between the independent variables? If there are many interactions, is there any method to reduce them? $\endgroup$ Commented Oct 21, 2023 at 19:07

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.