Entering edit mode
Hi Herve (and others),
in an attempt to update annotations for some probes, and after taking
your
workshop at Bioc2011, I have been playing with Biostrings, BSgenomes,
GRanges, et al, mostly to good effect (thank you and others for
writing and
documenting them, they're terrific once the learning curve flattens
out).
However, I rebuilt the SNPlocs package following instructions you
long ago
posted for Paul Shannon, and I think I have hit a wall in that
respect. I
put the offending package (and some others) up:
http://flaver.com/annotations/ (in case you need the actual package
for an
autopsy)
When I load and inject SNPlocs.Hsapiens.dbSNP.20101109, your package
of
build 132, everything goes as planned: I can inject the hard mask with
polymorphisms where they're supposed to be, and everything is hunkey
dorey.
When I do this after building (without apparent incident) the b134
based
SNPlocs.Hsapiens.dbSNP.20110815_0.99.3 , after editing the scripts to
look
for 'GRCh37.p2' instead of 'GRCh37' and rebuilding, I get the
following
error:
Error in x %in% SNPlocs_seqnames :
error in evaluating the argument 'table' in selecting a method for
function '%in%': Error in getSNPcount() : internal error: 'SNPcount'
data
set is broken.
Please contact the maintainer of the
SNPlocs.Hsapiens.dbSNP.20110815
package.
This seems odd, since the internal regression tests all passed when I
built
the package. But, obviously I must have missed something somewhere!
Any thoughts?
Thanks!
--
If people do not believe that mathematics is simple, it is only
because they
do not realize how complicated life is.
John von Neumann<http: www-groups.dcs.st-="" and.ac.uk="" ~history="" biographies="" von_neumann.html="">
[[alternative HTML version deleted]]