plrs.cb {plrs} | R Documentation

## 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.

### 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 g.g.r.leday@vu.nl

### 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)

[Package plrs version 1.26.0 Index]