plrs (version 1.12.0)

plrs.cb: Uniform confidence bands (CB) for plrs models

Description

Determine uniform confidence intervals for predicted values of a 'plrs' model.

Usage

plrs.cb(object, alpha=0.05, newcgh=NULL)

Arguments

object

An object of class plrs-class.

alpha

Significance level

newcgh

Vector of segmented values. Support for building CB.

Value

An object of class plrs-class that contains CB information.

Details

The input object of class plrs-class has to result from function plrs.test.

The problem of finding (at a given x) a confidence interval for the mean response is expressed as a semi-definite optimization problem and solved using function csdp of package Rcsdp.

References

Leday GGR, Van der Vaart AW, Van Wieringen WN, Van de Wiel MA. Modeling association between DNA copy number and gene expression with constrained piecewise linear regression splines. Accepted for publication. Ann Appl Stat. (2012).

See Also

plrs.test

Examples

# Simulate data
sim <- plrs.sim(n=80, states=4, sigma=0.5)

# Fit a model 
model <- plrs(expr=sim$expr, cghseg=sim$seg, cghcall=sim$cal)

# Confidence bands
model <- plrs.test(model)
model <- plrs.cb(model, alpha=0.05)
plot(model)