Package HGNChelper. August 29, 2017

Similar documents
This report is based on sampled data. Jun 1 Jul 6 Aug 10 Sep 14 Oct 19 Nov 23 Dec 28 Feb 1 Mar 8 Apr 12 May 17 Ju

Asks for clarification of whether a GOP must communicate to a TOP that a generator is in manual mode (no AVR) during start up or shut down.

Monthly SEO Report. Example Client 16 November 2012 Scott Lawson. Date. Prepared by

software.sci.utah.edu (Select Visitors)

Section 1.2: What is a Function? y = 4x

Polycom Advantage Service Endpoint Utilization Report

Polycom Advantage Service Endpoint Utilization Report

ICT PROFESSIONAL MICROSOFT OFFICE SCHEDULE MIDRAND

New Concept for Article 36 Networking and Management of the List

SCI - software.sci.utah.edu (Select Visitors)

HPE Security Data Security. HPE SecureData. Product Lifecycle Status. End of Support Dates. Date: April 20, 2017 Version:

Seattle (NWMLS Areas: 140, 380, 385, 390, 700, 701, 705, 710) Summary

Technical Manual Index

COURSE LISTING. Courses Listed. Training for Database & Technology with Modeling in SAP HANA. 20 November 2017 (12:10 GMT) Beginner.

Undergraduate Admission File

Annex A to the DVD-R Disc and DVD-RW Disc Patent License Agreement Essential Sony Patents relevant to DVD-RW Disc

BANGLADESH UNIVERSITY OF PROFESSIONALS ACADEMIC CALENDAR FOR MPhil AND PHD PROGRAM 2014 (4 TH BATCH) PART I (COURSE WORK)

ACTIVE MICROSOFT CERTIFICATIONS:

All King County Summary Report

CBERS-2. Attitude Control and its Effects on Image Geometric Correction. Follow up to TCM-06 INPE CBERS TEAM

Previous Intranet Initial intranet created in 2002 Created solely by Information Systems Very utilitarian i Created to permit people to access forms r

AIMMS Function Reference - Date Time Related Identifiers

MISO PJM Joint and Common Market Cross Border Transmission Planning

Programming for Engineers Structures, Unions

UAE PUBLIC TRAINING CALENDAR

DATE OF BIRTH SORTING (DBSORT)

NCC Cable System Order

Seattle (NWMLS Areas: 140, 380, 385, 390, 700, 701, 705, 710) Summary

MSRS Roadmap. As of January 15, PJM 2019

Seattle (NWMLS Areas: 140, 380, 385, 390, 700, 701, 705, 710) Summary

Excel Functions & Tables

COURSE LISTING. Courses Listed. with SAP Hybris Marketing Cloud. 24 January 2018 (23:53 GMT) HY760 - SAP Hybris Marketing Cloud

S-101. The New ENC Product Specification. Julia Powell S-100 Working Group Chair

3. EXCEL FORMULAS & TABLES

Grade 4 Mathematics Pacing Guide

Annex A to the MPEG Audio Patent License Agreement Essential Philips, France Telecom and IRT Patents relevant to DVD-Video Disc - MPEG Audio - general

CIMA Asia. Interactive Timetable Live Online

DAS LRS Monthly Service Report

Technical Manual Index

MHA LPI Date Correction Request Process

Package RcppBDT. August 29, 2016

2018 CALENDAR OF ACTIVITIES

More Binary Search Trees AVL Trees. CS300 Data Structures (Fall 2013)

Budget Transfers. To initiate a budget transfer, go to FGAJVCM (Journal Voucher Mass Entry). FGAJVCM

COURSE LISTING. Courses Listed. with Governance, Risk and Compliance (GRC) SAP BusinessObjects. 19 February 2018 (15:13 GMT) GRC100 -

More BSTs & AVL Trees bstdelete

Package taskscheduler

Embedded Indexing. Lucie Haskins.

