RXshrink (version 1.4.3)

plot.unr.biv: Plot method for unr.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 unr.ridge() Path from the OLS solution (m == 0) to the Minimum MSE Risk solution (m == mUnr) onto this bivariate-plane is also shown in "red".

Usage

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

Arguments

x

Output list object of class unr.ridge.

type

One of 2 options: "ellip" or "trace". The default option of "ellip" displays the Confidence Ellipse(s) specified by arguments to unr.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 <- unr.aug(unr.ridge(form4, sndr)) # compare ndr of 2 Intensive Care Units...
    plot(urobj <- unr.biv(usndra, 2, 4))
    urobj               # implicit print...
    # plot(urobj, type = "tr")
  
# }