calc_MW: Calculate molecular weights

In sybilccFBA: Cost Constrained Flux Balance Analysis (ccFBA): MetabOlic Modeling with ENzyme kineTics (MOMENT)

Description

Calculate Molecular weights of different proteins using the genome .faa file.

Usage

calc_MW(aa_fname = "aa.txt", ptt_fname = "test2.ptt", faa_fname = "NC_000913.faa",
 nchrm = 1)

Arguments

aa_fname	file name of file containing lisT of amino acid names
ptt_fname	file name of file containing gene names with gene code
faa_fname	file name of file containing gene code and sequence of amino acids
nchrm	the number of chromosomes in the genome

Value

generate a file containing gene name , length, and molecular weight

Author(s)

Abdelmoneim Amer Desouki

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.
## Not run: 
aa_fname <- system.file("extdata", "aa.txt", package="sybilccFBA")
ptt_fname <- system.file("extdata", "test2.ptt", package="sybilccFBA")
faa_fname <- system.file("extdata", "NC_000913.faa", package="sybilccFBA")

geneCnt <- calc_MW(aa_fname,ptt_fname,faa_fname)
write.csv(file="geneCnt.csv",geneCnt) 

## The function is currently defined as
"calc_MW"

## End(Not run)

sybilccFBA documentation built on Dec. 16, 2019, 1:34 a.m.