inibestcubecov: Preliminary parameter estimates for CUBE models with...

Description Usage Arguments Value See Also Examples

Description

Compute preliminary parameter estimates for a CUBE model with covariates for all the three parameters. These estimates are set as initial values to start the E-M algorithm within maximum likelihood estimation.

Usage

1
inibestcubecov(m,ordinal,Y,W,Z)

Arguments

m

Number of ordinal categories

ordinal

Vector of ordinal responses (factor type)

Y

Matrix of selected covariates to explain the uncertainty parameter

W

Matrix of selected covariates to explain the feeling parameter

Z

Matrix of selected covariates to explain the overdispersion parameter

Value

A vector (inibet, inigama, inialpha) of preliminary estimates of parameter vectors for π = π(\bold{β}), ξ=ξ(\bold{γ}), φ=φ(\bold{α}), respectively, of a CUBE model with covariates for all the three parameters. In details, inibet, inigama and inialpha have length equal to NCOL(Y)+1, NCOL(W)+1 and NCOL(Z)+1, respectively, to account for an intercept term for each component

See Also

inibestcube, inibestcubecsi, inibestgama

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
data(relgoods)
attach(relgoods)
m<-10
naord<-which(is.na(Tv))
nacovpai<-which(is.na(Gender))
nacovcsi<-which(is.na(year.12))
nacovphi<-which(is.na(EducationDegree))
na<-union(union(naord,nacovpai),union(nacovcsi,nacovphi))
ordinal<-Tv[-na]
Y<-Gender[-na]
W<-year.12[-na]
Z<-EducationDegree[-na]
ini<-inibestcubecov(m,ordinal,Y,W,Z)
p<-NCOL(Y)
q<-NCOL(W)
inibet<-ini[1:(p+1)]               # Preliminary estimates for uncertainty 
inigama<-ini[(p+2):(p+q+2)]        # Preliminary estimates for feeling 
inialpha<-ini[(p+q+3):length(ini)] # Preliminary estimates for overdispersion


Search within the CUB package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.