Description Usage Arguments Details Value Author(s) References Examples
The function returns the list of features in common using the two suggested rules hmax and h2 (Bayesian model) and additional ones defined by the user.
1 | extractFeatures.R(output.ratio, output.bay, feat.names, h = NULL)
|
output.ratio |
|
output.bay |
|
feat.names |
|
h |
|
To select a list of interesting features from the Bayesian model we suggest two decision rules in the paper: 1) the maximum of Median(R(h)) only for the subset of credibility intervals which do not include 1; 2) the largest threshold h for which the ratio R(h) is bigger than 2.
The first one is pointing out the strongest deviation from independence, whilst the second is the largest threshold where the number of features called in common at least doubles the number of features in common under independence. The user can define additional thresholds of interest and obtain the list of associated features.
The function returns an object of the class list. Each element is a matrix where the first column contains the name of the features while the other columns contain the p-values* from the experiments. It also saves a .csv file with the same information.
* instead of the p-values any other measure used to rank the features in the experiments can be used
max |
The list of features of interest selected on the basis of the threshold associated to R(hmax) |
rule2 |
The list of features of interest selected on the basis of the threshold associated to R(h2) |
User |
The list of features of interest selected on the basis of the additional thresholds selected by the user |
Alberto Cassese, Marta Blangiardo
1. M.Blangiardo and S.Richardson (2007) Statistical tools for synthesizing lists of differentially expressed features in related experiments, Genome Biology, 8, R54
1 2 3 4 5 6 7 | data = simulation(n=500,GammaA=1,GammaB=1,r1=0.5,r2=0.8,
DEfirst=300,DEsecond=200,DEcommon=100)
Th<- ratio(data=data$Pval)
Rh<- baymod(iter=100,output.ratio=Th)
feat.names = data$names
feat.lists <- extractFeatures.R(output.ratio=Th,output.bay=Rh,
feat.names=feat.names,h=NULL)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.