LPS (version 1.0.4)

plot.LPS: Plot method for LPS objects

Description

This function plots the distributions of the LPS scores in each group for a fitted LPS object.

Usage

## S3 method for class 'LPS':
plot(x, y, method=c("Wright", "Radmacher", "exact"), xlim, yaxt = "n", xlab = "LPS",
    ylab, ...)

Arguments

x

An object of class "LPS", as returned by LPS.

y

Single character value defining y axis : "density" or (bayesian) "probability".

method

Single character value, the method to use for predictions. See predict.LPS.

xlim

To be passed to plot, see plot.default.

yaxt

To be passed to plot, see par.

xlab

To be passed to plot, see plot.default.

ylab

To be passed to plot, see plot.default.

...

Further arguments to be passed to plot or par.

See Also

LPS

Examples

# Data with features in columns
  data(rosenwald)
  group <- rosenwald.cli$group
  expr <- t(rosenwald.expr)
  
  # NA imputation (feature's mean to minimize impact)
  f <- function(x) { x[ is.na(x) ] <- round(mean(x, na.rm=TRUE), 3); x }
  expr <- apply(expr, 2, f)
  
  # Coefficients
  coeff <- LPS.coeff(data=expr, response=group)
  
  # 10 best features model
  m <- LPS(data=expr, coeff=coeff, response=group, k=10)
  
  
  # Distributions of scores in each group
  plot(m, "density")
  
  # Probability for each group along the score axis
  plot(m, "probability", yaxt="s")