RXshrink (version 1.7)

plot.eff.biv: Plot method for eff.biv objects

Description

Plot Bivariate Confidence Ellipse(s) or a TRACE display for unrestricted ridge shrinkage coefficients. The default is to display two ellipses at confidence levels 0.95 and 0.50 on a single plot for the first 2 of p >= 2 coefficients. The projection of the eff.ridge() Path from the OLS solution (m == 0) to the Minimum MSE Risk solution (m = mStar) onto this bivariate-plane is also shown in "red".

Usage

# S3 method for eff.biv
plot(x, type = "ellip", ...)

Arguments

x

Output list object of class eff.ridge.

type

One of 2 options: "ellip" or "trace". The default option of "ellip" displays the Confidence Ellipse(s) specified by arguments to eff.biv().

…

Optional argument(s) passed on to plot().

Value

NULL

Examples

# NOT RUN {
  # Cost-Effectiveness inferences using Linear Models and GRR...
  
# }
# NOT RUN {
    # library(ICEinfer)
    data(sepsis)
    ndr <- ICEpref(sepsis$qalypres, sepsis$totcost, lambda=50000, beta=0.1)
    sndr <- data.frame(cbind(ndr, sepsis)) # ndr: non-linear diminishing returns...
    form4 <- ndr ~ icu + age + orgfails + apache
    usndra <- eff.aug(eff.ridge(form4, sndr)) # compare ndr of 2 Intensive Care Units...
    plot(efobj <- eff.biv(usndra, 2, 4))
    efobj          # implicit print...
    # plot(efobj, type = "tr")
  
# }