User-defined dose-response function

duser(fun, beta.1 = "rel", beta.2 = "rel", beta.3 = "rel", beta.4 = "rel")

Arguments

fun: A formula specifying any relationship including dose and one/several of: beta.1, beta.2, beta.3, beta.4.
beta.1: Pooling for the 1st coefficient. Can take "rel", "common", "random" or be assigned a numeric value (see details).
beta.2: Pooling for the 2nd coefficient. Can take "rel", "common", "random" or be assigned a numeric value (see details).
beta.3: Pooling for the 3rd coefficient. Can take "rel", "common", "random" or be assigned a numeric value (see details).
beta.4: Pooling for the 4th coefficient. Can take "rel", "common", "random" or be assigned a numeric value (see details).

Value

An object of class("dosefun")

Dose-response parameters

Argument	Model specification
`"rel"`	Implies that relative effects should be pooled for this dose-response parameter separately for each agent in the network.
`"common"`	Implies that all agents share the same common effect for this dose-response parameter.
`"random"`	Implies that all agents share a similar (exchangeable) effect for this dose-response parameter. This approach allows for modelling of variability between agents.
`numeric()`	Assigned a numeric value, indicating that this dose-response parameter should not be estimated from the data but should be assigned the numeric value determined by the user. This can be useful for fixing specific dose-response parameters (e.g. Hill parameters in Emax functions) to a single value.

When relative effects are modelled on more than one dose-response parameter, correlation between them is automatically estimated using a vague inverse-Wishart prior. This prior can be made slightly more informative by specifying the scale matrix omega and by changing the degrees of freedom of the inverse-Wishart prior using the priors argument in mbnma.run().

References

There are no references for Rd macro \insertAllCites on this help page.

Examples


dr <- ~ beta.1 * (1/(dose+1)) + beta.2 * dose^2
duser(fun=dr,
  beta.1="common", beta.2="rel")
#> $name
#> [1] "user"
#> 
#> $fun
#> ~beta.1 * (1/(dose + 1)) + beta.2 * dose^2
#> <environment: 0x55c46b3da910>
#> 
#> $params
#> [1] "beta.1" "beta.2"
#> 
#> $nparam
#> [1] 2
#> 
#> $jags
#> [1] "s.beta.1[agent[i,k]] * (1/(dose[i,k] + 1)) + s.beta.2[agent[i,k]] * dose[i,k]^2"
#> 
#> $apool
#>   beta.1   beta.2 
#> "common"    "rel" 
#> 
#> $bname
#> [1] "beta.1" "beta.2"
#> 
#> attr(,"class")
#> [1] "dosefun"