Calculate molecular weights

Description

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

Usage

1
2
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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
##---- 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)