DoD Environmental Security Technology Certification Program (ESTCP) Tim Tetreault DoD August 15, 2017

Freedom of Information Act 2000 reference number RFI

NMOSE GPCD CALCULATOR

Town of Georgetown Comprehensive Plan

CIMA Certificate BA Interactive Timetable

Next Steps for WHOIS Accuracy Global Domains Division. ICANN June 2015

SF Current Cumulative PTF Package. I B M i P R E V E N T I V E S E R V I C E P L A N N I N G I N F O R M A T I O N

Peak Season Metrics Summary

FREQUENTLY ASKED QUESTIONS

Peak Season Metrics Summary

2

Peak Season Metrics Summary

Martin Purvis CEO. Presentation following the AGM. 26 April 2012

Database Programming with SQL

Yield Reduction Due to Shading:

2

2

San Francisco Housing Authority (SFHA) Leased Housing Programs October 2015

ADVANCED SPREADSHEET APPLICATIONS (07)

OBJECT-ORIENTED PROGRAMMING IN R: S3 & R6. Environments, Reference Behavior, & Shared Fields

HPE Aruba. Course Training Year 2017 By IT Green

WHOIS Accuracy Reporting System (ARS): Phase 2 Cycle 1 Results Webinar 12 January ICANN GDD Operations NORC at the University of Chicago

CIMA Asia. Interactive Timetable Live Online

RDAP Implementation. Francisco Arias & Gustavo Lozano 21 October 2015

SF Current Cumulative PTF Package. I B M i P R E V E N T I V E S E R V I C E P L A N N I N G I N F O R M A T I O N

The State of the Raven. Jon Warbrick University of Cambridge Computing Service

Certificate in Security Management

ACTIAN PRODUCTS by Platform - Vector, Vector in Hadoop as of October 18, 2017

Withdrawal of Text Processing (Business Professional)

Fluidity Trader Historical Data for Ensign Software Playback

SF Current Cumulative PTF Package. I B M i P R E V E N T I V E S E R V I C E P L A N N I N G I N F O R M A T I O N

Obtaining and Managing IP Addresses. Xavier Le Bris IP Resource Analyst - Trainer

Data Types. 9. Types. a collection of values and the definition of one or more operations that can be performed on those values

Peak Season Metrics Summary

LocatorHub Product Life Cycle Status

SD3-60 AIRCRAFT MAINTENANCE MANUAL REF. No 360/MM

NIR February JPNIC Updates. Hiroki Kawabata Japan Network Information Center (JPNIC) Copyright 2016 Japan Network Information Center

Imputation for missing observation through Artificial Intelligence. A Heuristic & Machine Learning approach

ACTIVE MICROSOFT CERTIFICATIONS:

Automatic Renewal Using DIY Technology to Create an Improved Patron Experience

Create Designs. How do you draw a design on a coordinate grid?

Turn Stacked Summarized Bars into One Bridged Bar. Document No: SB001-P3 Revision No: Last Updated: March 18, 2003

Mpoli Archive - FTP Site Statistics. Top 20 Directories Sorted by Disk Space

CMR India Monthly Mobile Handset Market Report. June 2017

Package assertive.data

INFORMATION TECHNOLOGY SPREADSHEETS. Part 1

COURSE LISTING. Courses Listed. Training for Cloud with SAP Cloud Platform in Development. 23 November 2017 (08:12 GMT) Beginner.

Dashboard. Jan 13, Jan 8, 2012 Comparing to: Site. 12,742 Visits % Bounce Rate. 00:05:26 Avg. Time on Site.

COURSE LISTING. Courses Listed. Training for Database & Technology with Development in SAP Cloud Platform. 1 December 2017 (22:41 GMT) Beginner

GWDG Software Archive - FTP Site Statistics. Top 20 Directories Sorted by Disk Space

Spiegel Research 3.0 The Mobile App Story

Tracking the Internet s BGP Table

Transcription:

