profile.hopit | R Documentation |
hopit
modelCalculate the log likelihood profile for the fitted hopit
model
## S3 method for class 'hopit'
profile(fitted, ..., scope = 0.15, steps = 101)
fitted |
a |
... |
unused now. |
scope |
a value (fraction) defining the plotting range for a coefficient.
The range is |
steps |
at how many equally spaced points the log likelihood function is calculated for each coefficient. |
Maciej J. Danko
plot.profile.hopit
, print.profile.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 the model:
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)
# check the fit using the profile function (at 51 points)
pr <- profile(model1, steps = 51)
print(pr, plotf = FALSE)
# plot profile
plot(pr, relative = FALSE)
# alternative plot
plot(pr, relative = TRUE)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.