plrs.cb {plrs} | R Documentation |
Determine uniform confidence intervals for predicted values of a 'plrs' model.
plrs.cb(object, alpha=0.05, newcgh=NULL)
object |
An object of class |
alpha |
Significance level |
newcgh |
Vector of segmented values. Support for building CB. |
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.
An object of class plrs-class that contains CB information.
Gwenael G.R. Leday g.g.r.leday@vu.nl
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).
# 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)