Package HGNChelper August 29, 2017 Maintainer Levi Waldron <lwaldron.research@gmail.com> Depends R (>= 2.10), methods, utils Author Version 0.3.6 Date 2017-08-29 License GPL (>=2.0) Title Handy Functions for Working with HGNC Gene Symbols and Affymetri Probeset Identifiers Contains functions for identifying and correcting HGNC gene symbols which have been converted to date format by Ecel, for reversibly converting between HGNC symbols and valid R names, identifying invalid HGNC symbols and correcting synonyms and outdated symbols which can be mapped to an official symbol. URL https://bitbucket.org/lwaldron/hgnchelper BugReports https://bitbucket.org/lwaldron/hgnchelper/issues R topics documented: HGNChelper-package.................................... 2 affytor........................................... 2 checkgenesymbols..................................... 3 findecelgenesymbols................................... 4 hgnc.table.......................................... 5 rtoaffy........................................... 5 rtosymbol......................................... 6 symboltor......................................... 7 Inde 8 1

2 affytor HGNChelper-package Handy functions for working with HGNC gene symbols and Affymetri probeset identifiers. Details Contains functions for identifying and correcting HGNC gene symbols which have been converted to date format by Ecel, for reversibly converting between HGNC symbols and valid R names, identifying invalid HGNC symbols and correcting synonyms and outdated symbols which can be mapped to an official symbol. Package: HGNChelper Maintainer: Levi Waldron <lwaldron@hsph.harvard.edu> Depends: R (>= 2.10) Author: Version: 0.3.6 Date: 2017-08-29 License: GPL (>2.0) Title: Handy functions for working with HGNC gene symbols and Affymetri probeset identifiers. URL: https://bitbucket.org/lwaldron/hgnchelper BugReports: https://bitbucket.org/lwaldron/hgnchelper affytor function to convert Affymetri probeset identifiers to valid R names This function simply prepends "affy." to the probeset IDs to create valid R names. Reverse operation is done by the rtoaffy function. affytor() vector of Affymetri probeset identifiers, or any identifier which may with a digit.

checkgenesymbols 3 a character vector that is simply with "affy." prepended to each value. checkgenesymbols function to identify outdated or Ecel-mogrified gene symbols This function identifies gene symbols which are outdated or may have been mogrified by Ecel or other spreadsheet programs. If output is assigned to a variable, it returns a data.frame of the same number of rows as the input, with a second column indicating whether the symbols are valid and a third column with a corrected gene list. checkgenesymbols(, unmapped.as.na=true, hgnc.table=null) Vector of gene symbols to check for mogrified or outdated values unmapped.as.na If TRUE, unmapped symbols will appear as NA in the Suggested.Symbol column. If FALSE, the original unmapped symbol will be kept when no correction can be found. hgnc.table If hgnc.table is a data.frame with colnames(hgnc.table) identical to c("symbol", "Approved.Symbol"), it is used to correct gene symbols in. Otherwise, the default table data("hgnc.table", package="hgnchelper") is used. The function used for creating this the hgnc.table dataframe can be found in the inst/hgnclookup.r file, in the source of this package. By default the hgnc.table dataframe shipped with this package is used (see?hgnc.table) The function will return a data.frame of the same number of rows as the input, with corrections possible from hgnc.table. See Also hgnc.table

