compute.secondary function

Computes secondary parameter values

compute.secondary is a secondary function called during estimations runs. It evaluates the code provided in the $SECONDARY block of the model file; all parameters defined in this block are considered as secondary parameters at the initial and the final estimates of the model parameters. compute.secondary is typically not called directly by users.

compute.secondary(subproblem = NULL,
                    x = NULL)

Arguments

• subproblem: A list containing the following levels:

  • code: A list of R code extracted from the model file. Depending on content of the model file, the levels of this list could be: template, derived, lags, ode, dde, output, variance, and/or secondary.

  • method: A character string, indicating the scale of the analysis. Should be 'population' or 'subject'.

  • init: A data.frame of parameter data with the following columns: 'names', 'type', 'value', 'isfix', 'lb', and 'ub'.

  • debugmode: Logical indicator of debugging mode.

  • modfun: Model function.

  • data: A list containing as many levels as there are treatment levels for the subject (or population) being evaluated, plus the trts

     level listing all treatments for this subject (or population), and the `id` level giving the identification number of the subject (or set to 1 if the analysis was run at the level of the population.
     
     Each treatment-specific level is a list containing the following levels:
     
     - **ana**: mij x 3 data.frame containing the times of observations of the dependent variables (extracted from the TIME variable), the indicators of the type of dependent variables (extracted from the CMT variable), and the actual dependent variable observations (extracted from the DV variable) for this particular treatment.
     - **cov**: mij x c data.frame containing the times of observations of the dependent variables (extracted from the TIME variable) and all the covariates identified for this particular treatment.
     - **bolus**: bij x 4 data.frame providing the instantaneous inputs for a treatment and individual.
     - **infusion**: fij x (4+c) data.frame providing the zero-order inputs for a treatment and individual.
     - **trt**: the particular treatment identifier. 
    

• x: The vector of p final parameter estimates.

Returns

Return a list of with the following elements:

• init: The vector of s secondary parameter estimates derived from initial structural model parameter estimates.
• estimates: The vector of s secondary parameter estimates derived from final structural model parameter estimates.
• names: The vector of s secondary parameter names.

Author(s)

Sebastien Bihorel (sb.pmlab@gmail.com )

See Also

scarabee.analysis, weighting, fitmle, get.secondary