Description Usage Arguments Value Author(s) References Examples
This function will accept metabolomics sample data for a single metabolite and includes sample identification. It will bin corresponding samples by intensity and calculate a measure for each bin's error. The intensity per bin versus this measure of error will be plotted, and a best fit will be found and plotted that consists of a additive and a multiplicative part. For more details see the references.
1 2 | FoM.Calculate(data, cols.id = 1, col.value, ids.per.bin, quiet = FALSE,
repeats.per.id = -1, fit.type = 1, alpha.steps = 100)
|
data |
a matrix for which each row contains a sample's id and metabolomics data. |
cols.id |
indicate which column or columns in the matrix contain information pertaining to each sample's id, in which repeats of the same sample have the same id and can be binned together. |
col.value |
column in the data matrix that contains the metabolomics data to be plotted and fitted. |
ids.per.bin |
number of sample rows per bin. |
quiet |
should text output be given on the status of the calculations. |
repeats.per.id |
max number of repeats each sample may contain. If lower than the actual number repeats found, these extra repeats are ignored. |
fit.type |
which type of fitting should be used. Can be either "1" for use of the default R method "lm", or "2" for the more robust fitting method "rlm" from the 'MASS' package. |
alpha.steps |
how many points on the X axis in the final graph should be tried as a separation between the additive and the multiplicative fitting parts for finding the best fit. |
In addition to plotting the bins and showing the best fit, a list is returned containing the following values:
best.fit |
a vector containing the three parameters that define the best fit: alpha value, additive coefficient, multiplicative coefficient, and the residual of the fit. |
alphas |
a vector containing all the alpha values used for fitting. |
tot.res.ssq.per.alpha |
a vector that gives the alpha values' corresponding residuals per fit. |
ad.coeff.per.alpha |
a vector that gives the alpha values' corresponding additive coefficient per fit. |
mu.coeff.per.alpha |
a vector that gives the alpha values' corresponding multiplicative coefficient per fit. |
Marcel van Batenburg and Tim Dorscheidt
New Figures of Merit for Comprehensive Functional Genomics Data: The Metabolomics Case. Van Batenburg MF, Coulier L, van Eeuwijk F, Smilde AK, Westerhuis JA. Analytical Chemistry, Volume 83(9) (2011), pages 3267-3274
1 2 3 |
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.