4 findecelgenesymbols Eamples library(hgnchelper) = c("fn1", "TP53", "UNKNOWNGENE","7-Sep", "9/7", "1-Mar", "Oct4", "4-Oct", "OCT4-PG4", "C19ORF71", "C19orf71") res <- checkgenesymbols() res if (interactive()){ ##Run checkgenesymbols with a brand-new map downloaded from HGNC: source(system.file("hgnclookup.r", package = "HGNChelper")) ##You should save this if you are going to use it multiple times, ##then load it from file rather than burdening HGNC's servers. ## save(hgnc.table, file="hgnc.table.rda", compress="bzip2") ## load("hgnc.table.rda") res <- checkgenesymbols(, hgnc.table=hgnc.table) } findecelgenesymbols function to identify Ecel-mogrified gene symbols This function identifies gene symbols which may have been mogrified by Ecel or other spreadsheet programs. If output is assigned to a variable, it returns a vector of the same length where symbols which could be mapped have been mapped. Note that this function is superceded by checkgenesymbols, which corrects Ecel-mogrified gene symbols as well as aliases and outdated symbols. findecelgenesymbols(, mog.map = read.csv(system.file("etdata/mog_map.csv", package = "HGNChelper"), as.is = TRUE), rege = "[0-9]\\-(JAN FEB MAR APR MAY JUN JUL AUG SEP OCT NOV DEC) [0-9]\\.[0-9] [0-9]E\\+[[0-9][0-9]") mog.map rege Vector of gene symbols to check for mogrified values Map of known mogrifications. A default map is available with this package by data(mog.map), but any map may be used. This should be a dataframe with two columns: original and mogrified, containing the correct and incorrect symbols, respectively. Regular epression, recognized by the base::grep function which is called with ignore.case=true, to identify mogrified symbols. It is not necessary for all matches to have a corresponding entry in mog.map$mogrified; values in which are matched by this rege but are not found in mog.map$mogrified simply will not be corrected. This rege is based that provided by Zeeberg et al., Mistaken Identifiers: Gene name errors can be introduced inadvertently when using Ecel in bioinformatics. BMC Bioinformatics 2004, 5:80.

hgnc.table 5 if the return value of the function is assigned to a variable, the function will return a vector of the same length as the input, with corrections possible from mog.map made. See Also checkgenesymbols hgnc.table All current and withdrawn HGNC symbols. All current and withdrawn symbols from genenames.org. Format A dataframe with the first column providing a gene symbol or known alias (including withdrawn symbols), second column providing the approved HGNC symbol. Source Etracted from genenames.org prior to package build. Eamples data("hgnc.table", package="hgnchelper", envir=environment()) head(hgnc.table) rtoaffy function to convert the output of affytor back to the original Affymetri probeset identifiers. This function simply strips the "affy." added by the affytor function. rtoaffy() the character vector returned by the affytor function.

6 rtosymbol a character vector of Affymetri probeset identifiers. rtosymbol function to reverse the conversion made by symboltor This function reverses the actions of the symboltor function. rtosymbol() the character vector returned by the symboltor function. a character vector of HGNC gene symbols, which are not in general valid R names. Eamples library(hgnchelper) data("hgnc.table", envir=environment()) hgnc.symbols <- as.character(na.omit(unique(hgnc.table[,2]))) if(!identical(all.equal(hgnc.symbols, rtosymbol(make.names(symboltor(hgnc.symbols)))), TRUE)) stop("hgnc mapping was not reversible.")

symboltor 7 symboltor function to *reversibly* convert HGNC gene symbols to valid R names. This function reversibly converts HGNC gene symbols to valid R names by prepending "symbol.", and making the following substitutions: "-" to "hyphen", "@" to "ampersand", and "/" to "forwardslash". symboltor() vector of HGNC symbols a vector of valid R names, of the same length as, which can be converted to the same HGNC symbols using the rtosymbol function. Eamples library(hgnchelper) data("hgnc.table", envir=environment()) hgnc.symbols <- as.character(na.omit(unique(hgnc.table[,2]))) if(!identical(all.equal(hgnc.symbols, rtosymbol(make.names(symboltor(hgnc.symbols)))), TRUE)) stop("hgnc mapping was not reversible.")

Inde Topic datasets hgnc.table, 5 Topic package HGNChelper-package, 2 affytor, 2, 5 checkgenesymbols, 3 findecelgenesymbols, 4 hgnc.table, 5 HGNChelper-package, 2 rtoaffy, 2, 5 rtosymbol, 6 symboltor, 7 8