Description Usage Arguments Details Value Author(s) See Also Examples
Calculates the value of the Bivariate Response to Additive Interacting Doses (BRAID) surface model for the given concentration pairs.
1 | evalBRAIDrsm(DA, DB, parv)
|
DA |
vector of concentrations of drug A |
DB |
vector of concentrations of drug B |
parv |
ten-element vector specifying the full set of parameters for the BRAID surface (see Details below) |
The full ten-parameter BRAID model, which we refer to as the extended BRAID or eBRAID model is defined as:
E(DA,DB) = E0 + (Ef-E0)/(1+<DAB>^(-δ*√(na*nb)))
<DAB> = <DA>^(1/(δ*√(na*nb))) +<DB>^(1/(δ*√(na*nb)))+κ*√((<DA>^(1/(δ*√(na*nb))))*(<DB>^(1/(δ*√(na*nb)))))
<DA> = (((EfA-E0)/(Ef-E0))*((DA/IDMA)^na))/(1+(1-(EfA-E0)/(Ef-E0))*((DA/IDMA)^na))
<DB> (((EfB-E0)/(Ef-E0))*((DB/IDMB)^nb))/(1+(1-(EfB-E0)/(Ef-E0))*((DB/IDMB)^nb))
The parameters of this equation must satisfy na>0, nb>0, δ>0, κ> -2, sign(Ef-E0)=sign(EfA-E0)=sign(EfB-E0), |Ef-E0|≥|EfA-E0|, and |Ef-E0|≥|EfB-E0|. With this definition, the ten-element parameter vector is [IDMA, IDMB, na, nb, δ, κ, E0, EfA, EfB, Ef]. The simpler standard BRAID model, as described in Twarog et al. is obtained by setting δ equal to 1 and setting Ef such that |Ef-E0| is equal to the maximum of |EfA-E0| and |EfB-E0|. Assuming that this sets Ef equal to EfA, this causes the equation for <DA> to simplify to
<DA> = (DA/IDMA)^na
A vector of response values corresponding to the pairs of concentrations in DA
and DB
Nathaniel R. Twarog
1 2 3 4 5 6 7 8 9 10 11 12 13 | conc1 <- rep(seq(0,3*10^-6,length=50),each=50)
conc2 <- rep(seq(0,3*10^-6,length=50),times=50)
# Additive surface
act <- evalBRAIDrsm(conc1,conc2,parv=c(10^-6,10^-6,1.5,1.5,1,0,0,100,100,100))
# A BRAID additive surface is not a Loewe additive surface
act <- evalBRAIDrsm(conc1,conc2,parv=c(10^-6,10^-6,1,3,1,0,0,100,100,100))
# BRAID antagonism
act <- evalBRAIDrsm(conc1,conc2,parv=c(10^-6,10^-6,1.5,1.5,1,-1,0,100,100,100))
# delta-BRAID synergy
act <- evalBRAIDrsm(conc1,conc2,parv=c(10^-6,10^-6,1.5,1.5,1.75,0,0,100,100,100))
# Differing final effects
act <- evalBRAIDrsm(conc1,conc2,parv=c(10^-6,10^-6,1.5,1.5,1,0,0,75,100,100))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.