# input the dulxparx data of Obenchain et al. (2000).
# Effectiveness = idb, Cost = ru, trtm = dulx where
# dulx = 1 ==> Duloxetine and dulx = 0 ==> Paroxetine
# Display of Lambda => Shadow Price Summary Statistics...
ICEscale(dulxparx, dulx, idb, ru)
ICEscale(dulxparx, dulx, idb, ru, lambda=0.26)

# Bootstrap ICE Uncertainty calculations with R=25000 can be time consuming...
dpunc <- ICEuncrt(dulxparx, dulx, idb, ru, lambda=0.26, R=5000)

# Display the Bootstrap ICE Uncertainty Distribution...

dpwdg <- ICEwedge(dpunc)

# Compute VAGR Acceptability and ALICE Curves...
dpacc <- ICEalice(dpwdg)
plot(dpacc, show="VAGR")
plot(dpacc, show="Alice")

# Color Interior of Confidence Wedge with LINEAR Economic Preferences...
dpcol <- ICEcolor(dpwdg, gamma=1)
plot(dpcol, show="RBOW")
plot(dpcol, show="Hist")

# Increase Lambda and Recolor Confidence Wedge with NON-Linear Preferences...
dpcol <- ICEcolor(dpwdg, lfact=10)
plot(dpcol, show="RBOW")
plot(dpcol, show="Hist")

Try the ICEinfer package in your browser

Any scripts or data that you put into this service are public.

ICEinfer documentation built on Oct. 23, 2020, 8:31 p.m.