nrmlepln {pln} | R Documentation |
Full information maximum likelihood and bivariate composite likelihood estimation for polytomous logit-normit (graded logistic) model
Description
Full information maximum likelihood and bivariate composite likelihood estimation for polytomous logit-normit and Rasch models, via Newton Raphson iterations.
Usage
nrmlepln(x, ncat, nitem=NULL, alphas=NULL, betas=NULL, abound=c(-10,10),
bbound=c(-1,10), nq=48, mxiter=200, m2=TRUE, iprint=FALSE)
nrmlerasch(x, ncat, nitem=NULL, alphas=NULL, abound=c(-10,10),
bbound=c(-1,10), nq=48, mxiter=200, m2=TRUE, iprint=FALSE)
nrbcpln(x, ncat, nitem=NULL, alphas=NULL, betas=NULL, abound=c(-10,10),
bbound=c(-1,10), nq=48, mxiter=200, se=TRUE, iprint=FALSE)
Arguments
x |
A data matrix. Data can be in one of two formats: 1) raw data
where the number of rows corresponds to an individual's response and
each column represents an item, and 2) a matrix of dimensions
|
ncat |
Number of ordinal categories for each item, coded as
0,...,( |
nitem |
Number of items. If omitted, it is assumed that |
alphas |
A vector of length |
betas |
A vector of length |
abound |
Vector of length 2 that sets upper and lower bounds on parameter estimation for alphas. Currently experimental; changing defaults it not recommended. Estimation problems are more likely solved by changing starting values. |
bbound |
Vector of length 2 that sets upper and lower bounds on parameter estimation for betas. Currently experimental; changing defaults it not recommended. Estimation problems are more likely solved by changing starting values. |
nq |
Number of quadrature points to use during estimation. This argument is currently experimental. It is recommended to use the default of 48. |
mxiter |
Maximum number of iterations for estimation. |
m2 |
Logical. If |
iprint |
Logical. Enables debugging / diagnostic information from C code that conducts estimation. |
se |
Logical. If |
Details
Estimation of graded logistic models is performed under the following parameterization:
Where the items are , and response categories are
.
is the latent trait,
is the logistic distribution function,
is an intercept (cutpoint) parameter, and
is a slope parameter. When the number of categories for the items is 2, this reduces to the 2PL parameterization:
Both nrmlepln
and nrbcpln
perform estimation under these parameterizations, via Newton Raphson iterations, using full information maximum likelihood (nrmlepln
) and bivariate composite likelihood (nrbcpln
). See Maydeu-Olivares and Joe (2005, 2006) for more information on bivariate composite likelihood estimation (see also Varin, Reid, and Firth, 2011). Under nrmlerasch
a common parameter is estimated for all items.
Value
A list containing the following slots.
Slots
alphas
A vector of parameter estimates for alphas. Length is
nitem
(
ncat
-1). Estimates are in order by item, e.g., all alphas for item 1, followed by all alphas for item 2, and so on.betas
A vector of parameter estimates for betas. Length is
nitem
.nllk
Negative (composite) log-likelihood for polytomous logit-normit (or Rasch) model.
conv
Integer indicating whether estimation converged. Currently only returned for composite likelihood estimation.
sealphas
A vector of standard errors for the alpha estimates.
sebetas
A vector of standard errors for the beta estimates.
invhes
Inverse Hessian matrix for the MLE estimates.
vcov
Asymptotic covariance matrix for the composite likelihood estimates.
teststat
Value of
.
df
Degrees of freedom for
.
pval
P-value for
.
Author(s)
Carl F. Falk cffalk@gmail.com, Harry Joe
References
Bartholomew, D., Knott, M., and Moustaki, I. (2011). Latent Variable Models and Factor Analysis: A Unified Approach, 3rd Edition. Wiley.
Maydeu-Olivares, A., and Joe, H. (2005). Limited and full information estimation
and goodness-of-fit testing in contingency tables: A unified framework.
Journal of the American Statistical Association, 100, 1009-1020.
Maydeu-Olivares, A., and Joe, H. (2006). Limited information and goodness-of-fit testing in multidimensional contingency tables. Psychometrika, 71, 713-732.
Varin, C., Reid, N. and Firth, D. (2011). An overview of composite likelihood methods. Statistica Sinica, 21, 5-42.
See Also
Examples
### Matrix of response patterns and frequencies
data(item5fr)
## ML estimation
nrmleplnout<-nrmlepln(item5fr, ncat=3, nitem=5)
print(nrmleplnout)
## BCL estimation
nrbcplnout<-nrbcpln(item5fr, ncat=3, nitem=5)
print(nrbcplnout)
## ML Rasch estimation
nrmleraschout<-nrmlerasch(item5fr, ncat=3, nitem=5)
print(nrmleraschout)
### Raw data
data(item9cat5)
## ML estimation
nrmleplnout<-nrmlepln(item9cat5, ncat=5)
print(nrmleplnout)
## BCL estimation
nrbcplnout<-nrbcpln(item9cat5, ncat=5, se=FALSE)
print(nrbcplnout)
## ML Rasch estimation
nrmleraschout<-nrmlerasch(item9cat5, ncat=5)
print(nrmleraschout)