Remove Observation from PML models

Usage

remove_Observation(PMLParametersSets, ObservationName, PMLStructures = NULL)

Arguments

PMLParametersSets: A list of PML parameters sets (PMLModels class instance).
ObservationName: A character string giving the name of the Observation.
PMLStructures: Character or character vector specifying names of PML structures from which the observation will be removed. For the naming convention of PMLStructures, see Details section of create_ModelPK() for PK models and create_ModelPD() for PD models.

Value

An updated list of PML models (PMLModels class instance) matching the specified options.

Details

The current functionality does not support modifying custom observations that are defined within the PML code of custom model spaces.

Examples

PMLParametersSets <-
  create_ModelPK(
    CompartmentsNumber = c(2, 3),
    Parameterization = "Micro",
    Absorption = c("First-Order", "Gamma"),
    ByVector = TRUE,
    ClosedForm = TRUE,
    EliminationCpt = TRUE)

remove_Observation(PMLParametersSets,
                   ObservationName = "A0Obs",
                   PMLStructures = "PK3GME")
#> PK2FOME 
#>  test() {
#> 	deriv(Aa = - Ka * Aa)
#> 	deriv(A1 = Ka * Aa - Ke * A1 - K12 * A1 + K21 * A2)
#> 	deriv(A2 = K12 * A1 - K21 * A2)
#> 	urinecpt(A0 = Ke * A1)
#> 	C = A1 / V
#> 	dosepoint(Aa, idosevar = AaDose, infdosevar = AaInfDose, infratevar = AaInfRate)
#> 	error(A0Eps = 0.1)
#> 	observe(A0Obs = A0 * (1 + A0Eps))
#> 	error(CEps = 0.1)
#> 	observe(CObs = C * (1 + CEps))
#> 	
#> 	stparm(Ka = tvKa * exp( nKa ))
#> 	fixef(tvKa= c(, 1, ))
#> 	ranef(diag(nKa) = c(1))
#> 	stparm(Ke = tvKe * exp( nKe ))
#> 	fixef(tvKe= c(, 1, ))
#> 	ranef(diag(nKe) = c(1))
#> 	stparm(K12 = tvK12 * exp( nK12 ))
#> 	fixef(tvK12= c(, 1, ))
#> 	ranef(diag(nK12) = c(1))
#> 	stparm(K21 = tvK21 * exp( nK21 ))
#> 	fixef(tvK21= c(, 1, ))
#> 	ranef(diag(nK21) = c(1))
#> 	stparm(V = tvV * exp( nV ))
#> 	fixef(tvV= c(, 1, ))
#> 	ranef(diag(nV) = c(1))
#> 
#> } 
#> PK3GME 
#>  test() {
#> 	delayInfCpt(A1, MeanDelayTime, ShapeParamMinusOne, out = - Ke * A1 - K12 * A1 + K21 * A2 - K13 * A1 + K31 * A3, dist = Gamma)
#> 	deriv(A2 = K12 * A1 - K21 * A2)
#> 	deriv(A3 = K13 * A1 - K31 * A3)
#> 	urinecpt(A0 = Ke * A1)
#> 	C = A1 / V
#> 	dosepoint(A1, idosevar = A1Dose, infdosevar = A1InfDose, infratevar = A1InfRate)
#> 	error(CEps = 0.1)
#> 	observe(CObs = C * (1 + CEps))
#> 	
#> 	stparm(MeanDelayTime = tvMeanDelayTime * exp( nMeanDelayTime ))
#> 	fixef(tvMeanDelayTime= c(, 1, ))
#> 	ranef(diag(nMeanDelayTime) = c(1))
#> 	stparm(ShapeParamMinusOne = tvShapeParamMinusOne * exp( nShapeParamMinusOne ))
#> 	fixef(tvShapeParamMinusOne= c(, 1, ))
#> 	ranef(diag(nShapeParamMinusOne) = c(1))
#> 	stparm(Ke = tvKe * exp( nKe ))
#> 	fixef(tvKe= c(, 1, ))
#> 	ranef(diag(nKe) = c(1))
#> 	stparm(K12 = tvK12 * exp( nK12 ))
#> 	fixef(tvK12= c(, 1, ))
#> 	ranef(diag(nK12) = c(1))
#> 	stparm(K21 = tvK21 * exp( nK21 ))
#> 	fixef(tvK21= c(, 1, ))
#> 	ranef(diag(nK21) = c(1))
#> 	stparm(K13 = tvK13 * exp( nK13 ))
#> 	fixef(tvK13= c(, 1, ))
#> 	ranef(diag(nK13) = c(1))
#> 	stparm(K31 = tvK31 * exp( nK31 ))
#> 	fixef(tvK31= c(, 1, ))
#> 	ranef(diag(nK31) = c(1))
#> 	stparm(V = tvV * exp( nV ))
#> 	fixef(tvV= c(, 1, ))
#> 	ranef(diag(nV) = c(1))
#> 
#> }