compareAIC {pulsedSilac} | R Documentation

## Calculates the probability of each model from a set of models.

### Description

For a given set of AIC from models, the probability of each model relative to the rest of the models of the set is calculated using the following formula:

### Usage

    compareAIC(...)

### Arguments

`...` | a `list` with the model metrics, the output from `modelTurnover` and `calculateAIC`.

### Value

a `list` with a matrix for each experiment condition. The matrix contains the probabilities of each model (columns) for each protein/peptide (rows).

### See Also

`calculateAIC`, `modelTurnover`.

### Examples

    data('wormsPE')
    wormsPE <- calculateIsotopeFraction(wormsPE, ratioAssay = 'ratio')

    modelList1 <- modelTurnover(x = wormsPE[1:10],
                               assayName = 'fraction',
                               formula = 'fraction ~ 1 - exp(-k*t)',
                               start = list(k = 0.02),
                               mode = 'protein',
                               robust = FALSE,
                               returnModel = TRUE)

    modelList1 <- calculateAIC(modelList1, smallSampleSize = TRUE)

    modelList2 <- modelTurnover(x = wormsPE[1:10],
                               assayName = 'fraction',
                               formula = 'fraction ~ 1 - exp(-k*t) + b',
                               start = list(k = 0.02, b = 0),
                               mode = 'protein',
                               robust = FALSE,
                               returnModel = TRUE)

    modelList2 <- calculateAIC(modelList2, smallSampleSize = TRUE)

    modelProbabilities <- compareAIC(modelList1, modelList2)

[Package pulsedSilac version 1.9.1 Index]