Preparing the Final Data Set
|
|
- Donna Daniels
- 5 years ago
- Views:
Transcription
1 Preparing the Final Data Set Kevin R. Coombes 17 March 2011 Contents 1 Executive Summary Introduction Aims/Objectives Methods Description of Data Statistical Methods Results Conclusions Details Load All Segment Information Appendix 5 1 Executive Summary 1.1 Introduction This report describes the analysis of a data set from Lynn Barron, a member of the laboratory of Lynne V. Abruzzo. This dataset was acquired using Illumina 610K SNP chips. The main goal of the study is to identify genetic abnormalities that are associated with clinical outcome (including overall survival and time-to-treatment). This is the thirteenth report in a series of related reports Aims/Objectives The objective of this report is to use the segment calls (as produced by Report 6) to produce a final data set containing the copy number and LOH status of each segment by patient and by chromosome. 1
2 13-finalDataSet Methods Description of Data The dataset contains measurements on 176 previously untreated patients with CLL. Extensive clinical followup is available Statistical Methods Raw data were processed in BeadStudio to yield genotype calls, log R ratios (LRR), and B allele frequencies (BAF) for each SNP in each patient. Since the study does not include matched normal DNA, the computations were performed relative to the pool of 120 HapMap samples run by Illumina. In Report 2, we applied the circular binary segmentation (CBS) algorithm to the intensity (log R ratio; LRR) data for each patient and each chromosome. CBS was first described by Olshan et al. [Biostatistics 2004; 23:657 63]; we use the implementation of CBS from the R package DNAcopy. In Report 3, we computed the odds ratio for LOH versus no LOH in windows of width 40 along each chromosome. In Report 4, we applied the CBS algorithm to transformed B allele frequency (BAF) values on each chromosome of each patient sample. In all three of those reports, we saved the segmentation results in per-patient files. In Report 6, we pooled the segment data from the different algorithms for each patient. We also computed summary statistics along each resulting segment, including the LRR mean and standard deviation, a summary of the genotpes for the SNPs in the region, and the best fit for modeling the BAF as a mixture of multiple components. In Reports 7 through 10, we assigned a meaningful call or interpretation to each segment. In this report, we combine all of these interpretations into the final data set. 1.3 Results We generate A COOL FILE Conclusions It would be nice to have some... 2 Details 2.1 Load All Segment Information We start by loading the segmentation data produced in Report 6. > load("rewind.rda") > shortnames <- sort(levels(rewind$samid)) We remove the sampels that failed QC or that should not be combined into the final data set.
3 13-finalDataSet 3 > removeme <- which(shortnames %in% c(.failedqc,.do.not.combine)) > if (length(removeme) > 0) { + combinenames <- shortnames[-removeme] else { + combinenames <- shortnames > save(combinenames, file = "combinenames.rda") > source("00rnw/snp-utils.r") Now we put the code in place to iterate over all chromosomes. In order to do this, we have to load the data for each chromosome (for at least one sample) in order to extract the positions of the SNP probes on the Illumina array. > curchr <- loadsnpdata(.egsample, chrname) > posn <- curchr$position Next, we read the positions and use the unique values from the segment data to indicate the start locations of pooled segments. > data <- rewind[rewind$chrom == chrname, ] > myholder <- data.frame(chrom = chrname, loc.start = sort(unique(data$loc.start))) Next, we compute the probe locations of the ends of the segments. > padded <- c(myholder$loc.start, 1 + posn[length(posn)]) > temp <- sapply(padded, function(i) { + n <- sum(posn < i) + ifelse(n < 0, NA, posn[n]) ) > myholder$loc.end <- temp[2:length(temp)] > rm(temp, padded) Then we count the number of SNP markers contained in each segment. > myholder$num.mark <- apply(myholder, 1, function(arow) { + sum((posn >= as.numeric(arow[2])) & (posn <= as.numeric(arow[3]))) ) The final step is to expand the calls for each patient from their own segments to the list of pooled segments > lev <- levels(rewind$call) > for (cid in combinenames) { + calls <- sapply(myholder$loc.start, function(loc, mydata) { + n <- sum(mydata$loc.start <= loc)
4 13-finalDataSet 4 + if (n < 1) { + warning("train") + return("nocall") + as.character(mydata[n, "Call"]), mydata = data[data$samid == cid, ]) + myholder[, cid] <- factor(calls, levels = lev) > if (!file.exists("final")) dir.create("final") Here is the main loop that puts the data set together. > memory.limit(2048) [1] > for (chrname in c(1:22, "X")) { + cat(paste("working on chromosome", chrname, "\n"), file = stderr()) + curchr <- loadsnpdata(.egsample, chrname) + posn <- curchr$position + data <- rewind[rewind$chrom == chrname, ] + myholder <- data.frame(chrom = chrname, loc.start = sort(unique(data$loc.start))) + padded <- c(myholder$loc.start, 1 + posn[length(posn)]) + temp <- sapply(padded, function(i) { + n <- sum(posn < i) + ifelse(n < 0, NA, posn[n]) ) + myholder$loc.end <- temp[2:length(temp)] + rm(temp, padded) + myholder$num.mark <- apply(myholder, 1, function(arow) { + sum((posn >= as.numeric(arow[2])) & (posn <= as.numeric(arow[3]))) ) + lev <- levels(rewind$call) + for (cid in combinenames) { + calls <- sapply(myholder$loc.start, function(loc, mydata) { + n <- sum(mydata$loc.start <= loc) + if (n < 1) { + warning("train") + return("nocall") + as.character(mydata[n, "Call"]), mydata = data[data$samid == cid, ]) + myholder[, cid] <- factor(calls, levels = lev)
5 13-finalDataSet 5 + write.table(myholder, file = file.path("final", paste("chr", chrname, + ".tsv", sep = "")), row.names = FALSE, col.names = TRUE, sep = "\t") + if (exists("finaldata")) { + finaldata <- rbind(finaldata, myholder) + else { + finaldata <- myholder > save(finaldata, file = "finaldata.rda") > write.table(finaldata, file = "finaldata.tsv", row.names = FALSE, col.names = TRUE, + sep = "\t") 3 Appendix This analysis was run in the following directory: > getwd() [1] "c:/mystuff/snp-cll/aa" Note that \\mdadqsfs02 is the standard insititutional location for storing data and analyses; N: is the name given to that location on this machine. This analysis was run in the following software environment: > sessioninfo() R version ( ) Platform: x86_64-pc-mingw32/x64 (64-bit) locale: [1] LC_COLLATE=English_United States.1252 LC_CTYPE=English_United States.1252 [3] LC_MONETARY=English_United States.1252 LC_NUMERIC=C [5] LC_TIME=English_United States.1252 attached base packages: [1] stats graphics grdevices utils datasets methods base loaded via a namespace (and not attached): [1] RColorBrewer_1.0-2
Reorganizing the data by sample
Reorganizing the data by sample Kevin R. Coombes 17 March 2011 Contents 1 Executive Summary 1 1.1 Introduction......................................... 1 1.1.1 Aims/Objectives..................................
More informationPreliminary Figures for Renormalizing Illumina SNP Cell Line Data
Preliminary Figures for Renormalizing Illumina SNP Cell Line Data Kevin R. Coombes 17 March 2011 Contents 1 Executive Summary 1 1.1 Introduction......................................... 1 1.1.1 Aims/Objectives..................................
More informationReorganizing the data by sample
Reorganizing the data by sample Kevin R. Coombes 23 March 2011 Contents 1 Executive Summary 1 1.1 Introduction......................................... 1 1.1.1 Aims/Objectives..................................
More informationPlotting Segment Calls From SNP Assay
Plotting Segment Calls From SNP Assay Kevin R. Coombes 17 March 2011 Contents 1 Executive Summary 1 1.1 Introduction......................................... 1 1.1.1 Aims/Objectives..................................
More informationMain Results. Kevin R, Coombes. 10 September 2011
Main Results Kevin R, Coombes 10 September 2011 Contents 1 Executive Summary 1 1.1 Introduction......................................... 1 1.1.1 Aims/Objectives.................................. 1 1.2
More informationPreprocessing and Genotyping Illumina Arrays for Copy Number Analysis
Preprocessing and Genotyping Illumina Arrays for Copy Number Analysis Rob Scharpf September 18, 2012 Abstract This vignette illustrates the steps required prior to copy number analysis for Infinium platforms.
More informationHow to use CNTools. Overview. Algorithms. Jianhua Zhang. April 14, 2011
How to use CNTools Jianhua Zhang April 14, 2011 Overview Studies have shown that genomic alterations measured as DNA copy number variations invariably occur across chromosomal regions that span over several
More informationRAPIDR. Kitty Lo. November 20, Intended use of RAPIDR 1. 2 Create binned counts file from BAMs Masking... 1
RAPIDR Kitty Lo November 20, 2014 Contents 1 Intended use of RAPIDR 1 2 Create binned counts file from BAMs 1 2.1 Masking.................................................... 1 3 Build the reference 2 3.1
More informationPackage RLMM. March 7, 2019
Version 1.44.0 Date 2005-09-02 Package RLMM March 7, 2019 Title A Genotype Calling Algorithm for Affymetrix SNP Arrays Author Nusrat Rabbee , Gary Wong
More informationPooling Segments. Kevin R. Coombes. 17 March 2011
Pooling Segments Kevin R. Coombes 17 March 2011 Contents 1 Executive Summary 2 1.1 Introduction......................................... 2 1.1.1 Aims/Objectives.................................. 2 1.2
More informationPackage icnv. R topics documented: March 8, Title Integrated Copy Number Variation detection Version Author Zilu Zhou, Nancy Zhang
Title Integrated Copy Number Variation detection Version 1.2.1 Author Zilu Zhou, Nancy Zhang Package icnv March 8, 2019 Maintainer Zilu Zhou Integrative copy number variation
More informationPackage lodgwas. R topics documented: November 30, Type Package
Type Package Package lodgwas November 30, 2015 Title Genome-Wide Association Analysis of a Biomarker Accounting for Limit of Detection Version 1.0-7 Date 2015-11-10 Author Ahmad Vaez, Ilja M. Nolte, Peter
More informationgenocn: integrated studies of copy number and genotype
genocn: integrated studies of copy number and genotype Sun, W., Wright, F., Tang, Z., Nordgard, S.H., Van Loo, P., Yu, T., Kristensen, V., Perou, C. February 22, 2010 1 Overview > library(genocn) This
More informationPackage GEM. R topics documented: January 31, Type Package
Type Package Package GEM January 31, 2018 Title GEM: fast association study for the interplay of Gene, Environment and Methylation Version 1.5.0 Date 2015-12-05 Author Hong Pan, Joanna D Holbrook, Neerja
More informationAffymetrix GeneChip DNA Analysis Software
Affymetrix GeneChip DNA Analysis Software User s Guide Version 3.0 For Research Use Only. Not for use in diagnostic procedures. P/N 701454 Rev. 3 Trademarks Affymetrix, GeneChip, EASI,,,, HuSNP, GenFlex,
More informationsurvsnp: Power and Sample Size Calculations for SNP Association Studies with Censored Time to Event Outcomes
survsnp: Power and Sample Size Calculations for SNP Association Studies with Censored Time to Event Outcomes Kouros Owzar Zhiguo Li Nancy Cox Sin-Ho Jung Chanhee Yi June 29, 2016 1 Introduction This vignette
More informationHow to use cghmcr. October 30, 2017
How to use cghmcr Jianhua Zhang Bin Feng October 30, 2017 1 Overview Copy number data (arraycgh or SNP) can be used to identify genomic regions (Regions Of Interest or ROI) showing gains or losses that
More informationQuality control of array genotyping data with argyle Andrew P Morgan
Quality control of array genotyping data with argyle Andrew P Morgan 2015-10-08 Introduction Proper quality control of array genotypes is an important prerequisite to further analysis. Genotype quality
More informationPackage SimGbyE. July 20, 2009
Package SimGbyE July 20, 2009 Type Package Title Simulated case/control or survival data sets with genetic and environmental interactions. Author Melanie Wilson Maintainer Melanie
More informationRelease Notes. JMP Genomics. Version 4.0
JMP Genomics Version 4.0 Release Notes Creativity involves breaking out of established patterns in order to look at things in a different way. Edward de Bono JMP. A Business Unit of SAS SAS Campus Drive
More informationPolymorphism and Variant Analysis Lab
Polymorphism and Variant Analysis Lab Arian Avalos PowerPoint by Casey Hanson Polymorphism and Variant Analysis Matt Hudson 2018 1 Exercise In this exercise, we will do the following:. 1. Gain familiarity
More informationPRSice: Polygenic Risk Score software - Vignette
PRSice: Polygenic Risk Score software - Vignette Jack Euesden, Paul O Reilly March 22, 2016 1 The Polygenic Risk Score process PRSice ( precise ) implements a pipeline that has become standard in Polygenic
More informationGenetic Analysis. Page 1
Genetic Analysis Page 1 Genetic Analysis Objectives: 1) Set up Case-Control Association analysis and the Basic Genetics Workflow 2) Use JMP tools to interact with and explore results 3) Learn advanced
More informationImporting and Merging Data Tutorial
Importing and Merging Data Tutorial Release 1.0 Golden Helix, Inc. February 17, 2012 Contents 1. Overview 2 2. Import Pedigree Data 4 3. Import Phenotypic Data 6 4. Import Genetic Data 8 5. Import and
More informationAffymetrix Genotyping Console 3.0 User Manual
Affymetrix Genotyping Console 3.0 User Manual For research use only. Not for use in diagnostic procedures. Trademarks Affymetrix,, GeneChip, HuSNP, GenFlex, Flying Objective, CustomExpress, CustomSeq,
More informationExamples of implementation of pre-processing method described in paper with R code snippets - Electronic Supplementary Information (ESI)
Electronic Supplementary Material (ESI) for Analyst. This journal is The Royal Society of Chemistry 2015 Examples of implementation of pre-processing method described in paper with R code snippets - Electronic
More informationcrlmm to downstream data analysis
crlmm to downstream data analysis VJ Carey, B Carvalho March, 2012 1 Running CRLMM on a nontrivial set of CEL files To use the crlmm algorithm, the user must load the crlmm package, as described below:
More information500K Data Analysis Workflow using BRLMM
500K Data Analysis Workflow using BRLMM I. INTRODUCTION TO BRLMM ANALYSIS TOOL... 2 II. INSTALLATION AND SET-UP... 2 III. HARDWARE REQUIREMENTS... 3 IV. BRLMM ANALYSIS TOOL WORKFLOW... 3 V. RESULTS/OUTPUT
More informationTutorial. Identification of Variants Using GATK. Sample to Insight. November 21, 2017
Identification of Variants Using GATK November 21, 2017 Sample to Insight QIAGEN Aarhus Silkeborgvej 2 Prismet 8000 Aarhus C Denmark Telephone: +45 70 22 32 44 www.qiagenbioinformatics.com AdvancedGenomicsSupport@qiagen.com
More informationRecalling Genotypes with BEAGLECALL Tutorial
Recalling Genotypes with BEAGLECALL Tutorial Release 8.1.4 Golden Helix, Inc. June 24, 2014 Contents 1. Format and Confirm Data Quality 2 A. Exclude Non-Autosomal Markers......................................
More informationGBS Bioinformatics Pipeline(s) Overview
GBS Bioinformatics Pipeline(s) Overview Getting from sequence files to genotypes. Pipeline Coding: Ed Buckler Jeff Glaubitz James Harriman Presentation: Terry Casstevens With supporting information from
More informationPackage conumee. February 17, 2018
Package conumee February 17, 2018 Title Enhanced copy-number variation analysis using Illumina DNA methylation arrays Version 1.13.0 Author Volker Hovestadt, Marc Zapatka Maintainer Address Division of
More informationAn Introduction to the genoset Package
An Introduction to the genoset Package Peter M. Haverty April 4, 2013 Contents 1 Introduction 2 1.1 Creating Objects........................................... 2 1.2 Accessing Genome Information...................................
More informationGxE.scan. October 30, 2018
GxE.scan October 30, 2018 Overview GxE.scan can process a GWAS scan using the snp.logistic, additive.test, snp.score or snp.matched functions, whereas snp.scan.logistic only calls snp.logistic. GxE.scan
More informationPrepare input data for CINdex
1 Introduction Prepare input data for CINdex Genomic instability is known to be a fundamental trait in the development of tumors; and most human tumors exhibit this instability in structural and numerical
More informationKaryoStudio v1.4 User Guide
KaryoStudio v1.4 User Guide FOR RESEARCH USE ONLY ILLUMINA PROPRIETARY Part # 11328837 Rev. C June 2011 Notice This document and its contents are proprietary to Illumina, Inc. and its affiliates ("Illumina"),
More informationBEAGLECALL 1.0. Brian L. Browning Department of Medicine Division of Medical Genetics University of Washington. 15 November 2010
BEAGLECALL 1.0 Brian L. Browning Department of Medicine Division of Medical Genetics University of Washington 15 November 2010 BEAGLECALL 1.0 P a g e i Contents 1 Introduction... 1 1.1 Citing BEAGLECALL...
More informationUsing crlmm for copy number estimation and genotype calling with Illumina platforms
Using crlmm for copy number estimation and genotype calling with Illumina platforms Rob Scharpf November, Abstract This vignette illustrates the steps necessary for obtaining marker-level estimates of
More informationPackage saascnv. May 18, 2016
Version 0.3.4 Date 2016-05-10 Package saascnv May 18, 2016 Title Somatic Copy Number Alteration Analysis Using Sequencing and SNP Array Data Author Zhongyang Zhang [aut, cre], Ke Hao [aut], Nancy R. Zhang
More informationBICF Nano Course: GWAS GWAS Workflow Development using PLINK. Julia Kozlitina April 28, 2017
BICF Nano Course: GWAS GWAS Workflow Development using PLINK Julia Kozlitina Julia.Kozlitina@UTSouthwestern.edu April 28, 2017 Getting started Open the Terminal (Search -> Applications -> Terminal), and
More informationBuilding an R package
Division of Biostatistics School of Public Health, University of Minnesota Steps Prepare your functions, example data sets Build package structure using package.skeleton() Edit DESCRIPTION file Edit NAMESPACE
More informationMACAU User Manual. Xiang Zhou. March 15, 2017
MACAU User Manual Xiang Zhou March 15, 2017 Contents 1 Introduction 2 1.1 What is MACAU...................................... 2 1.2 How to Cite MACAU................................... 2 1.3 The Model.........................................
More informationMaximizing Public Data Sources for Sequencing and GWAS
Maximizing Public Data Sources for Sequencing and GWAS February 4, 2014 G Bryce Christensen Director of Services Questions during the presentation Use the Questions pane in your GoToWebinar window Agenda
More informationPROPER: PROspective Power Evaluation for RNAseq
PROPER: PROspective Power Evaluation for RNAseq Hao Wu [1em]Department of Biostatistics and Bioinformatics Emory University Atlanta, GA 303022 [1em] hao.wu@emory.edu October 30, 2017 Contents 1 Introduction..............................
More informationR.devices. Henrik Bengtsson. November 19, 2012
R.devices Henrik Bengtsson November 19, 2012 Abstract The R.devices package provides utility methods that enhance the existing graphical device functions already available in R for the purpose of simplifying
More informationPackage snpstatswriter
Type Package Package snpstatswriter February 20, 2015 Title Flexible writing of snpstats objects to flat files Version 1.5-6 Date 2013-12-05 Author Maintainer Write snpstats
More informationPackage LSDinterface
Type Package Title Reading LSD Results (.res) Files Version 0.3.1 Date 2017-11-24 Author Package LSDinterface November 27, 2017 Maintainer Interfaces R with LSD. Reads object-oriented
More informationPackage PTE. October 10, 2017
Type Package Title Personalized Treatment Evaluator Version 1.6 Date 2017-10-9 Package PTE October 10, 2017 Author Adam Kapelner, Alina Levine & Justin Bleich Maintainer Adam Kapelner
More informationReproducible Homerange Analysis
Reproducible Homerange Analysis (Sat Aug 09 15:28:43 2014) based on the rhr package This is an automatically generated file with all parameters and settings, in order to enable later replication of the
More informationAnalysis of two-way cell-based assays
Analysis of two-way cell-based assays Lígia Brás, Michael Boutros and Wolfgang Huber April 16, 2015 Contents 1 Introduction 1 2 Assembling the data 2 2.1 Reading the raw intensity files..................
More informationImage Analysis with beadarray
Mike Smith October 30, 2017 Introduction From version 2.0 beadarray provides more flexibility in the processing of array images and the extraction of bead intensities than its predecessor. In the past
More informationSIBER User Manual. Pan Tong and Kevin R Coombes. May 27, Introduction 1
SIBER User Manual Pan Tong and Kevin R Coombes May 27, 2015 Contents 1 Introduction 1 2 Using SIBER 1 2.1 A Quick Example........................................... 1 2.2 Dealing With RNAseq Normalization................................
More informationPeter Schweitzer, Director, DNA Sequencing and Genotyping Lab
The instruments, the runs, the QC metrics, and the output Peter Schweitzer, Director, DNA Sequencing and Genotyping Lab Overview Roche/454 GS-FLX 454 (GSRunbrowser information) Evaluating run results Errors
More informationAgilent Genomic Workbench 7.0
Agilent Genomic Workbench 7.0 Workflow User Guide For Research Use Only. Not for use in diagnostic procedures. Agilent Technologies Notices Agilent Technologies, Inc. 2012, 2015 No part of this manual
More informationRnBeadsDJ A Quickstart Guide to the RnBeads Data Juggler
RnBeadsDJ A Quickstart Guide to the RnBeads Data Juggler Fabian Müller, Yassen Assenov, Pavlo Lutsik Contact: rnbeads@mpi-inf.mpg.de Package version: 1.12.2 September 25, 2018 RnBeads is an R package for
More informationRLMM - Robust Linear Model with Mahalanobis Distance Classifier
RLMM - Robust Linear Model with Mahalanobis Distance Classifier Nusrat Rabbee and Gary Wong June 13, 2018 Contents 1 Introduction 1 2 Instructions for Genotyping Affymetrix Mapping 100K array - Xba set
More informationIntroduction to GE Microarray data analysis Practical Course MolBio 2012
Introduction to GE Microarray data analysis Practical Course MolBio 2012 Claudia Pommerenke Nov-2012 Transkriptomanalyselabor TAL Microarray and Deep Sequencing Core Facility Göttingen University Medical
More informationBuilding Map File and Mapping Object
Building Map File and Mapping Object Lori Shepherd and Jonathan Dare July 10, 2007 Contents Please download buildmapexamples. The acghplus package requires that a mapping object be built off a mapping
More informationAn Introduction to the methylumi package
An Introduction to the methylumi package Sean Davis and Sven Bilke October 13, 2014 1 Introduction Gene expression patterns are very important in understanding any biologic system. The regulation of gene
More informationGenetic Programming. Charles Chilaka. Department of Computational Science Memorial University of Newfoundland
Genetic Programming Charles Chilaka Department of Computational Science Memorial University of Newfoundland Class Project for Bio 4241 March 27, 2014 Charles Chilaka (MUN) Genetic algorithms and programming
More informationMPG NGS workshop I: Quality assessment of SNP calls
MPG NGS workshop I: Quality assessment of SNP calls Kiran V Garimella (kiran@broadinstitute.org) Genome Sequencing and Analysis Medical and Population Genetics February 4, 2010 SNP calling workflow Filesize*
More informationPackage ridge. R topics documented: February 15, Title Ridge Regression with automatic selection of the penalty parameter. Version 2.
Package ridge February 15, 2013 Title Ridge Regression with automatic selection of the penalty parameter Version 2.1-2 Date 2012-25-09 Author Erika Cule Linear and logistic ridge regression for small data
More informationSTAT 540 Computing in Statistics
STAT 540 Computing in Statistics Introduces programming skills in two important statistical computer languages/packages. 30-40% R and 60-70% SAS Examples of Programming Skills: 1. Importing Data from External
More informationsegmentseq: methods for detecting methylation loci and differential methylation
segmentseq: methods for detecting methylation loci and differential methylation Thomas J. Hardcastle October 30, 2018 1 Introduction This vignette introduces analysis methods for data from high-throughput
More informationBioinformatics - Homework 1 Q&A style
Bioinformatics - Homework 1 Q&A style Instructions: in this assignment you will test your understanding of basic GWAS concepts and GenABEL functions. The materials needed for the homework (two datasets
More informationConvert Dosages to Genotypes Author: Autumn Laughbaum, Golden Helix, Inc.
Convert Dosages to Genotypes Author: Autumn Laughbaum, Golden Helix, Inc. Overview This script converts allelic dosage values to genotypes based on user-specified thresholds. The dosage data may be in
More informationcbioportal https://www.ncbi.nlm.nih.gov/pubmed/ /5/401
cbioportal http://www.cbioportal.org/ https://www.ncbi.nlm.nih.gov/pubmed/23550210 http://cancerdiscovery.aacrjournals.org/content/ 2/5/401 Tutorials http://www.cbioportal.org/tutorial.jsp http://www.cbioportal.org/faq.jsp
More informationTriform: peak finding in ChIP-Seq enrichment profiles for transcription factors
Triform: peak finding in ChIP-Seq enrichment profiles for transcription factors Karl Kornacker * and Tony Håndstad October 30, 2018 A guide for using the Triform algorithm to predict transcription factor
More information1. Summary statistics test_gwas. This file contains a set of 50K random SNPs of the Subjective Well-being GWAS of the Netherlands Twin Register
Quality Control for Genome-Wide Association Studies Bart Baselmans & Meike Bartels Boulder 2017 Setting up files and directories To perform a quality control protocol in a Genome-Wide Association Meta
More informationThe Lander-Green Algorithm in Practice. Biostatistics 666
The Lander-Green Algorithm in Practice Biostatistics 666 Last Lecture: Lander-Green Algorithm More general definition for I, the "IBD vector" Probability of genotypes given IBD vector Transition probabilities
More informationProcessing of the Chang data
Processing of the Chang data Kevin R. Coombes, Jing Wang, and Keith A. Baggerly 16 March 2007 1 Description of the problem We need to load the Chang data, but we have two choices. We can start with the
More informationGenetic type 1 Error Calculator (GEC)
Genetic type 1 Error Calculator (GEC) (Version 0.2) User Manual Miao-Xin Li Department of Psychiatry and State Key Laboratory for Cognitive and Brain Sciences; the Centre for Reproduction, Development
More informationsegmentseq: methods for detecting methylation loci and differential methylation
segmentseq: methods for detecting methylation loci and differential methylation Thomas J. Hardcastle October 13, 2015 1 Introduction This vignette introduces analysis methods for data from high-throughput
More informationIntroduction to Statistics using R/Rstudio
Introduction to Statistics using R/Rstudio R and Rstudio Getting Started Assume that R for Windows and Macs already installed on your laptop. (Instructions for installations sent) R on Windows R on MACs
More informationPcAux. Kyle M. Lang, Jacob Curtis, Daniel E Bontempo Institute for Measurement, Methodology, Analysis & Policy at Texas Tech University May 5, 2017
PcAux Kyle M. Lang, Jacob Curtis, Daniel E Bontempo Institute for Measurement, Methodology, Analysis & Policy at Texas Tech University May 5, 2017 What is PcAux? PcAux is an R package that uses the Principal
More informationSNP HiTLink Manual. Yoko Fukuda 1, Hiroki Adachi 2, Eiji Nakamura 2, and Shoji Tsuji 1
SNP HiTLink Manual Yoko Fukuda 1, Hiroki Adachi 2, Eiji Nakamura 2, and Shoji Tsuji 1 1 Department of Neurology, Graduate School of Medicine, the University of Tokyo, Tokyo, Japan 2 Dynacom Co., Ltd, Kanagawa,
More informationRelease Note. Agilent Genomic Workbench 6.5 Lite
Release Note Agilent Genomic Workbench 6.5 Lite Associated Products and Part Number # G3794AA G3799AA - DNA Analytics Software Modules New for the Agilent Genomic Workbench SNP genotype and Copy Number
More informationPackage EMLRT. August 7, 2014
Package EMLRT August 7, 2014 Type Package Title Association Studies with Imputed SNPs Using Expectation-Maximization-Likelihood-Ratio Test LazyData yes Version 1.0 Date 2014-08-01 Author Maintainer
More informationPackage inversion. R topics documented: July 18, Type Package. Title Inversions in genotype data. Version
Package inversion July 18, 2013 Type Package Title Inversions in genotype data Version 1.8.0 Date 2011-05-12 Author Alejandro Caceres Maintainer Package to find genetic inversions in genotype (SNP array)
More informationData input & output. Hadley Wickham. Assistant Professor / Dobelman Family Junior Chair Department of Statistics / Rice University.
Data input & output Hadley Wickham Assistant Professor / Dobelman Family Junior Chair Department of Statistics / Rice University June 2012 1. Working directories 2. Loading data 3. Strings and factors
More informationCalling variants in diploid or multiploid genomes
Calling variants in diploid or multiploid genomes Diploid genomes The initial steps in calling variants for diploid or multi-ploid organisms with NGS data are the same as what we've already seen: 1. 2.
More informationWelcome to MAPHiTS (Mapping Analysis Pipeline for High-Throughput Sequences) tutorial page.
Welcome to MAPHiTS (Mapping Analysis Pipeline for High-Throughput Sequences) tutorial page. In this page you will learn to use the tools of the MAPHiTS suite. A little advice before starting : rename your
More informationGemTools Documentation
Literature: GemTools Documentation Bert Klei and Brian P. Kent February 2011 This software is described in GemTools: a fast and efficient approach to estimating genetic ancestry (in preparation) Klei L,
More informationPackage HMMASE. February 4, HMMASE R package
Package HMMASE February 4, 2014 Type Package Title HMMASE R package Version 1.0 Date 2014-02-04 Author Juan R. Steibel, Heng Wang, Ping-Shou Zhong Maintainer Heng Wang An R package that
More informationStep-by-Step Guide to Relatedness and Association Mapping Contents
Step-by-Step Guide to Relatedness and Association Mapping Contents OBJECTIVES... 2 INTRODUCTION... 2 RELATEDNESS MEASURES... 2 POPULATION STRUCTURE... 6 Q-K ASSOCIATION ANALYSIS... 10 K MATRIX COMPRESSION...
More informationGSEQ Software User s Guide for AccuID TM
GSEQ Software User s Guide for AccuID TM This protocol is majorly based on Affymetrix GeneChip Sequence Analysis Software User s Guide Version 4.1 guidebook and has modified it to fit GSEQ Software User
More informationVariant calling using SAMtools
Variant calling using SAMtools Calling variants - a trivial use of an Interactive Session We are going to conduct the variant calling exercises in an interactive idev session just so you can get a feel
More informationRelease Note. Agilent Genomic Workbench Standard Edition
Release te Agilent Genomic Workbench Standard Edition 5.0.14 A. New for the Agilent Genomic Workbench DNA Analytics becomes a component of Agilent Genomic Workbench DNA Analytics becomes a component of
More informationPackage globalgsa. February 19, 2015
Type Package Package globalgsa February 19, 2015 Title Global -Set Analysis for Association Studies. Version 1.0 Date 2013-10-22 Author Natalia Vilor, M.Luz Calle Maintainer Natalia Vilor
More informationThe software comes with 2 installers: (1) SureCall installer (2) GenAligners (contains BWA, BWA-MEM).
Release Notes Agilent SureCall 3.5 Product Number G4980AA SureCall Client 6-month named license supports installation of one client and server (to host the SureCall database) on one machine. For additional
More informationSOLOMON: Parentage Analysis 1. Corresponding author: Mark Christie
SOLOMON: Parentage Analysis 1 Corresponding author: Mark Christie christim@science.oregonstate.edu SOLOMON: Parentage Analysis 2 Table of Contents: Installing SOLOMON on Windows/Linux Pg. 3 Installing
More informationDjork-Arné Clevert and Andreas Mitterecker. Institute of Bioinformatics, Johannes Kepler University Linz
Software Manual Institute of Bioinformatics, Johannes Kepler University Linz cn.farms: a latent variable model to detect copy number variations in microarray data with a low false discovery rate Manual
More informationThe analysis of acgh data: Overview
The analysis of acgh data: Overview JC Marioni, ML Smith, NP Thorne January 13, 2006 Overview i snapcgh (Segmentation, Normalisation and Processing of arraycgh data) is a package for the analysis of array
More informationToCatchAThief c ryan campbell & jenn coughlan 7/23/2018
ToCatchAThief c ryan campbell & jenn coughlan 7/23/2018 Welcome to the To Catch a Thief: With Data! walkthrough! https://bioconductor.org/packages/devel/ bioc/vignettes/snprelate/inst/doc/snprelatetutorial.html
More informationPackage semisup. March 10, Version Title Semi-Supervised Mixture Model
Version 1.7.1 Title Semi-Supervised Mixture Model Package semisup March 10, 2019 Description Useful for detecting SNPs with interactive effects on a quantitative trait. This R packages moves away from
More informationStep-by-Step Guide to Basic Genetic Analysis
Step-by-Step Guide to Basic Genetic Analysis Page 1 Introduction This document shows you how to clean up your genetic data, assess its statistical properties and perform simple analyses such as case-control
More informationImport GEO Experiment into Partek Genomics Suite
Import GEO Experiment into Partek Genomics Suite This tutorial will illustrate how to: Import a gene expression experiment from GEO SOFT files Specify annotations Import RAW data from GEO for gene expression
More informationThe software comes with 2 installers: (1) SureCall installer (2) GenAligners (contains BWA, BWA- MEM).
Release Notes Agilent SureCall 4.0 Product Number G4980AA SureCall Client 6-month named license supports installation of one client and server (to host the SureCall database) on one machine. For additional
More informationQTL DETECTION WITH GS3. A Legarra, 19/06/2013
QTL DETECTION WITH GS3. A Legarra, 19/06/2013 Hélène Gilbert simulated a data set with 2 chromosomes and 5 QTLs at positions 20 cl, 73.14 cm on chromosome 1 and 20, 33.24 and 59.60 in chromosome 2. Each
More informationPackage cghmcr. March 6, 2019
Version 1.40.0 Package cghmcr March 6, 2019 Title Find chromosome regions showing common gains/losses Author J. Zhang and B. Feng Maintainer J. Zhang Depends methods, DNAcopy,
More information