A function to model foldchange variance along several studies This function calculate the REM-summary fold-change — remodel • MetaVolcanoR

A function to model foldchange variance along several studies This function calculate the REM-summary fold-change

remodel(gene, foldchangecol, vcol)

Arguments

gene: named vector with foldchanges and variances <vector<
foldchangecol: the column name of the foldchange variable <string>
vcol: name of the fold change variance variable <string>

Value

data.frame with REM results for a gene

Examples

g <- data.frame('Symbol'="XGENE", 'Log2FC_1'=1.2, 'Log2FC'=0.8, 
                'vi_1'=0.01, 'vi_2'=0.1)
remodel(g, 'Log2FC', 'vi')
#>   signcon ntimes randomSummary randomCi.lb randomCi.ub      randomP   het_QE
#> 1       2      2        1.1125   0.7884014    1.436599 1.722858e-11 1.454545
#>   het_QEp   het_QM      het_QMp error
#> 1  0.2278 45.26286 1.722858e-11 FALSE