indParm: Generate individual parameters

In REMixed: Regularized Estimation in Mixed Effects Model

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)

REMixed documentation built on Jan. 19, 2026, 5:07 p.m.