Description Usage Arguments Value References See Also Examples
This function fits the multinomial mixed model with two independent random effects
for each category of the response variable: one domain random effect and another correlated time and domain random effect (Model 3). The formulation is described in Lopez-Vizcaino et al. (2013).
The fitting algorithm combine the penalized quasi-likelihood method (PQL) for estimating
and predicting the fixed and random effects, respectively, with the residual maximun likelihood method (REML)
for estimating the variance components. This function uses as initial values the output of the function
initial.values
.
1 |
d |
number of areas. |
t |
number of time periods. |
pp |
vector with the number of the auxiliary variables per category. |
Xk |
list of matrices with the auxiliary variables per category obtained from |
X |
list of matrices with the auxiliary variables obtained from |
Z |
design matrix of random effects obtained from |
initial |
output of the function |
y |
matrix with the response variable obtained from |
M |
vector with the area sample sizes. |
MM |
vector with the population sample sizes. |
b |
parameter that indicates the bootstrap. |
A list containing the following components.
Estimated.probabilities |
matrix with the estimated probabilities for the categories of response variable. |
Fisher.information.matrix.phi |
Fisher information matrix of phi. |
Fisher.information.matrix.beta |
Fisher information matrix of beta. |
u1 |
matrix with the estimated first random effect. |
u2 |
matrix with the estimated second random effect. |
mean |
matrix with the estimated mean of the response variable. |
warning1 |
0=OK,1=The model could not be fitted. |
warning2 |
0=OK,1=The value of the variance component is negative: the initial value is taken. |
beta.Stddev.p.value |
matrix with the estimated fixed effects, its standard deviations and its p-values. |
phi.Stddev.p.value |
matrix with the estimated variance components, its standard deviations and its p-values. |
rho |
estimated correlation parameter. |
rho.Stddev.p.value |
matrix with the estimated correlation parameter, its standard deviations and its p-values. |
Lopez-Vizcaino, ME, Lombardia, MJ and Morales, D (2013). Small area estimation of labour force indicators under a multinomial mixed model with correlated time and area effects. Submitted for review.
data.mme
, initial.values
,
wmatrix
, phi.mult.ct
,
prmu.time
, phi.direct.ct
,
sPhikf.ct
, omega
, ci
,
Fbetaf.ct
, msef.ct
,
mseb
1 2 3 4 5 6 7 8 9 10 11 12 | ## Not run:
k=3 #number of categories of the response variable
pp=c(1,1) #vector with the number of auxiliary variables in each category
mod=3 #type of model
data(simdata3) #data
datar=data.mme(simdata3,k,pp,mod)
##Model fit
result=modelfit3(datar$d,datar$t,pp,datar$Xk,datar$X,datar$Z,datar$initial,datar$y[,1:(k-1)],
datar$n,datar$N,0)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.