R: Fisher's Exact g Test for Multiple (Genetic) Time Series
fisher.g.test R Documentation
Fisher's Exact g Test for Multiple (Genetic) Time Series
Description
fisher.g.test
calculates the p-value(s) according to Fisher's
exact g test for one or more time series. This test is useful to detect hidden
periodicities of unknown frequency in a data set. For an application to
microarray data see Wichert, Fokianos, and Strimmer (2004).
Usage
fisher.g.test(x, ...)
Arguments
x
vector or matrix with time series data (one time series per column).
...
arguments passed to periodogram
Details
Fisher (1929) devised an exact procedure to test the null hypothesis of Gaussian
white noise against the alternative of an added deterministic periodic component
of unspecified frequency. The basic idea behind the test is to reject the
null hypothesis if the periodogram contains a value significantly larger
than the average value (cf. Brockwell and Davis, 1991).
This test is useful in the context of microarray genetic time series
analysis as a gene selection method - see Wichert, Fokianos and Strimmer (2004)
for more details. Note that in the special case of a constant time series
the p-value returned by fisher.g.test
is
exactly 1 (i.e. the null hypothesis is not rejected).
Value
A vector of p-values (one for each time series). Multiple testing
may then be done using the the false discover rate approach
(function fdrtool
).
Author(s)
Konstantinos Fokianos (http://www.ucy.ac.cy/~fokianos/ ) and
Korbinian Strimmer (http://strimmerlab.org ).
References
Fisher, R.A. (1929). Tests of significance in harmonic analysis. Proc. Roy. Soc. A , 125 , 54–59.
Brockwell, P.J., and Davis, R.A. (1991). Time Series: Theory and Methods (2nd ed). Springer Verlag.
(the g-test is discussed in section 10.2).
Wichert, S., Fokianos, K., and Strimmer, K. (2004). Identifying periodically expressed transcripts
in microarray time series data. Bioinformatics 20 :5-20.
See Also
fdrtool
.
Examples
# load GeneCycle library
library("GeneCycle")
# load data set
data(caulobacter)
# how many samples and and how many genes?
dim(caulobacter)
# p-values from Fisher's g test
pval.caulobacter <- fisher.g.test(caulobacter)
pval.caulobacter
# compute Fdr and fdr values
fdr.out <- fdrtool(pval.caulobacter, statistic="pvalue")
# how many significant?
sum(fdr.out$qval < 0.05) # tail area-based Fdr
sum(fdr.out$lfdr < 0.2) # density-based local fdr
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(GeneCycle)
Loading required package: MASS
Loading required package: longitudinal
Loading required package: corpcor
Loading required package: fdrtool
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/GeneCycle/fisher.g.test.Rd_%03d_medium.png", width=480, height=480)
> ### Name: fisher.g.test
> ### Title: Fisher's Exact g Test for Multiple (Genetic) Time Series
> ### Aliases: fisher.g.test
> ### Keywords: htest
>
> ### ** Examples
>
>
> # load GeneCycle library
> library("GeneCycle")
>
> # load data set
> data(caulobacter)
>
> # how many samples and and how many genes?
> dim(caulobacter)
[1] 11 1444
>
> # p-values from Fisher's g test
> pval.caulobacter <- fisher.g.test(caulobacter)
> pval.caulobacter
[1] 6.957582e-01 6.572720e-04 3.877791e-04 3.409797e-04 3.808679e-03
[6] 4.153409e-03 1.033611e-03 1.806916e-04 3.492938e-01 7.050289e-01
[11] 7.780124e-02 7.650025e-01 8.918939e-01 6.649296e-01 2.865021e-01
[16] 9.292695e-01 1.919834e-01 2.166660e-03 8.005379e-01 8.951159e-01
[21] 6.200447e-01 9.659602e-01 2.458986e-01 6.380772e-01 5.083606e-01
[26] 8.201076e-01 2.373970e-01 7.290474e-01 4.540155e-01 2.125520e-01
[31] 1.377189e-01 5.476873e-01 4.402326e-01 8.459343e-01 2.900579e-01
[36] 1.990583e-02 5.511967e-01 3.620677e-01 8.208651e-02 4.125649e-01
[41] 8.274811e-01 9.306406e-01 4.932974e-01 6.265290e-01 3.981707e-01
[46] 3.629268e-01 6.651511e-01 4.142419e-01 2.754448e-01 3.448970e-01
[51] 7.857757e-01 4.431577e-01 9.677192e-01 7.312011e-02 1.506115e-01
[56] 4.591084e-01 7.548248e-01 1.780073e-02 2.524271e-01 9.146059e-02
[61] 8.749949e-02 6.158725e-01 6.634993e-01 5.154622e-02 1.528813e-01
[66] 6.909086e-01 6.993570e-01 4.265096e-01 5.204751e-01 4.868765e-01
[71] 6.406963e-01 3.104530e-01 9.702399e-01 9.016301e-01 2.198965e-01
[76] 9.040238e-02 9.695813e-01 5.733719e-01 1.184392e-02 1.086777e-01
[81] 6.909210e-02 9.349006e-01 2.896606e-01 6.247639e-01 6.732243e-01
[86] 7.469375e-01 3.673823e-01 8.049709e-01 7.100718e-01 9.403190e-01
[91] 8.265153e-01 8.644746e-01 2.197912e-02 7.905798e-01 3.501532e-01
[96] 1.118495e-01 4.600084e-01 8.255402e-01 7.700308e-01 1.501435e-01
[101] 1.312772e-01 8.588145e-02 2.666295e-01 6.406726e-01 3.616478e-01
[106] 1.172896e-02 8.566740e-01 9.602021e-01 9.858587e-01 9.268793e-01
[111] 9.165017e-01 3.588045e-01 2.338072e-01 1.554536e-01 8.868706e-01
[116] 6.680362e-01 8.220551e-01 3.622606e-01 1.729644e-01 7.473866e-01
[121] 2.547502e-01 3.150158e-01 8.029990e-01 7.218624e-01 4.988579e-01
[126] 8.920214e-01 9.914074e-01 5.229729e-01 9.348177e-01 4.150144e-01
[131] 4.167467e-01 1.345287e-01 3.742887e-01 1.185967e-04 3.472107e-04
[136] 2.667554e-02 4.563386e-01 9.931242e-01 4.555756e-01 5.272976e-01
[141] 2.970034e-01 7.997943e-01 1.458062e-02 5.993556e-01 3.598636e-01
[146] 3.087767e-01 5.488368e-01 2.088606e-02 2.983366e-03 5.645703e-01
[151] 4.013581e-01 1.475373e-01 3.631304e-01 2.960130e-01 7.458307e-01
[156] 5.426550e-01 3.416952e-02 2.517090e-01 6.155889e-02 4.925830e-01
[161] 1.643476e-03 1.856964e-01 1.268309e-01 7.667609e-01 1.875144e-01
[166] 4.573246e-01 8.311574e-01 5.658619e-01 5.614408e-01 8.675004e-01
[171] 5.273037e-01 1.866561e-01 2.681174e-01 4.251116e-01 8.902975e-01
[176] 8.807892e-01 3.876443e-01 2.594724e-01 4.193321e-01 7.738691e-01
[181] 3.219802e-01 4.009737e-01 1.682982e-01 7.316735e-01 5.633408e-01
[186] 8.885674e-01 3.003461e-01 6.717160e-01 9.021131e-01 7.861132e-01
[191] 2.439295e-02 5.423629e-01 7.982583e-01 4.236666e-01 2.233259e-01
[196] 7.621355e-02 3.359703e-01 8.311387e-02 5.880570e-01 4.345749e-01
[201] 8.321142e-01 5.766262e-01 6.402952e-01 7.894946e-01 9.695539e-01
[206] 4.138271e-01 2.821174e-02 4.760360e-01 2.985410e-03 9.865482e-01
[211] 1.280038e-01 6.870965e-03 3.423217e-01 6.882136e-01 2.568232e-01
[216] 4.385322e-01 8.677993e-01 8.643668e-01 1.946444e-02 2.987561e-01
[221] 3.245677e-03 9.439477e-01 3.371625e-01 6.143802e-01 6.377270e-01
[226] 6.390556e-01 4.672238e-01 8.078267e-01 4.155513e-01 1.729735e-01
[231] 2.866205e-01 4.468063e-01 4.661044e-01 2.627919e-02 7.450308e-01
[236] 2.273614e-02 9.825899e-01 7.858125e-01 5.303018e-02 7.220298e-01
[241] 8.986970e-02 7.056805e-01 2.710809e-01 7.927084e-01 9.900654e-02
[246] 6.928949e-01 2.941077e-01 7.003644e-01 5.356449e-01 7.219664e-01
[251] 1.021960e-04 7.384524e-01 5.259256e-01 4.640287e-03 2.701425e-03
[256] 2.793628e-03 4.634926e-03 1.325413e-01 4.227630e-02 9.692349e-01
[261] 7.589140e-01 9.175020e-01 8.945046e-01 1.141227e-02 2.378904e-01
[266] 2.041241e-02 4.127767e-03 8.166171e-01 9.822976e-01 4.050244e-01
[271] 8.927769e-02 1.928945e-01 8.979106e-01 2.327370e-02 8.032586e-01
[276] 2.281262e-01 7.619742e-01 1.936227e-01 4.453233e-02 6.147091e-01
[281] 1.716166e-01 6.999905e-02 9.991640e-04 2.963258e-01 9.473392e-01
[286] 3.136379e-01 4.238604e-03 8.421128e-01 4.965446e-02 6.473588e-01
[291] 9.362357e-01 3.836964e-01 1.208302e-01 8.977753e-01 3.031602e-01
[296] 8.311371e-01 3.749187e-01 7.282960e-02 3.556775e-01 5.053673e-01
[301] 3.885984e-01 3.815538e-01 6.454507e-02 3.114086e-02 3.288292e-01
[306] 9.162602e-01 8.355281e-01 4.570446e-01 9.537891e-03 7.469562e-02
[311] 5.900778e-01 9.852590e-01 7.889871e-01 1.454706e-01 9.362889e-01
[316] 8.670954e-01 6.054716e-01 4.472269e-01 8.596157e-01 3.823569e-01
[321] 9.441047e-01 9.592196e-01 9.930065e-01 7.386568e-01 1.623137e-01
[326] 5.596849e-01 7.278630e-01 1.266092e-01 3.675676e-01 9.643305e-01
[331] 2.500090e-01 3.104657e-03 4.223752e-02 1.601321e-02 1.281211e-01
[336] 2.171780e-02 7.258995e-01 1.871228e-01 2.255636e-01 1.218652e-01
[341] 1.518490e-01 8.796307e-01 4.276765e-01 2.012813e-01 2.588466e-02
[346] 3.917236e-01 8.382023e-02 5.520956e-03 1.594891e-01 9.903471e-01
[351] 7.455499e-01 2.133748e-01 7.308912e-01 1.312181e-04 4.024348e-03
[356] 3.848581e-01 8.202119e-01 9.960704e-01 4.645862e-01 2.426024e-03
[361] 3.411462e-01 2.353550e-01 9.999875e-01 8.191420e-01 2.100525e-01
[366] 7.375736e-01 1.011863e-03 3.131536e-02 3.041193e-01 1.201271e-01
[371] 2.379065e-01 9.381795e-01 8.325070e-01 2.674887e-01 1.866361e-01
[376] 9.629465e-01 1.556878e-01 3.866703e-01 9.244960e-05 1.906270e-03
[381] 2.503285e-01 2.478306e-01 2.774533e-01 2.199961e-01 7.110286e-01
[386] 5.513217e-01 2.312851e-01 7.328042e-01 3.111535e-01 4.267100e-01
[391] 8.178665e-01 5.029359e-01 7.187750e-01 5.830873e-01 3.180274e-01
[396] 8.205862e-01 8.768745e-01 1.631300e-01 9.297346e-01 6.866389e-01
[401] 1.741614e-01 5.270236e-01 3.732738e-01 4.161881e-01 8.454710e-01
[406] 9.231293e-01 3.804401e-01 8.657298e-01 9.728014e-01 9.315583e-01
[411] 9.850426e-01 4.498170e-01 1.515954e-01 3.012194e-01 9.957844e-02
[416] 9.288190e-01 5.421034e-03 5.281743e-01 9.803985e-02 4.546400e-01
[421] 4.719984e-01 5.415663e-02 3.881779e-01 6.406709e-01 9.364433e-02
[426] 9.915986e-01 9.466315e-01 1.606686e-01 6.388864e-02 2.443890e-01
[431] 7.959988e-01 3.113006e-01 4.425099e-01 5.969156e-03 1.516459e-01
[436] 8.507428e-01 6.439749e-01 2.801349e-01 4.110360e-01 2.522978e-01
[441] 8.590671e-01 7.805582e-01 4.331353e-01 2.047385e-01 7.136079e-01
[446] 9.915615e-01 2.503229e-01 8.391188e-01 1.024957e-01 2.345758e-03
[451] 4.203264e-02 4.919907e-01 2.860194e-02 6.856095e-01 6.359977e-01
[456] 4.183271e-01 1.568529e-01 6.175610e-01 9.803484e-01 6.289148e-03
[461] 1.087317e-02 1.646749e-01 6.276610e-01 4.947336e-01 3.223862e-01
[466] 5.568415e-01 1.345485e-02 3.595046e-01 1.284787e-01 7.326520e-01
[471] 3.995319e-02 3.397330e-04 6.204629e-02 3.471587e-02 7.342529e-04
[476] 2.085420e-01 9.270975e-01 8.650647e-01 1.599654e-01 7.240588e-01
[481] 4.611222e-01 9.941928e-01 8.702079e-01 2.357011e-01 1.000782e-01
[486] 9.775334e-02 3.250739e-01 8.474362e-01 2.120962e-01 5.178397e-02
[491] 1.613238e-02 2.914256e-01 4.438390e-01 1.516367e-01 4.605520e-01
[496] 5.809246e-01 5.441381e-01 7.100079e-02 5.240143e-01 2.267033e-01
[501] 6.301787e-01 9.989272e-01 5.849209e-01 7.497115e-01 7.727085e-01
[506] 4.175159e-01 8.307857e-01 7.253600e-01 9.434394e-01 5.963877e-01
[511] 4.543537e-02 5.128698e-01 9.552839e-01 3.456565e-02 3.834347e-01
[516] 9.102414e-01 6.435882e-02 8.893139e-01 5.548900e-01 3.602259e-01
[521] 3.142481e-01 7.708980e-01 8.553345e-01 1.339819e-01 9.184362e-01
[526] 4.702357e-01 2.477462e-01 4.931256e-04 4.778655e-03 1.075658e-01
[531] 7.097352e-01 8.319253e-01 6.013804e-01 4.704843e-01 7.509122e-01
[536] 8.006773e-01 6.570155e-01 6.427713e-01 4.106667e-02 6.000031e-01
[541] 3.816553e-01 7.083117e-01 1.862461e-02 6.607500e-01 1.987338e-01
[546] 2.335917e-01 2.435867e-02 2.494450e-01 2.733832e-01 4.620150e-01
[551] 1.299559e-01 9.319400e-01 9.581400e-01 2.639612e-01 5.691451e-02
[556] 6.769699e-02 4.133026e-04 7.688968e-01 3.267121e-01 3.525769e-01
[561] 9.949599e-01 1.620463e-01 2.921626e-01 2.208596e-01 7.252941e-01
[566] 4.023193e-01 9.415503e-01 5.927569e-01 4.512701e-01 6.809320e-01
[571] 9.093878e-01 2.774012e-01 5.143119e-01 4.616561e-01 3.488988e-01
[576] 4.492857e-01 3.513496e-01 9.844882e-01 5.131194e-01 4.077239e-01
[581] 6.675256e-01 8.009013e-01 9.847941e-03 5.399430e-01 2.570232e-01
[586] 6.591305e-01 8.264347e-01 2.125855e-01 3.707441e-01 9.853167e-01
[591] 5.514630e-02 1.076950e-01 1.484733e-01 4.777273e-02 2.174058e-02
[596] 6.763464e-01 5.357211e-01 2.997573e-01 8.678383e-01 2.076426e-02
[601] 5.040043e-02 6.792124e-02 2.129812e-01 2.082991e-01 7.441106e-01
[606] 6.669969e-01 9.936190e-01 2.673133e-01 1.059217e-01 4.306246e-01
[611] 7.291241e-01 8.154711e-01 3.806358e-03 1.874063e-01 2.739549e-01
[616] 5.105968e-01 1.581880e-01 7.161223e-02 2.554165e-01 8.479717e-01
[621] 9.626166e-01 1.053702e-04 6.524699e-01 6.943743e-02 8.744291e-01
[626] 7.532835e-01 4.445215e-01 1.209160e-01 6.946932e-01 9.364102e-01
[631] 1.362558e-01 6.335665e-01 8.249599e-01 9.379926e-01 6.775424e-02
[636] 5.122305e-01 1.037373e-01 8.802731e-02 4.500645e-01 7.653056e-01
[641] 2.231518e-02 2.520391e-03 3.218609e-01 9.480053e-01 3.968503e-01
[646] 6.591328e-01 8.368015e-01 6.963739e-01 1.635446e-02 1.758017e-01
[651] 1.238093e-02 4.136834e-01 6.513239e-01 5.943220e-01 6.795014e-01
[656] 6.899170e-01 1.788828e-01 1.556662e-01 2.699104e-01 8.150919e-01
[661] 1.720421e-01 6.499139e-01 1.313071e-02 2.448135e-04 4.560735e-04
[666] 8.990290e-01 2.868369e-01 9.515209e-03 1.953028e-02 6.453108e-01
[671] 9.857320e-01 9.718811e-01 5.812481e-01 9.814259e-01 9.987588e-01
[676] 6.444588e-01 4.201170e-01 2.530293e-04 9.722392e-01 2.699816e-01
[681] 4.821607e-03 5.351362e-01 9.474984e-01 3.050693e-01 2.593710e-01
[686] 9.816021e-01 9.969266e-01 7.046633e-01 2.427893e-01 1.260618e-01
[691] 3.278107e-02 4.056404e-01 3.383441e-01 9.616519e-01 3.231719e-02
[696] 7.951003e-01 9.065626e-01 9.445137e-01 9.367120e-01 9.597458e-01
[701] 2.797956e-01 6.804222e-01 1.370337e-01 5.286277e-01 2.250303e-01
[706] 8.697296e-01 6.062550e-02 2.550944e-02 4.651819e-03 1.690289e-04
[711] 2.579959e-01 3.492557e-01 4.389570e-01 9.201090e-01 2.027734e-01
[716] 9.241783e-02 1.211306e-01 6.918357e-01 3.186303e-03 2.023917e-01
[721] 1.659317e-01 6.602783e-01 4.616060e-01 3.419425e-01 3.014077e-01
[726] 1.954159e-01 7.633865e-03 4.141055e-01 9.718490e-01 9.246178e-02
[731] 5.721654e-03 1.179999e-03 8.443594e-01 4.997234e-01 5.871057e-01
[736] 9.074745e-03 9.906487e-01 9.369051e-01 5.646390e-01 9.605563e-01
[741] 4.938036e-01 9.828271e-01 7.625638e-03 6.220690e-03 1.092072e-02
[746] 9.529104e-01 3.939461e-01 6.915874e-01 1.940584e-03 1.211826e-02
[751] 3.458169e-01 1.816189e-01 8.858322e-02 2.916807e-01 8.119717e-02
[756] 4.300965e-01 8.382959e-01 5.135425e-01 7.126161e-01 1.125676e-01
[761] 4.283487e-01 2.726035e-01 1.734516e-01 1.366409e-01 1.402843e-01
[766] 6.186259e-02 9.923443e-01 4.106757e-01 1.233845e-03 2.880205e-02
[771] 8.455883e-01 6.189510e-01 7.869625e-01 2.658458e-01 1.774566e-01
[776] 5.523409e-01 4.920584e-01 7.194589e-01 9.383944e-01 9.438998e-01
[781] 5.703510e-01 2.308209e-01 5.890828e-01 2.162061e-01 6.739332e-01
[786] 1.690648e-01 5.566905e-01 3.195920e-01 6.974147e-01 2.974540e-02
[791] 1.248976e-02 3.959422e-01 9.609203e-01 4.528756e-01 7.468850e-01
[796] 3.233588e-01 2.059090e-01 8.950656e-02 2.186007e-01 1.082131e-02
[801] 9.686566e-02 5.609556e-01 2.535308e-01 3.289491e-01 3.684633e-01
[806] 4.661023e-01 4.519616e-01 2.606708e-01 6.077917e-01 9.348124e-01
[811] 7.123743e-01 9.125393e-01 8.091288e-01 5.761414e-01 1.933778e-01
[816] 2.816370e-01 3.568191e-02 1.805090e-01 5.523296e-01 5.813616e-01
[821] 5.498313e-01 6.442930e-01 2.290860e-01 5.353114e-01 2.605694e-01
[826] 1.588067e-01 7.802966e-01 4.812482e-01 5.924484e-01 7.172330e-01
[831] 1.442753e-01 3.906199e-01 8.875678e-01 6.567467e-01 2.252117e-01
[836] 1.799217e-01 4.511542e-01 7.566618e-01 3.623843e-01 7.104871e-01
[841] 4.729367e-01 6.802957e-01 8.803732e-01 8.815120e-01 2.155963e-02
[846] 7.692687e-01 6.814754e-01 1.698704e-03 7.007347e-01 2.348412e-02
[851] 2.604989e-02 2.772953e-02 3.547123e-01 4.708619e-01 7.728479e-02
[856] 6.763311e-02 1.275990e-02 2.030273e-01 3.313563e-02 3.327029e-01
[861] 9.897549e-01 4.532453e-01 2.177856e-01 3.637008e-01 6.545863e-01
[866] 1.673984e-04 3.655036e-01 4.251163e-01 1.917127e-01 1.917994e-02
[871] 5.667644e-01 9.574902e-05 3.570130e-02 3.719551e-01 1.671766e-01
[876] 1.860429e-02 1.165818e-01 6.254788e-02 2.437254e-01 1.084871e-02
[881] 4.156172e-01 8.447736e-02 8.123455e-01 9.806683e-01 5.701668e-01
[886] 3.409775e-01 9.020057e-01 9.944149e-01 1.204878e-01 9.804804e-01
[891] 4.234213e-01 3.768866e-02 2.507265e-03 2.215405e-01 1.664054e-01
[896] 3.287257e-01 1.897209e-01 4.374497e-01 3.457080e-04 2.034930e-05
[901] 4.650868e-01 9.766230e-01 4.107538e-01 4.567503e-01 4.652241e-01
[906] 8.667126e-02 9.522607e-01 7.294360e-01 2.682115e-05 6.456708e-03
[911] 5.337513e-01 7.174880e-01 9.642602e-01 3.507976e-01 6.591116e-01
[916] 8.490271e-02 3.379324e-01 5.187127e-01 7.063862e-01 9.983893e-01
[921] 7.074883e-01 8.250495e-01 5.039576e-01 3.335575e-01 1.368796e-01
[926] 2.777701e-01 4.783495e-01 6.222744e-01 1.346031e-01 5.914629e-01
[931] 9.454357e-01 6.988868e-01 3.831617e-01 5.836115e-01 9.881672e-01
[936] 4.846436e-01 4.352237e-01 3.384386e-01 2.801484e-01 5.076308e-01
[941] 2.673246e-01 1.207577e-01 7.804471e-03 5.928348e-01 6.321820e-01
[946] 6.978468e-01 6.093792e-01 4.535714e-01 8.992148e-01 1.136191e-03
[951] 8.162549e-03 2.045287e-01 3.464981e-01 8.306309e-01 7.726178e-01
[956] 1.684516e-01 9.996171e-01 7.590415e-01 9.602712e-01 2.274324e-01
[961] 2.521684e-02 1.229460e-02 3.520556e-01 3.173910e-01 6.577730e-01
[966] 3.472595e-01 3.757442e-01 7.851130e-01 9.299585e-01 6.032045e-01
[971] 1.924160e-01 3.501480e-01 3.344731e-01 8.751086e-01 7.910403e-03
[976] 3.161428e-01 9.078679e-04 6.733678e-01 9.372217e-01 4.958274e-01
[981] 5.046865e-01 2.407771e-01 4.253970e-01 8.106489e-02 8.376135e-01
[986] 6.981286e-01 3.656011e-01 1.006340e-01 2.136966e-03 6.397280e-02
[991] 3.816677e-02 3.793893e-01 1.312685e-01 6.768749e-01 5.879401e-01
[996] 5.946420e-01 2.200591e-01 1.917649e-01 2.046407e-01 2.688565e-01
[1001] 9.927945e-01 7.827219e-01 8.216932e-01 7.639140e-01 4.331255e-01
[1006] 8.536059e-01 9.854338e-03 5.315666e-01 2.339500e-01 1.689562e-01
[1011] 7.458535e-01 1.192675e-03 6.015991e-02 9.362453e-01 9.719672e-01
[1016] 3.407929e-01 6.035872e-01 8.779567e-01 5.418943e-01 3.296956e-01
[1021] 1.157567e-01 3.269984e-01 8.370632e-01 7.939335e-01 2.418308e-01
[1026] 4.429760e-01 8.868583e-01 2.996306e-01 9.445083e-01 3.132587e-02
[1031] 1.382059e-01 5.831313e-01 5.301108e-01 5.219493e-01 4.237601e-01
[1036] 9.123520e-01 9.353518e-01 9.996158e-01 6.501239e-01 1.532589e-01
[1041] 4.894652e-02 3.088557e-01 5.654534e-01 9.177715e-01 2.206476e-01
[1046] 6.046674e-01 9.320981e-01 7.321181e-01 9.333605e-01 1.312468e-03
[1051] 2.096150e-01 1.451235e-01 9.899671e-01 8.070048e-03 6.894543e-01
[1056] 1.240755e-02 8.955552e-01 1.383742e-01 1.402303e-01 5.447412e-03
[1061] 3.371312e-01 1.763426e-01 6.471613e-03 1.060774e-01 5.844932e-01
[1066] 7.436239e-01 9.984208e-01 6.785675e-01 8.963255e-01 1.913970e-01
[1071] 1.952804e-01 2.424176e-01 2.373751e-01 6.028071e-02 9.587592e-01
[1076] 2.404029e-01 5.205871e-01 3.591339e-01 5.008585e-01 8.238015e-02
[1081] 6.326477e-01 3.042068e-01 4.632327e-01 1.716829e-01 1.157641e-02
[1086] 5.763309e-01 8.710321e-01 8.431919e-01 1.905977e-01 1.171836e-01
[1091] 9.786325e-01 4.112923e-01 3.232382e-01 7.087520e-05 8.007405e-01
[1096] 1.468459e-01 9.646273e-01 4.857101e-01 9.223585e-02 6.834072e-01
[1101] 7.630949e-01 8.706647e-01 3.394838e-01 5.202456e-01 9.406846e-02
[1106] 3.601428e-01 2.774861e-02 7.452649e-01 3.591567e-01 6.355553e-01
[1111] 1.795070e-03 2.898256e-01 4.080636e-01 5.369867e-01 4.866465e-01
[1116] 1.799685e-02 4.571740e-01 9.103782e-01 7.775382e-01 3.503138e-01
[1121] 6.000318e-03 5.173214e-03 9.032226e-01 1.962604e-01 1.855214e-01
[1126] 3.689710e-01 5.909792e-01 5.057802e-01 5.113507e-01 7.409818e-01
[1131] 2.381879e-01 3.047865e-01 6.874998e-01 4.012675e-01 7.012597e-01
[1136] 8.578441e-01 2.611423e-01 7.798505e-02 1.918213e-02 3.294593e-01
[1141] 2.262600e-02 7.366264e-01 1.756298e-01 9.999396e-01 2.473014e-01
[1146] 2.716313e-02 6.011133e-01 9.219414e-03 7.160083e-04 8.744555e-01
[1151] 6.794000e-01 8.299704e-01 5.097563e-03 2.943119e-04 3.343434e-02
[1156] 9.257059e-02 6.064518e-01 9.259609e-01 6.565869e-01 5.610267e-01
[1161] 7.204515e-01 7.339376e-01 2.239772e-01 8.675129e-01 8.315208e-01
[1166] 8.288535e-01 1.135566e-01 2.901186e-01 5.196995e-02 3.191262e-04
[1171] 5.553692e-01 8.841875e-01 9.747981e-01 7.050407e-01 1.726083e-01
[1176] 9.139211e-01 1.531392e-03 3.084814e-01 2.351344e-01 3.428809e-01
[1181] 9.146032e-01 8.140290e-01 6.227279e-01 1.979116e-01 1.136669e-01
[1186] 4.896677e-01 8.994918e-01 6.355617e-01 2.439593e-01 7.818651e-01
[1191] 9.338323e-01 1.053343e-01 4.214656e-01 1.512838e-01 5.459431e-01
[1196] 3.294818e-01 2.981122e-01 5.259052e-01 2.634187e-01 4.003674e-01
[1201] 7.480913e-01 1.983166e-02 6.829663e-01 3.962474e-01 5.267314e-01
[1206] 4.600526e-01 9.891912e-02 5.297584e-01 2.971309e-02 5.974698e-01
[1211] 1.793098e-01 4.173478e-01 8.980057e-03 1.392042e-03 2.852270e-01
[1216] 4.509667e-01 7.959867e-01 4.457942e-01 7.013342e-01 6.792209e-01
[1221] 6.812030e-01 2.002745e-01 5.902089e-01 1.851516e-03 5.375253e-01
[1226] 2.442172e-01 2.087469e-01 4.975146e-01 1.666559e-01 5.975315e-01
[1231] 3.474704e-01 1.938561e-01 3.459574e-01 1.816959e-04 8.303703e-01
[1236] 7.429790e-01 2.333124e-01 1.938687e-02 2.271341e-01 4.170813e-01
[1241] 1.898019e-01 5.945289e-01 8.953346e-01 6.802865e-01 9.942728e-01
[1246] 8.850626e-01 2.022030e-02 8.182462e-01 2.339688e-01 2.704030e-02
[1251] 1.226598e-01 8.201587e-01 5.387029e-01 3.011918e-01 8.494211e-02
[1256] 1.559728e-01 9.177945e-01 8.018606e-01 4.501169e-01 3.544321e-01
[1261] 9.348873e-01 7.737535e-03 7.618657e-01 1.801546e-01 9.658815e-01
[1266] 4.418612e-01 4.367811e-01 8.790757e-01 1.425815e-01 1.730006e-01
[1271] 2.591742e-01 1.036088e-01 1.023781e-03 1.226695e-01 5.770939e-02
[1276] 4.759576e-01 3.889099e-01 6.462880e-01 1.941950e-01 5.735553e-01
[1281] 9.287005e-01 7.904743e-01 6.520779e-03 1.502478e-03 5.734128e-01
[1286] 6.951488e-02 5.481882e-01 2.277598e-01 8.254309e-01 6.291286e-01
[1291] 3.642264e-01 3.699068e-01 8.852641e-01 1.121432e-01 1.423599e-01
[1296] 8.734297e-01 8.447231e-01 2.491936e-01 9.448490e-01 7.271265e-01
[1301] 7.736852e-01 8.639188e-01 9.836097e-01 3.558641e-01 7.968994e-01
[1306] 4.374252e-01 9.695788e-01 3.775602e-01 8.468990e-01 9.795808e-01
[1311] 4.545902e-01 9.838040e-01 9.832613e-02 8.120928e-01 9.238701e-01
[1316] 7.437641e-01 2.941839e-02 5.553886e-01 9.999447e-01 2.266798e-01
[1321] 3.399635e-01 9.510012e-01 7.068365e-01 2.291417e-01 2.758068e-01
[1326] 1.592216e-01 4.715313e-01 7.522982e-01 5.599509e-01 2.848102e-01
[1331] 2.067470e-04 3.633856e-02 3.779812e-01 9.095188e-01 6.807155e-01
[1336] 5.981583e-01 3.885862e-03 2.533911e-01 2.288072e-01 4.535321e-01
[1341] 5.880058e-01 5.418773e-02 8.804519e-01 9.159724e-01 9.809949e-01
[1346] 2.563099e-02 3.741069e-01 3.745867e-01 6.929856e-01 5.389837e-01
[1351] 6.447368e-01 6.129741e-01 2.237609e-01 1.709164e-01 7.717357e-01
[1356] 8.496802e-03 1.965783e-01 1.975437e-01 4.953667e-03 4.948682e-01
[1361] 5.219853e-03 1.222875e-01 4.661485e-01 6.063303e-01 6.558735e-01
[1366] 9.610701e-01 6.024654e-01 6.926055e-01 9.816451e-01 3.883081e-01
[1371] 9.294555e-01 1.229579e-01 1.393033e-01 1.544243e-03 2.982837e-01
[1376] 1.362086e-01 1.076612e-01 5.315674e-03 4.620034e-03 3.069538e-01
[1381] 7.020441e-01 8.834083e-01 5.904733e-01 3.214284e-01 6.673249e-05
[1386] 9.926557e-01 2.109068e-04 7.285113e-01 5.496351e-01 8.463240e-01
[1391] 2.172265e-01 1.617920e-01 7.824169e-01 1.150808e-01 5.162716e-01
[1396] 2.833101e-01 1.594076e-01 1.755835e-01 6.752280e-01 1.874966e-01
[1401] 1.569779e-02 9.969462e-01 9.909602e-01 7.247149e-01 2.679907e-01
[1406] 5.863567e-02 2.216278e-01 4.268352e-01 4.400036e-01 7.094671e-03
[1411] 5.801707e-03 3.153221e-01 3.769857e-02 4.017904e-01 8.760763e-01
[1416] 8.170560e-01 5.442216e-01 8.252882e-01 8.964636e-01 9.393250e-01
[1421] 6.860821e-01 5.231100e-01 3.764019e-01 8.569157e-01 1.105108e-02
[1426] 2.780476e-02 7.239196e-01 4.412772e-01 3.224004e-01 2.278989e-01
[1431] 8.012767e-01 6.630372e-01 1.852143e-01 1.836574e-01 3.927131e-01
[1436] 4.204495e-01 1.195382e-01 3.495134e-02 8.055733e-01 8.717563e-01
[1441] 8.854306e-01 9.796768e-01 7.382750e-01 2.952665e-03
>
> # compute Fdr and fdr values
> fdr.out <- fdrtool(pval.caulobacter, statistic="pvalue")
Step 1... determine cutoff point
Step 2... estimate parameters of null distribution and eta0
Step 3... compute p-values and estimate empirical PDF/CDF
Step 4... compute q-values and local fdr
Step 5... prepare for plotting
>
> # how many significant?
> sum(fdr.out$qval < 0.05) # tail area-based Fdr
[1] 52
> sum(fdr.out$lfdr < 0.2) # density-based local fdr
[1] 94
>
>
>
>
>
>
> dev.off()
null device
1
>