Print simulation output from emaxsim
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
## 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)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.