Last data update: 2014.03.03

R: Calculates the Fisher's method p-value for each tested...
FM.chi.pvalueR Documentation

Calculates the Fisher's method p-value for each tested pathway

Description

Calculates the Fisher's method p-value for a set of p-values. It returns both the p-value and the test statistic value of the Fisher's product method.

Usage

FM.chi.pvalue(x)

Arguments

x

A vector of p-values. These p-values can be either gene or SNP p-values of a tested pathway

Value

FMstatistic

Fisher's product method test statistic

FMpvalue

Fisher's method p-value, computed using the exact distribution of the Fisher's method test statistic which is a Chi^2 distribution with degrees of freedom twice the size of vector x

References

Evangelou M, Rendon A, Ouwehand WH, Wernisch L, Dudbridge F (2012) Comparison of Methods for Competitive Tests of Pathway Analysis. PLoS ONE 7(7): e41018. doi:10.1371/journal.pone.0041018

See Also

pathways, snps.to.pathways

Examples

FM.chi.pvalue(x=c(0.05,0.1))

Results


R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(PAGWAS)
Loading required package: lars
Loaded lars 1.2

Loading required package: foreach
Loading required package: mnormt
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/PAGWAS/FM.chi.pvalue.Rd_%03d_medium.png", width=480, height=480)
> ### Name: FM.chi.pvalue
> ### Title: Calculates the Fisher's method p-value for each tested pathway
> ### Aliases: FM.chi.pvalue
> 
> ### ** Examples
> 
> FM.chi.pvalue(x=c(0.05,0.1))
$FMstatistic
[1] 10.59663

$FMpvalue
[1] 0.03149159

> 
> 
> 
> 
> dev.off()
null device 
          1 
>