indParm: 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)

Value

a list with phi_i and psi_i parameters.

Arguments

theta

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

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.

Examples

Run this code

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)