profilelike.polr {ProfileLikelihood} | R Documentation |
Profile Likelihood for Proportional Odds Models
Description
This function provides values for a profile likelihood and a normalized profile likelihood for a parameter of interest in a proportional odds model.
Usage
profilelike.polr(formula, data, profile.theta, method = "logistic",
lo.theta = NULL, hi.theta = NULL, length = 300, round = 2,
subset = NULL, weights = NULL, offset = NULL, ...)
Arguments
formula |
see corresponding documentation in |
data |
a data frame. See corresponding documentation in |
profile.theta |
a parameter of interest, theta; must be a numeric variable. |
method |
see corresponding documentation in |
lo.theta |
lower bound for a parameter of interest to obtain values for a profile likelihood. |
hi.theta |
upper bound for a parameter of interest to obtain values for a profile likelihood. |
length |
length of numerical grid values for a parameter of interest to obtain values for a profile likelihood. |
round |
the number of decimal places for |
subset |
should not be provided. |
weights |
should not be provided. |
offset |
should not be provided. |
... |
further arguments passed to or from other methods. |
Details
This function provides values for a profile likelihood and a normalized profile likelihood for a parameter of interest in a proportional odds model. Users must define a parameter of interest in a proportional odds model. This function can be used for proportional odds models comparable with the polr
function. However, arguments weights, subset, and offset should not be provided. A normalized profile likelihood is obtained by a profile likelihood being divided by the maximum value of the profile likelihood so that a normalized profile likelihood ranges from 0 to 1.
Value
theta |
numerical grid values for a parameter of interest in a specified range (between lower and upper bounds). |
profile.lik |
numerical values for a profile likelihood corresponding to theta in a specified range (between lower and upper bounds). |
profile.lik.norm |
numerical values for a normalized profile likelihood ranging from 0 to 1. |
Warning
Arguments weights, subset, and offset in the polr
function are not comparable.
Missing values should be removed.
Author(s)
Leena Choi <naturechoi@gmail.com>
See Also
profilelike.plot
, profilelike.summary
, profilelike.lm
, profilelike.glm
, profilelike.gls
, profilelike.lme
, polr
Examples
data(datapolr)
datapolr$y <- as.factor(datapolr$y)
xx <- profilelike.polr(y ~ x1 + x2, data=datapolr, profile.theta="group",
method="logistic", lo.theta=-2, hi.theta=2.5, length=500)
profilelike.plot(theta=xx$theta, profile.lik.norm=xx$profile.lik.norm, round=2)