normexp.fit {limma} | R Documentation

## Fit Normal+Exp Convolution Model to Observed Intensities

### Description

Fit normal+exponential convolution model to observed intensities. The normal part represents the background and the exponential represents the signal intensities. This function is called by `backgroundCorrect` and is not normally called directly by the user.

### Usage

    normexp.fit(foreground,background,trace=0)

### Arguments

`foreground` | numeric vector of foreground intensities
`background` | vector of background intensity values
`trace` | integer value passed to `optim`. If positive then tracing information on the progress of the optimization is given. Higher values give more information.

### Details

These functions use maximum likelihood estimation to fit a model to the foreground and background intensities. The model is a $normal(μ,σ^2)+exponential(α)$ convolution model for the foreground intensities in which $μ$ is a linear function of the background. `normexp.fit` uses a BFGS modified Newton-Raphson iteration maxmize the likelihood. `normexp.fit` assumes that $μ$ is a constant plus the observed background.

### Value

A list containing the components

`beta` | numeric vector of estimated regression coefficients
`sigma` | numeric scalar giving estimated value of $σ$
`alpha` | numeric scalar giving estimated value of $α$
`m2loglik` | numeric scalar giving minus twice the log-likelihood
`convergence` | integer code indicating successful convergence or otherwise of the optimization. See `optim`.

### Author(s)

Jeremy Silver and Gordon Smyth

### See Also

An overview of background correction functions is given in `04.Background`.

### Examples

    f <- c(2,3,1,10,3,20,5,6)
    b <- c(2,2,2,2,2,2,2,2)
    out <- normexp.fit(f,b)

[Package limma version 2.2.0 Index]