You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
HLA_model <- glm(disease ~ term + sex + age, data = HLA, family = binomial)
HLA_results <- runMiDAS(
object = HLA_model,
experiment = "hla_alleles",
inheritance_model = "dominant",
lower_frequency_cutoff = 0.02,
upper_frequency_cutoff = 0.98,
correction = "bonferroni",
exponentiate = TRUE
)
Alternatively, an omnibus test for multiple HLA alleles
HLA_AA_omnibus_results <- runMiDAS(
HLA_model,
experiment = "hla_aa",
inheritance_model = "dominant",
conditional = FALSE,
omnibus = TRUE,
lower_frequency_cutoff = 0.02,
upper_frequency_cutoff = 0.98,
correction = "bonferroni"
)
In both cases, the output includes only the p-values or odds ratios (ORs) for the HLA alleles. My question is: Is it possible to output the p-values and ORs for the covariates (e.g., sex, age) as well?
Thanks!
The text was updated successfully, but these errors were encountered:
In the following examples,
HLA_model <- glm(disease ~ term + sex + age, data = HLA, family = binomial)
HLA_results <- runMiDAS(
object = HLA_model,
experiment = "hla_alleles",
inheritance_model = "dominant",
lower_frequency_cutoff = 0.02,
upper_frequency_cutoff = 0.98,
correction = "bonferroni",
exponentiate = TRUE
)
Alternatively, an omnibus test for multiple HLA alleles
HLA_AA_omnibus_results <- runMiDAS(
HLA_model,
experiment = "hla_aa",
inheritance_model = "dominant",
conditional = FALSE,
omnibus = TRUE,
lower_frequency_cutoff = 0.02,
upper_frequency_cutoff = 0.98,
correction = "bonferroni"
)
In both cases, the output includes only the p-values or odds ratios (ORs) for the HLA alleles. My question is: Is it possible to output the p-values and ORs for the covariates (e.g., sex, age) as well?
Thanks!
The text was updated successfully, but these errors were encountered: