print.emaxsim {clinDR}

Print simulation output from emaxsim

Description

Prints key summary variables of Emax estimation peformance for each simulation. Can be used to identify simulated data sets yielding problems with common estimation methods.

Usage

## S3 method for class 'emaxsim'
print(x, 
        nprint = min(length(x$fitType), 20), 
        id = x$idmax, 
        digits = 3, ...)

Arguments

x

Output of emaxsim

nprint

Number of simulations to print. If a vector of length 2, nprint is the range of simulations to print.

id

Output includes the stdBias for the dose with index id vs placebo

digits

Number of decimal digits to print for Z and p-values

...

Other print parameters (none currently implemented)

Value

Printed output returned as invisible matrix.

Note

The stdBias printed is the difference between the estimated dose response at the dose with index id and its population value. The difference is divided by the SE of the estimator computed using the delta method.

Author(s)

Neal Thomas

See Also

emaxsim, summary.emaxsim, plot.emaxsim

Examples

## Not run: 
## emaxsim changes the random number seed
nsim<-50
idmax<-5
doselev<-c(0,5,25,50,100)
n<-c(78,81,81,81,77)

### population parameters for simulation
e0<-2.465375 
ed50<-67.481113 

dtarget<-100
diftarget<-9.032497
emax<-solveEmax(diftarget,dtarget,log(ed50),1,e0)

sdy<-7.967897
pop.parm<-c(log(ed50),emax,e0)    
meanlev<-emaxfun(doselev,pop.parm)  

###FixedMean is specialized constructor function for emaxsim
gen.parm<-FixedMean(n,doselev,meanlev,sdy)  
D1 <- emaxsim(nsim,gen.parm)

print(D1,c(31,50),digits=2,id=4)

print(D1,c(1,20))

D1  ### implicitly calls print with default parameter settings

## End(Not run)