Skip to contents

Create an instance of Sigmas class.

Sigmas.Rd

This function creates a new instance of different error models object to be applied. 0s are treated as no values.

Usage

Sigmas(
  Additive = 0,
  LogAdditive = 0,
  Proportional = 0.1,
  AdditiveMultiplicative = list(PropPart = 0, AddPart = 0),
  MixRatio = list(PropPart = 0, AddPart = 0),
  Power = list(PowerPart = 0, StdevPart = 0),
  ObservationName = ""
)

Arguments

Additive

The additive error sigma value.

LogAdditive

The log-additive error sigma value.

Proportional

The proportional error sigma value.

AdditiveMultiplicative

A list specifying the additive and multiplicative parts for the additive-multiplicative error model. The list should have elements PropPart and AddPart. Alternatively the proportional part (PropPart) could be presented as StParm, see StParm().

MixRatio

A list specifying the proportional and additive parts for the mix-ratio error model. The list should have elements PropPart and AddPart. Alternatively the proportional part (PropPart) could be presented as StParm, see StParm().

Power

A numeric vector specifying the standard deviation and power parts for the power error model. The vector should have names StdevPart and PowerPart.

ObservationName

A character string giving the name of the Observation.

See also

Functions used for Observation specification: Observation(), ObservationCustom(), create_ModelPD(), create_ModelPK(), modify_Observation(), remove_Observation()

Examples

RSE_CObs <-
  Observation(SigmasChosen =
    Sigmas(MixRatio = list(PropPart = 2,
                           AddPart = 0.01),
           Proportional = 0))
models <-
  create_ModelPK(CompartmentsNumber = 2,
                 CObs = RSE_CObs)
print(models)
#> PK2IVC 
#>  test() {
#> 	cfMicro(A1, Cl / V, Cl2 / V, Cl2 / V2)
#> C = A1 / V
#> 	dosepoint(A1, idosevar = A1Dose, infdosevar = A1InfDose, infratevar = A1InfRate)
#> 	error(CEps = 0.01)
#> 	observe(CObs = C + CEps*(1 + C*CMixRatio))
#> 	stparm(CMixRatio = tvCMixRatio )
#> 	fixef(tvCMixRatio= c(, 2, ))
#> 	
#> 	
#> 	stparm(Cl = tvCl * exp( nCl ))
#> 	fixef(tvCl= c(, 1, ))
#> 	ranef(diag(nCl) = c(1))
#> 	stparm(V = tvV * exp( nV ))
#> 	fixef(tvV= c(, 1, ))
#> 	ranef(diag(nV) = c(1))
#> 	stparm(Cl2 = tvCl2 * exp( nCl2 ))
#> 	fixef(tvCl2= c(, 1, ))
#> 	ranef(diag(nCl2) = c(1))
#> 	stparm(V2 = tvV2 * exp( nV2 ))
#> 	fixef(tvV2= c(, 1, ))
#> 	ranef(diag(nV2) = c(1))
#> 
#> }