R: Display Scatter for a possibly Transformed Bootstrap...
Display Scatter for a possibly Transformed Bootstrap Distribution of ICE Uncertainty
Assuming x is an output list object of class ICEuncrt, the default invocation of plot(x)
graphically displays the bootstrap distrib of ICE uncertainty currently stored in x. An
invocation of the form x10 <- plot(x, lfact=10) increases the value of x item lambda by a factor
of 10, displays that transformed bootstrap distribution, and stores it in object x10. When
the x item unit is cost, an invocation of the form xs <- plot(x, swu=TRUE) displays the bootstrap
distribution stored in x using effe units and stores the transformed distribution in object
## S3 method for class 'ICEuncrt'
plot(x, lfact = 1, swu = FALSE, alibi = FALSE, ...)
Output list object of class ICEuncrt.
Positive factor multiplying the stored value of x item lambda.
Logical value of TRUE or FALSE to control switching the stored value of
x item unit between the 2 possibilities, cost and effe.
Logical value of TRUE or FALSE to control scaling of axes. alibi = FALSE
produces the default alias graphic in which points in the bootstrap uncertainty scatter are
held fixed in space, and changes in lambda merely change the scaling (tick marks) along either
the horizontal axis of a cost unit display or else along the vertical axis of an effe unit
display. alibi = TRUE produces an alibi graphic in which the scaling (and range) is the
same along both axes, and changes in lambda cause the points in the bootstrap uncertainty
scatter to literally move either left or right in a cost unit display or else up or dowm in
an effe unit display.
Optional argument(s) passed on to plot().
After a single call to ICEuncrt() for an initial value of the Shadow Price of Health,
lambda, and an initial choice of display unit (cost or effe), multiple calls to
plot.ICEuncrt() are usually made. Alternative economic choices for lambda can be suggested
by making calls to ICEscale() with different values for lambda. Because the Bootstrap
Distribution of ICE Uncertainty is equivariant under changes in lambda, it is much faster to
transform an existing bootstrap distribution than to generate a new one for a different value
of lambda. The print.ICEuncrt() and plot.ICEuncrt() functions thus have 2 special parameters,
lfact and swa, that can change lambda and switch the units of measurement, respectively, without
actually regenerating the bootstrap distribution via a call to ICEuncrt().
Object of class ICEuncrt containing a possibly TRANSFORMED output list with items:
Saved value of the name of the data.frame in the original call to ICEuncrt().
Possibly changed, positive value of lfact * (x item lambda).
Possibly switched value of x item unit, cost or effe.
Saved integer value for number of bootstrap replications input to ICEuncrt.
Saved name of the treatment indicator within the df data.frame.
Saved name of the treatment effectiveness variable within the df data.frame.
Saved name of the treatment cost variable within the df data.frame.
Saved value of the sorted 3-variable (trtm,effe,cost) data.frame.
Observed value of (DeltaEffe, DeltaCost) when each patient is included exactly once.
R x 2 matrix of values of (DeltaEffe, DeltaCost) computed by transformation.
Saved value of the seed used to start pseudo random number generation.
Bob Obenchain <firstname.lastname@example.org>
Obenchain RL. Issues and algorithms in cost-effectiveness inference. Biopharmaceutical
Reports 1997; 5(2): 1-7. Washington, DC: American Statistical Association.
Obenchain RL. ICEplane: a windows application for incremental cost-effectiveness (ICE)
statistical inference. Copyright (c) Pharmaceutical Research and Manufacturers of America
(PhRMA.) http://members.iquest.net/~softrx/ 1997–2007.
Obenchain RL. Resampling and multiplicity in cost-effectiveness inference. Journal of
Biopharmaceutical Statistics 1999; 9(4): 563–582.
Cook JR, Heyse JF. Use of an angular transformation for ratio estimation in cost-effectiveness
analysis. Statistics in Medicine 2000; 19: 2989-3003.
ICEuncrt, ICEscale and ICEwedge.
# Transformation of a bootstrap distribution is fast.
dpuncs <- plot(dpunc, swu=TRUE)
R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
Loading required package: lattice
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/ICEinfer/plot.ICEuncrt.Rd_%03d_medium.png", width=480, height=480)
> ### Name: plot.ICEuncrt
> ### Title: Display Scatter for a possibly Transformed Bootstrap
> ### Distribution of ICE Uncertainty
> ### Aliases: plot.ICEuncrt
> ### Keywords: methods hplot
> ### ** Examples
Incremental Cost-Effectiveness (ICE) Bivariate Bootstrap Uncertainty
Shadow Price = Lambda = 0.26
Bootstrap Replications, R = 25000
Effectiveness variable Name = idb
Cost variable Name = ru
Treatment factor Name = dulx
New treatment level is = 1 and Standard level is = 0
Cost and Effe Differences are both expressed in cost units
Observed Treatment Diff = 1.6
Mean Bootstrap Trtm Diff = 1.576
Observed Cost Difference = -2.899
Mean Bootstrap Cost Diff = -2.915
> # Transformation of a bootstrap distribution is fast.
> dpuncs <- plot(dpunc, swu=TRUE)
Created & Maintained by Osamu Ogasawara (email@example.com) and