Description Usage Arguments Details Value See Also Examples

`LimmaFn`

uses functions in `limma`

package to easily compute the moderated t-statistics and p-values from differential gene expression tests comparing between different phenotypes even when sample size is small.

1 | ```
LimmaFn(pheno.v, data.m)
``` |

`pheno` |
A vector of sample phenotypes. Sample phenotype in a scientific research could be treatment/control, normal/cancer or smoker/non-smoker. Different phenotypes could each be encoded as 0/1 when inputting to |

`data` |
A matrix, the unnormalized gene expression dataset, should be a numeric matrix, with rows referring to genes and columns to samples. In this matrix you could use any type of gene IDs, like Entrez ID, Ensembl ID, HUGO gene symbol... But make sure to use the same gene annotation through out your analysis. |

This function computes the moderated t-statistic for users using empirical Bayes method, it is especially useful when the sample size is too small to perform parametric tests.

Given a gene expression data matrix and a vector indicating sample phenotype, `LimmaFn`

first fits a linear model using `lmFit`

, then it refits the model and do comparisons between any two different phenotypes with `contrasts.fit`

, finally it estimates moderated t-statistics for each comparison from the fitted model using empirical Bayes method (`eBayes`

) and output the result from the `topTable`

function.

Note that doing the `contrasts.fit`

step will not make a difference if you do comparison between two different sample status (treatment/control). However, When there are more than two sample status in your data set, this step will do comparison between every two status. And resulted summary tables will be stored in a list.

A table with rows for all genes (ranked by significance) and columns of log2 fold-change, average expression, moderated t-statistic, p-value, adjusted p-value (Benjaminiâ€“Hochberg procedure). The table is the output of `topTable`

function.

`lmFit`

for fitting a linear model, `contrasts.fit`

for refitting, `eBayes`

for Bayes method, `topTable`

for the output table.

1 2 3 4 5 6 |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.