Using Importance Sampling to Improve Simulation in Linkage Analysis
-
Lars Ängquist
und Ola Hössjer
In this article we describe and discuss implementation of a weighted simulation procedure, importance sampling, in the context of nonparametric linkage analysis. The objective is to estimate genome-wide p-values, i.e. the probability that the maximal linkage score exceeds given thresholds under the null hypothesis of no linkage. In order to reduce variance of the estimate for large thresholds, we simulate linkage scores under a distribution different from the null with an artificial disease locus positioned somewhere along the genome. To compensate for the fact that we simulate under the wrong distribution, the simulated scores are reweighted using a certain likelihood ratio. If the sampling distribution are properly chosen the variance of the corresponding estimate is reduced. This results in accurate genome-wide p-value estimates for a wide range of large thresholds with a substantially smaller cost adjusted relative efficiency with respect to standard unweighted simulation.We illustrate the performance of the method for several pedigree examples, discuss implementation including the amount of variance reduction and describe some possible generalizations.
©2011 Walter de Gruyter GmbH & Co. KG, Berlin/Boston
Artikel in diesem Heft
- Article
- Using Alpha Wisely: Improving Power to Detect Multiple QTL
- Relating HIV-1 Sequence Variation to Replication Capacity via Trees and Forests
- Linear Models and Empirical Bayes Methods for Assessing Differential Expression in Microarray Experiments
- Asymptotic Optimality of Likelihood-Based Cross-Validation
- Using Importance Sampling to Improve Simulation in Linkage Analysis
- Model-Based Assignment and Inference of Protein Backbone Nuclear Magnetic Resonances
- Error-Rate and Decision-Theoretic Methods of Multiple Testing: Which Genes Have High Objective Probabilities of Differential Expression?
- Evaluation of Multiple Models to Distinguish Closely Related Forms of Disease Using DNA Microarray Data: an Application to Multiple Myeloma
- Saturation and Quantization Reduction in Microarray Experiments using Two Scans at Different Sensitivities
- Combining Nearest Neighbor Classifiers Versus Cross-Validation Selection
- Multiple Testing. Part I. Single-Step Procedures for Control of General Type I Error Rates
- Multiple Testing. Part II. Step-Down Procedures for Control of the Family-Wise Error Rate
- Augmentation Procedures for Control of the Generalized Family-Wise Error Rate and Tail Probabilities for the Proportion of False Positives
- Calculating the Statistical Significance of Changes in Pathway Activity From Gene Expression Data
- A Family-Based Association Test for Repeatedly Measured Quantitative Traits Adjusting for Unknown Environmental and/or Polygenic Effects
- Deletion/Substitution/Addition Algorithm in Learning with Applications in Genomics
- Classifying Gene Expression Profiles from Pairwise mRNA Comparisons
- Hierarchical Bayesian Neural Network for Gene Expression Temporal Patterns
- A Mixed Model Approach to Identify Yeast Transcriptional Regulatory Motifs via Microarray Experiments
- Mammalian Genomes Ease Location of Human DNA Functional Segments but Not Their Description
- On the Dependence Structure of Sequence Alignment Scores Calculated with Multiple Scoring Matrices
- Increasing Power for Tests of Genetic Association in the Presence of Phenotype and/or Genotype Error by Use of Double-Sampling
- A Method for Evaluating the Impact of Individual Haplotypes on Disease Incidence in Molecular Epidemiology Studies
- Statistical Methods for Identifying Conserved Residues in Multiple Sequence Alignment
- MergeMaid: R Tools for Merging and Cross-Study Validation of Gene Expression Data
- Sparse Inverse of Covariance Matrix of QTL Effects with Incomplete Marker Data
- Maximum Likelihood for Genome Phylogeny on Gene Content
- Confidence Levels for the Comparison of Microarray Experiments
- PLS Dimension Reduction for Classification with Microarray Data
- Statistical Analysis of Genomic Tag Data
- Statistical Analysis of Adsorption Models for Oligonucleotide Microarrays
- Statistical Significance Threshold Criteria For Analysis of Microarray Gene Expression Data
- A Compendium to Ensure Computational Reproducibility in High-Dimensional Classification Tasks
- Validation and Discovery in Markov Models of Genetics Data
- Making Sense of High-Throughput Protein-Protein Interaction Data
- Reader's Reaction
- Reader Reaction
- Response to Foulkes and De Gruttola
- Software Communication
- BayesMendel: an R Environment for Mendelian Risk Prediction
- Letter to the Editor
- Concerns About Unreliable Data from Spotted cDNA Microarrays Due to Cross-Hybridization and Sequence Errors
Artikel in diesem Heft
- Article
- Using Alpha Wisely: Improving Power to Detect Multiple QTL
- Relating HIV-1 Sequence Variation to Replication Capacity via Trees and Forests
- Linear Models and Empirical Bayes Methods for Assessing Differential Expression in Microarray Experiments
- Asymptotic Optimality of Likelihood-Based Cross-Validation
- Using Importance Sampling to Improve Simulation in Linkage Analysis
- Model-Based Assignment and Inference of Protein Backbone Nuclear Magnetic Resonances
- Error-Rate and Decision-Theoretic Methods of Multiple Testing: Which Genes Have High Objective Probabilities of Differential Expression?
- Evaluation of Multiple Models to Distinguish Closely Related Forms of Disease Using DNA Microarray Data: an Application to Multiple Myeloma
- Saturation and Quantization Reduction in Microarray Experiments using Two Scans at Different Sensitivities
- Combining Nearest Neighbor Classifiers Versus Cross-Validation Selection
- Multiple Testing. Part I. Single-Step Procedures for Control of General Type I Error Rates
- Multiple Testing. Part II. Step-Down Procedures for Control of the Family-Wise Error Rate
- Augmentation Procedures for Control of the Generalized Family-Wise Error Rate and Tail Probabilities for the Proportion of False Positives
- Calculating the Statistical Significance of Changes in Pathway Activity From Gene Expression Data
- A Family-Based Association Test for Repeatedly Measured Quantitative Traits Adjusting for Unknown Environmental and/or Polygenic Effects
- Deletion/Substitution/Addition Algorithm in Learning with Applications in Genomics
- Classifying Gene Expression Profiles from Pairwise mRNA Comparisons
- Hierarchical Bayesian Neural Network for Gene Expression Temporal Patterns
- A Mixed Model Approach to Identify Yeast Transcriptional Regulatory Motifs via Microarray Experiments
- Mammalian Genomes Ease Location of Human DNA Functional Segments but Not Their Description
- On the Dependence Structure of Sequence Alignment Scores Calculated with Multiple Scoring Matrices
- Increasing Power for Tests of Genetic Association in the Presence of Phenotype and/or Genotype Error by Use of Double-Sampling
- A Method for Evaluating the Impact of Individual Haplotypes on Disease Incidence in Molecular Epidemiology Studies
- Statistical Methods for Identifying Conserved Residues in Multiple Sequence Alignment
- MergeMaid: R Tools for Merging and Cross-Study Validation of Gene Expression Data
- Sparse Inverse of Covariance Matrix of QTL Effects with Incomplete Marker Data
- Maximum Likelihood for Genome Phylogeny on Gene Content
- Confidence Levels for the Comparison of Microarray Experiments
- PLS Dimension Reduction for Classification with Microarray Data
- Statistical Analysis of Genomic Tag Data
- Statistical Analysis of Adsorption Models for Oligonucleotide Microarrays
- Statistical Significance Threshold Criteria For Analysis of Microarray Gene Expression Data
- A Compendium to Ensure Computational Reproducibility in High-Dimensional Classification Tasks
- Validation and Discovery in Markov Models of Genetics Data
- Making Sense of High-Throughput Protein-Protein Interaction Data
- Reader's Reaction
- Reader Reaction
- Response to Foulkes and De Gruttola
- Software Communication
- BayesMendel: an R Environment for Mendelian Risk Prediction
- Letter to the Editor
- Concerns About Unreliable Data from Spotted cDNA Microarrays Due to Cross-Hybridization and Sequence Errors