Last data update: 2014.03.03

R: Display Scatter for a possibly Transformed Bootstrap...
plot.ICEuncrtR Documentation

Display Scatter for a possibly Transformed Bootstrap Distribution of ICE Uncertainty

Description

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 xs.

Usage

 ## S3 method for class 'ICEuncrt'
plot(x, lfact = 1, swu = FALSE, alibi = FALSE, ...)

Arguments

x

Output list object of class ICEuncrt.

lfact

Positive factor multiplying the stored value of x item lambda.

swu

Logical value of TRUE or FALSE to control switching the stored value of x item unit between the 2 possibilities, cost and effe.

alibi

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().

Details

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().

Value

Object of class ICEuncrt containing a possibly TRANSFORMED output list with items:

df

Saved value of the name of the data.frame in the original call to ICEuncrt().

lambda

Possibly changed, positive value of lfact * (x item lambda).

unit

Possibly switched value of x item unit, cost or effe.

R

Saved integer value for number of bootstrap replications input to ICEuncrt.

trtm

Saved name of the treatment indicator within the df data.frame.

xeffe

Saved name of the treatment effectiveness variable within the df data.frame.

ycost

Saved name of the treatment cost variable within the df data.frame.

effcst

Saved value of the sorted 3-variable (trtm,effe,cost) data.frame.

t1

Observed value of (DeltaEffe, DeltaCost) when each patient is included exactly once.

tb

R x 2 matrix of values of (DeltaEffe, DeltaCost) computed by transformation.

seed

Saved value of the seed used to start pseudo random number generation.

Author(s)

Bob Obenchain <wizbob@att.net>

References

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.

See Also

ICEuncrt, ICEscale and ICEwedge.

Examples

  data(dpunc)
  dpunc
  # Transformation of a bootstrap distribution is fast.
  dpuncs <- plot(dpunc, swu=TRUE)

Results


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.

> library(ICEinfer)
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
> 
>   data(dpunc)
>   dpunc

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)
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>