[BioC] Protein/peptide mass

[john seers (IFR)]

Hi Sean

That works fine for me from home so I just need to solve why I cannot
get round the proxy server at work.

Thanks very much for your help.

This thread has been a help in solving a couple of side issues as well.
I had not heard of Rcurl either so I will have a look at that.


Regards


John Seers





 
---




 x <- url(' 
http://ca.expasy.org/cgi-bin/pi_tool?protein=MKWVTFISLLFLFSSAYS&resoluti
on=m
onoisotopic')
 res <- readLines(x)
 

You'll notice that one of the lines of res is something like:

[68] "Theoretical pI/Mw: 8.34 / 2139.11"

You can then use typical R tools like gsub/grep to find what you need
like
so:

 as.numeric(gsub('.*/ ','',res[grep('Theoretical pI/Mw:',res)]))

Which will return:

 2139.11