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

In plrs: Piecewise Linear Regression Splines (PLRS) for the association between DNA copy number and gene expression

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.

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.

Value

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

Author(s)

Gwenael G.R. Leday [e-mail]

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)

plrs documentation built on April 28, 2020, 6:09 p.m.