In building a graphical model, accuracy in edge detection for the model structure is crucial for the quality of the model. We explored methods for improvement of false discovery rate(FDR) by devising an estimation procedure which is more data sensitive under some condition. The estimation is made by applying an EM method where the parameters include the density function under the null hypothesis (no edge) and the location parameters of the density functions under the alternative hypothesis (presence of edge). Our method is compared favorably with a most popular FDR tool in numerical experiments. We applied our method for analysing gene data of 800 genes and built a network of vector autoregressive model for the data.
AubertJ.Bar-HenA.DaudinJ.J. and RobinS., Determination of the differentially expressed genes in microarray experiments using local FDR, BMC Bioinformatics5 (2004), 125.
2.
BayS.D.ShragerJ.PohorilleA. and LangleyP., Revising regulatory networks: From expression data to linear causal models, J. Biomed. Informatics35 (2002), 298–497.
3.
BenjaminiY., Discovering the false discovery rate, J. Roy. Statist. Soc. B72(4) (2010), 405–416.
4.
BenjaminiY. and HochbergY., Controling the false discovery rate: A practical and powerful approach to multiple testing, J. Roy. Statist. Soc. B57 (1995), 289–300.
5.
BickelP.J. and DoksumK.A., in: Mathematical Statistics. 2nd ed. New-York, CRC Press (Talor and Francis Group), 2015.
6.
CraigonD.J.JamesN.OkyereJ.HigginsJ.JothamJ. and MayS., NASCArrays: A repository for microarray data generated by NASC’s transcriptomics service, Nucleic Acids Research32 (2004), D575–D577.
7.
DalmassoC.BröretP. and MoreauT., A simple procedure fir estimating the false discovery rate, Bioinformatics21 (2005), 660–668.
8.
DempsterA.P.LairdM. and RubinB., Maximum likelihood from incomplete data via the EM algorithm, Journal of the Royal Statistical Society B39 (1977), 1–38.
9.
DempsterA.P., The direct use of likelihood for significance testing, Statistics and Computing7 (1997), 247–252.
DunnetC.W., Multiple comparisons procedure for comparing several treatments with a control, J. Am. Statist. Assoc.50 (1955), 1096–1121.
12.
DunnetC.W., New tables for multiple comparisons with a control, Biometrics20 (1964), 482–491.
13.
EdwardsD., in: Introduction to Graphical Modelling, New York, Springer, 1995.
14.
EfronB., Robbins, empirical Bayes and microarrays, The Annals of Statistics31 (2003), 366–378.
15.
EfronB., Large-scale simultaneous hypothesis testing: The choice of a null hypothesis, Journal of the American Statistical Association99 (2004), 96–104.
16.
EfronB., Local False Discovery Rates, National Institute of Health grant 8R01 EB002784 and National Science Foundation grant DMS-0072360, Technical Report N. 2005-20B (234), 2005.
17.
EfronB., Correlation and large-scale simultaneous significance testing, J. Amer. Statist. Assoc.102 (2007), 93–103.
18.
EfronB. and TibshiraniR., Empirical Bayes methods and false discovery rates for microarrays, Genetic Epidemiology23 (2002), 70–86.
19.
EfronB.TibshiraniR.StoreyJ. and TusherV., Empirical Bayes analysis of a microarray experiment, J. Amer. Statist. Assoc.96 (2001), 1151–1160.
20.
FisherR.A., Frequency distribution of the values of the correlation coefficient in samples of an indefinitely large population, Biometrika10(4) (1915), 507–521.
21.
FisherR.A., On the ‘probable error’ of a coefficient of correlation deduced from a small sample, Metron1 (1921), 3–32.
22.
FriedmanJ., Regularized discriminant analysis, Journal of the American Statistical Association84 (1989), 165–175.
23.
HastieT.TibshiraniR. and FriedmanJ., The Elements of Statistical Learning, Springer, NY, 2001.
24.
HotellingH., New light on the correlation coefficient and its transforms, J. R. Statist. Soc. B15 (1953), 193–232.
25.
LaurizenS.L., Graphical Models, New York, Oxford University Press, 1996.
26.
LeeN.KimA.-K.ParkC. and KimS.-H., An improvement on local FDR analysis applied to functional MRI data, J. Neuroscience Methods267 (2016), 115–125.
27.
LeeN.ChoiH. and KimS.-H., Bayes shrinkage estimation for high-dimensional VAR models with scale mixture of normal distributions for noise, Computational Statistics and Data Analysis101 (2016), 250–276.
28.
McLachlanG. and PeelD., Finite mixture models, New York, John Wiley Sons, 2000.
29.
McLachlanG. and KrishnanT., The EM Algorithm and Extensions, New York, John Wiley Sons, 1997.
Opgen-RheinR. and StrimmerK., Learning causal networks from systems biology time course data: An effective model selection procedure for the vector autoregressive process, BMC Bioinformatics8(Suppl 2) (2007), S3.
32.
PaulsoE., On the comparison of several experimental categories with a control, Ann. Math. Statist.23 (1952), 239–246.
33.
PoundsS. and MorrisS.W., Estimaing the occurrence of false positives and false negatives in microarray studies by approximating and partitioning the empirical distribution of p-values, Bioinformatics19 (2003), 1236–1242.
34.
SchaferJ.Opgen-RheinR. and StrimmerK., GeneNet: Modeling and inferring gene networks, R Package Version 1.2.11, 2014.
35.
SchaferJ. and StrimmerK., An empirical Bayes approach to inferring large-scale gene association networks, Bioinformatics21 (2005), 754–764.
36.
SmithS.M.FultonD.C.ChiaT.ThorneycroftD.ChappleA.DunstanH.HyltonC.ZeemanS.C. and SmithA.M., Diurnal changes in the transcriptom encoding enzymes of starch metabolism provide evidence for both transcriptional and posttranscriptional regulation of starch metabolism in Arabidopsis leaves, Plant Physiology136 (2004), 2687–2699.
37.
StoreyJ.D., A direct approach to false discovery rates, J. R. Statist. Soc. B64 (2002), 479–498.
38.
StoreyJ.D. and TibshiraniR., Statistical significance fir genomewide studies, Proc. Natl. Acad. Sci.100 (2003), 9440–9445.
39.
StrimmerK., Fdrtool: A versatile R package for estimating local and tail area-based false discovery rates, Bioinformatics24 (2008), 1461–1462.
40.
SunW. and CaiT.T., Large-scale multiple testing under dependence, J. R. Statist. Soc. B71 (2009), 393–424.
41.
TohH. and HorimotoK., Inference of a genetic network by a combined approach of cluster analysis and graphical Gaussian modeling, Bioinformatics18 (2002), 287–297.
42.
WaddellP.J. and KishinoH., Cluster inferences methods and graphical models evaluated on NCI60 microarray gene expression data, Genome Informatics11 (2000), 129–140.
43.
WangJ.MyklebostO. and HovigE., MGraph: Graphical model for microarray data analysis, Bioinformatics19 (2003), 2210–2211.
44.
WhittakerJ., Graphical Models in Applied Multivariate Statistics, New York, John Wiley Sons, 1990.
45.
WuX.YeY. and SubramanianK.R., Interactive analysis of gene interactions using graphical Gaussian model, in: Proceedings of the ACM SIGKDD Workshop on Data Mining in Bioinformatics, Vol. 3, 2003, pp. 63–69.
46.
ZweigM.H. and CampbellG., Receiver-operating characteristic (ROC) plots: A fundamental evaluation tool in clinical medicine, Clinical Chemistry39 (1993), 561–577.
