Last data update: 2014.03.03

R: Calculate molecular weight
calcMWR Documentation

Calculate molecular weight

Description

Function to calculate molecular weight of a fragment generated by the MassCLEAVE assay for either the T or C cleavage reactions

Usage

calcMW(sequence, extra = c("5OH", "5PPP-3P", "5PPP-3OH"), adduct = c("", "Na", "K"), rxn = c("T", "C"))

Arguments

sequence

Nucleotide sequence input

extra

One of "5OH" (default), "5PPP-3P", or "5PPP-3OH" describing 5' and/or 3' modifications of the fragment

adduct

One of 'Na', 'K', or ” (default) specifying whether or not the molecular weight should be calculated for a salt adduct

rxn

One of 'T' or 'C' indicating which cleavage reaction is employed to generate the fragment

Value

Returns a numerical output corresponding to the molecular weight (in Da) of sequence input. Note that the output may actually represent multiple molecular weights if/whenever the input sequence contains one or more degenerate bases (e.g. R or Y).

Author(s)

Reid F. Thompson (rthompso@aecom.yu.edu), John M. Greally (jgreally@aecom.yu.edu)

Examples

calcMW("AAATCCC")
calcMW("AARTYCC")

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(MassArray)
> png(filename="/home/ddbj/snapshot/RGM3/R_BC/result/MassArray/calcMW.Rd_%03d_medium.png", width=480, height=480)
> ### Name: calcMW
> ### Title: Calculate molecular weight
> ### Aliases: calcMW
> ### Keywords: character arith
> 
> ### ** Examples
> 
> calcMW("AAATCCC")
[1] 2180.377
> calcMW("AARTYCC")
[1] 2197.361 2180.377 2213.360 2196.376
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>