Hi there ...
I was looking to work with some of the Bioconductor annotation packages
directly in Python and came across Rpy2. The first problem is just one of
inexperience, I'm probably doing things all wrong and that could be the
source of my issue :) I did just notice something odd:
t = robjects.r('get("1000_at", hgu95av2SYMBOL)')
>>> t
<RVector - Python:0x2aaab0e9ca70 / R:0x233c208>
>>> t[0]
'MAPK3'
So far so good.
t = robjects.r('get("1017_at", hgu95av2SYMBOL)')
## note that in R, this returns NA
>>> t
<RVector - Python:0x2aaab0e9cb90 / R:0x207dfc8>
>>> t[0]
True
Is this expected behavior (True replacing NA)?
Also, while I'm here, is there a way to skip the intermediate step, ie
have the singular Python object returned instead of an RVector?
Thanks
-J
------------------------------------------------------------------------------
This SF.net email is sponsored by:
SourcForge Community
SourceForge wants to tell your story.
http://p.sf.net/sfu/sf-spreadtheword
_______________________________________________
rpy-list mailing list
rpy-list@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/rpy-list
