lrt.hopit | R Documentation |
Likelihood ratio test for a pair of models
lrt.hopit(full, nested)
full, nested |
models to be compared. |
a vector with the results of the test.
Maciej J. Danko
print.lrt.hopit
, anova.hopit
,
hopit
.
# DATA data(healthsurvey) # the order of response levels decreases from the best health to # the worst health; hence the hopit() parameter decreasing.levels # is set to TRUE levels(healthsurvey$health) # Example 1 --------------------- # fitting two nested models model1 <- hopit(latent.formula = health ~ hypertension + high_cholesterol + heart_attack_or_stroke + poor_mobility + very_poor_grip + depression + respiratory_problems + IADL_problems + obese + diabetes + other_diseases, thresh.formula = ~ sex + ageclass + country, decreasing.levels = TRUE, control = list(trace = FALSE), data = healthsurvey) # model with an interaction between hypertension and high_cholesterol model2 <- hopit(latent.formula = health ~ hypertension * high_cholesterol + heart_attack_or_stroke + poor_mobility + very_poor_grip + depression + respiratory_problems + IADL_problems + obese + diabetes + other_diseases, thresh.formula = ~ sex + ageclass + country, decreasing.levels = TRUE, control = list(trace = FALSE), data = healthsurvey) # Likelihood ratio test lrt1 <- lrt.hopit(full = model2, nested = model1) lrt1 # print the results in a shorter form print(lrt1, short = TRUE) # equivalently print(anova(model2, model1), short = TRUE)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.