scatterCompareModels: Scatter plot of two conditions for a model metric.

In marcpaga/pulsedSilac: Analysis of pulsed-SILAC quantitative proteomics data

Description

Scatter plot of two conditions/replicates for a selected metric of a model. For example to compare turnover rates, model errors... facet_wrap.

Usage

scatterCompareModels(
  modelList,
  conditions,
  value = "param_values",
  returnDataFrame = FALSE
)

Arguments

modelList	A list containing all the model objects, this should be the output of modelTurnover with returnModel as TRUE.
conditions	A character of length 2 indicating which 2 conditions should be compared.
value	A character indicating which metric to plot. Check names(modelList) for available options. (Default = 'param_values')
returnDataFrame	A logical indicating if the data.frame used for the plot should be returned instead.

Value

A ggplot object or the data.frame that would be used instead in the plot.

Examples

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

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

scatterCompareModels(modelList = modelList,
                     conditions = c('OW40', 'OW450'),
                     value = 'param_values')

marcpaga/pulsedSilac documentation built on March 11, 2020, 8:49 p.m.