User-defined dose-response function
duser(fun, beta.1 = "rel", beta.2 = "rel", beta.3 = "rel", beta.4 = "rel")
An object of class("dosefun")
A formula specifying any relationship including dose
and
one/several of: beta.1
, beta.2
, beta.3
, beta.4
.
Pooling for the 1st coefficient. Can take "rel"
, "common"
, "random"
or be
assigned a numeric value (see details).
Pooling for the 2nd coefficient. Can take "rel"
, "common"
, "random"
or be
assigned a numeric value (see details).
Pooling for the 3rd coefficient. Can take "rel"
, "common"
, "random"
or be
assigned a numeric value (see details).
Pooling for the 4th coefficient. Can take "rel"
, "common"
, "random"
or be
assigned a numeric value (see details).
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()
.
dr <- ~ beta.1 * (1/(dose+1)) + beta.2 * dose^2
duser(fun=dr,
beta.1="common", beta.2="rel")
Run the code above in your browser using DataLab