indParm {REMixed}R Documentation

Generate individual parameters

Description

Generate the individual parameters of indivual whose covariates are covariates and random effects eta_i.

Usage

indParm(theta, covariates, eta_i, transfo, transfo.inv)

Arguments

theta

list with at least phi_pop, psi_pop, gamma, beta (named ; corresponding to the model parameter \phi_{pop}, \psi_{pop}, \gamma, \beta ) :

  • phi_pop named vector of population parameters without r.e ;

  • psi_pop named vector of population parameters with r.e ;

  • gamma named list of vector of covariates effects for phi_pop parameters, if NULL no covariates effect on parameters. ;

  • beta named list of vector of covariates effects for each psi_pop, if NULL no covariates effect on parameters.

covariates

line data.frame of individual covariates ;

eta_i

named vector of random effect for each psi parameter ;

transfo

named list of transformation functions (h_l)_{l\leq m} and (s_k)_{k\leq K} for the individual parameter model (names must be consistent with phi_pop and psi_pop, missing entries are set by default to the identity function).

transfo.inv

amed list of inverse transformation functions for the individual parameter model (names must be consistent withphi_pop andpsi_pop).

Details

The models used for the parameters are :

h_l(\psi_{li}) = h_l(\psi_{lpop})+X_i\beta_l + \eta_{li}

with h_l the transformation, \beta_l the vector of covariates effect and with \eta_i the random effects associated \psi_l parameter ;

g_k(\phi_{ki}) = g_k(\phi_{kpop})+X_i \gamma_l

with g_k the transformation and \gamma_k the vector of covariates effect associated \phi_k parameter.

Value

a list with phi_i and psi_i parameters.

See Also

model.clairon, model.pasin.

Examples

phi_pop = c(delta_S = 0.231, delta_L = 0.000316)
psi_pop = c(delta_Ab = 0.025,phi_S = 3057, phi_L = 16.6)
gamma = NULL
covariates = data.frame(cAGE = runif(1,-15,15), G1 = rnorm(1), G2 = rnorm(1))
beta = list(delta_Ab=c(0,1.2,0),phi_S = c(0.93,0,0),phi_L=c(0,0,0.8))

theta=list(phi_pop = phi_pop,psi_pop = psi_pop,gamma = gamma, beta = beta)
eta_i = c(delta_Ab = rnorm(1,0,0.3),phi_S=rnorm(1,0,0.92),phi_L=rnorm(1,0,0.85))
transfo = list(delta_Ab=log,phi_S=log,phi_L=log)
transfo.inv = list(delta_Ab = exp,phi_S=exp,phi_L=exp)

indParm(theta,covariates,eta_i,transfo,transfo.inv)

[Package REMixed version 0.1.0 Index]