After applying matching on multiple imputed datasets. How do I set the reference levels for categorical variables in my final model. For example how do I set the reference level for KOA as 1 instead of 0?

library(MatchThem)

data("osteoarthritis")

library(mice)
imputed.datasets <- mice(osteoarthritis, m = 5)

table(imputed.datasets$data$OSP)

table(imputed.datasets$data$KOA)

matched.datasets <- matchthem(OSP ~ AGE + SEX + BMI + RAC + SMK,
                              datasets = imputed.datasets,
                              approach = 'within',
                              method = 'nearest',
                              caliper = 0.05,
                              ratio = 2)

library(survey)
matched.models <- with(matched.datasets,
                       svyglm(KOA ~ OSP, family = quasibinomial()),
                       cluster = TRUE)

After applying matching on multiple imputed datasets. How do I set the reference levels for categorical variables in my final model. For example how do I set the reference level for KOA as 1 instead of 0?

library(MatchThem)

data("osteoarthritis")

library(mice)
imputed.datasets <- mice(osteoarthritis, m = 5)

table(imputed.datasets$data$OSP)

table(imputed.datasets$data$KOA)

matched.datasets <- matchthem(OSP ~ AGE + SEX + BMI + RAC + SMK,
                              datasets = imputed.datasets,
                              approach = 'within',
                              method = 'nearest',
                              caliper = 0.05,
                              ratio = 2)

library(survey)
matched.models <- with(matched.datasets,
                       svyglm(KOA ~ OSP, family = quasibinomial()),
                       cluster = TRUE)

• r
• matching
• propensity-score-matching
• causal-inference

Share Improve this question Follow edited Mar 31 at 20:12 desertnaut 60.5k3232 gold badges155155 silver badges182182 bronze badges asked Nov 15, 2024 at 21:51 Science11Science11 88311 gold badge1111 silver badges2626 bronze badges

Add a comment  | 

1 Answer 1

Sorted by: Reset to default 2

You can just make the outcome KOA == 0, e.g.,

svyglm(KOA == 0 ~ OSP, family = quasibinomial())

Note this doesn't really have anything to do with matching or multiple imputation.