[BioC] an error from ShortRead
wang peter
wng.peter at gmail.com
Fri Mar 16 17:59:29 CET 2012
it works well on
qual <- SolexaQuality(quality(quality(reads)))
but
that is the error place:
myqual_mat <- matrix(charToRaw(as.character(unlist(qual))),nrow=length(qual),
byrow=TRUE)
Error in .Call(.NAME, ..., PACKAGE = PACKAGE) :
negative length vectors are not allowed
Calls: trimRead ... matrix -> charToRaw -> unlist -> unlist -> .Call2 -> .Call
Execution halted
>> R version 2.14.1 (2011-12-22)
>> Platform: x86_64-redhat-linux-gnu (64-bit)
>>
>> locale:
>> [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
>> [3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8
>> [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
>> [7] LC_PAPER=C LC_NAME=C
>> [9] LC_ADDRESS=C LC_TELEPHONE=C
>> [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
>>
>> attached base packages:
>> [1] stats graphics grDevices utils datasets methods base
>>
>> --
>> shan gao
>> Room 231(Dr.Fei lab)
>> Boyce Thompson Institute
>> Cornell University
>> Tower Road, Ithaca, NY 14853-1801
>> Office phone: 1-607-254-1267(day)
>> Official email:sg839 at cornell.edu
>> Facebook:http://www.facebook.com/profile.php?id=100001986532253
>>
>> _______________________________________________
>> Bioconductor mailing list
>> Bioconductor at r-project.org
>> https://stat.ethz.ch/mailman/listinfo/bioconductor
>> Search the archives:
>> http://news.gmane.org/gmane.science.biology.informatics.conductor
>
>
--
shan gao
Room 231(Dr.Fei lab)
Boyce Thompson Institute
Cornell University
Tower Road, Ithaca, NY 14853-1801
Office phone: 1-607-254-1267(day)
Official email:sg839 at cornell.edu
Facebook:http://www.facebook.com/profile.php?id=100001986532253
More information about the Bioconductor
mailing list