Description Usage Arguments Value References See Also Examples

This function calculates the score vector S and the inverse of the Fisher information
matrix for the fixed (beta) and the random effects (u1, u2) in Model 3. This model has two independet sets of random effects.
The first one contains independent random effects u1dk associated to each category and domain. The second set contains random effects
u2dkt associated to each category, domain and time period. Model 3 assumes that the u2dk are AR(1) correlated across time.
`modelfit3`

uses the output of this function to estimate the fixed and random effect by the PQL method.

1 | ```
Fbetaf.ct(sigmap, X, Z, phi1, phi2, y, mu, u1, u2, rho)
``` |

`sigmap` |
a list with the model variance-covariance matrices for each domain obtained from |

`X` |
list of matrices with the auxiliary variables obtained from |

`Z` |
design matrix of random effects. |

`phi1` |
vector with the values of the variance components for the first random effects obtained from |

`phi2` |
vector with the values of the variance components for the second random effects obtained from |

`y` |
matrix with the response variable, except the reference category. The rows are the domains and the columns are the categories of the response variable minus one. |

`mu` |
matrix with the estimated mean of the response variable. |

`u1` |
matrix with the values of the first random effect obtained from |

`u2` |
matrix with the values of the second random effect obtained from |

`rho` |
vector with the values of the correlation parameter obtained from |

A list containing the following components.

`F` |
the inverse of the Fisher information matrix of (beta, u1, u2). |

`S` |
(beta, u1, u2) score vectors |

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`

, `ci`

,
`modelfit3`

, `msef.ct`

, `omega`

,
`mseb`

1 2 3 4 5 6 7 8 9 10 11 12 | ```
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)
datar=data.mme(simdata3,k,pp,mod)
initial=datar$initial
mean=prmu.time(datar$n,datar$Xk,initial$beta.0,initial$u1.0,initial$u2.0)
sigmap=wmatrix(datar$n,mean$estimated.probabilities) #variance-covariance
##The inverse of the Fisher information matrix and the score matrix
Fisher.beta=Fbetaf.ct(sigmap,datar$X,datar$Z,initial$phi1.0,initial$phi2.0,
datar$y[,1:(k-1)],mean$mean,initial$u1.0,initial$u2.0,initial$rho.0)
``` |

mme documentation built on May 30, 2017, 3:38 a.m.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.