gain | R Documentation |
Compute the gain value with a given dose level, given a pseudo DLE model, a DLE sample, a pseudo Efficacy log-log model and a Efficacy sample
gain(dose, DLEmodel, DLEsamples, Effmodel, Effsamples, ...) ## S4 method for signature 'numeric,ModelTox,Samples,Effloglog,Samples' gain(dose, DLEmodel, DLEsamples, Effmodel, Effsamples, ...) ## S4 method for signature 'numeric,ModelTox,Samples,EffFlexi,Samples' gain(dose, DLEmodel, DLEsamples, Effmodel, Effsamples, ...) ## S4 method for signature 'numeric,ModelTox,missing,Effloglog,missing' gain(dose, DLEmodel, DLEsamples, Effmodel, Effsamples, ...)
dose |
the dose |
DLEmodel |
the |
DLEsamples |
the |
Effmodel |
the |
Effsamples |
the |
... |
unused |
gain(
dose = numeric,
DLEmodel = ModelTox,
DLEsamples = Samples,
Effmodel = EffFlexi,
Effsamples = Samples
)
: Compute the gain given a dose level, a pseduo DLE model, a DLE sample,
the pseudo EffFlexi model and an Efficacy sample
gain(
dose = numeric,
DLEmodel = ModelTox,
DLEsamples = missing,
Effmodel = Effloglog,
Effsamples = missing
)
: Compute the gain value given a dose level, a pseudo DLE model and a pseudo
efficacy model of Effloglog
class object without DLE and the efficacy sample
##Obtain the gain value for a given dose, a pseudo DLE model, a DLE sample, ## a pseudo efficacy model and an efficacy sample ##The DLE model must be from 'ModelTox' class (DLEmodel slot) emptydata<- DataDual(doseGrid=seq(25,300,25),placebo=FALSE) data<-emptydata DLEmodel<-LogisticIndepBeta(binDLE=c(1.05,1.8),DLEweights=c(3,3),DLEdose=c(25,300),data=data) DLEsamples <- mcmc(data, DLEmodel, McmcOptions(burnin=100,step=2,samples=200)) ##The efficacy model must be from 'ModelEff' class (Effmodel slot) ## The DLE and efficayc samples must be from 'Samples' class (DLEsamples and Effsamples slot) Effmodel<-Effloglog(Eff=c(1.223,2.513),Effdose=c(25,300),nu=c(a=1,b=0.025),data=data,c=0) Effsamples <- mcmc(data, Effmodel, McmcOptions(burnin=100,step=2,samples=200)) ## Given a dose level 75, gain(dose=75,DLEmodel=DLEmodel,DLEsamples=DLEsamples,Effmodel=Effmodel,Effsamples=Effsamples) ##Obtain the gain value for a given dose, a pseudo DLE model, a DLE sample, ## the 'EffFlexi' efficacy model and an efficacy sample ##The DLE model must be from 'ModelTox' class (DLEmodel slot) emptydata<- DataDual(doseGrid=seq(25,300,25)) data<-emptydata DLEmodel<-LogisticIndepBeta(binDLE=c(1.05,1.8),DLEweights=c(3,3),DLEdose=c(25,300),data=data) DLEsamples <- mcmc(data, DLEmodel, McmcOptions(burnin=100,step=2,samples=200)) ##The efficacy model must be from 'EffFlexi' class (Effmodel slot) ## The DLE and efficayc samples must be from 'Samples' class (DLEsamples and Effsamples slot) EffFleximodel <- EffFlexi(Eff=c(1.223, 2.513),Effdose=c(25,300), sigma2=c(a=0.1,b=0.1),sigma2betaW=c(a=20,b=50),smooth="RW2",data=data) Effsamples <- mcmc(data, EffFleximodel, McmcOptions(burnin=100,step=2,samples=200)) ## Given a dose level 75, gain(dose=75,DLEmodel=DLEmodel,DLEsamples=DLEsamples,Effmodel=EffFleximodel,Effsamples=Effsamples) ##Obtain the gain value for a given dose, a pseudo DLE model and a pseudo efficacy model ## without samples ##The DLE model must be from 'ModelTox' class (DLEmodel slot) emptydata<- DataDual(doseGrid=seq(25,300,25),placebo=FALSE) data<-Data(doseGrid=seq(25,300,25),placebo=FALSE) DLEmodel<-LogisticIndepBeta(binDLE=c(1.05,1.8),DLEweights=c(3,3),DLEdose=c(25,300),data=data) ##The efficacy model must be from 'Effloglog' class (Effmodel slot) Effmodel<-Effloglog(Eff=c(1.223,2.513),Effdose=c(25,300),nu=c(a=1,b=0.025),data=emptydata,c=0) ## Given a dose level 75, gain(dose=75,DLEmodel=DLEmodel,Effmodel=Effmodel)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.