print.emaxsim | R Documentation |
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.
## S3 method for class 'emaxsim'
print(x,
nprint = min(length(x$fitType), 20),
id = x$idmax,
digits = 3, ...)
x |
Output of |
nprint |
Number of simulations to print. If a vector of
length 2, |
id |
Output includes the stdBias for the dose with index |
digits |
Number of decimal digits to print for Z and p-values |
... |
Other print parameters (none currently implemented) |
Printed output returned as invisible matrix.
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.
Neal Thomas
emaxsim
, summary.emaxsim
,
plot.emaxsim
## 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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.