1

mclapply は CompressedRleList で動作しますか?

たとえば、葉巻 (a) のベクトルがあり、cigarToRleList はそのベクトルを指定して CompressedRleList を返します。

a = c("44M","44M","9S35M","44M","40M4S")
b = cigarToRleList(a)

b で mclapply を実行しようとすると、次のエラーが発生します。

mclapply(b,function(x) class(x))
Error in as.list.default(X) : 
  no method for coercing this S4 class to a vector

このリストの長さは 1 であるため、 b をリストに強制しても役に立たないようです。

mclapply(list(b),function(x) class(x))
[[1]]
[1] "CompressedRleList"
attr(,"package")
[1] "IRanges"

mclapply を CompressedRleList 構造で動作させる方法はありますか?

4

1 に答える 1

2

これはバイオコンダクターに関する質問なので、バイオコンダクターのメーリングリストに問い合わせてください。過去に問題があったため、の出力を提供しますsessionInfo()。私は現在のBiocでこれがうまくいくと信じています。それは

> library(GenomicRanges)
> library(parallel)
> a = c("44M","44M","9S35M","44M","40M4S")
> b = cigarToRleList(a)
> unlist(mclapply(b, class))
[1] "Rle" "Rle" "Rle" "Rle" "Rle"
> sessionInfo()
R version 3.0.2 RC (2013-09-23 r63972)
Platform: x86_64-unknown-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=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

attached base packages:
[1] parallel  stats     graphics  grDevices utils     datasets  methods  
[8] base     

other attached packages:
[1] GenomicRanges_1.12.5 IRanges_1.18.3       BiocGenerics_0.6.0  

loaded via a namespace (and not attached):
[1] stats4_3.0.2


library(GenomicRanges)
library(parallel)
> unlist(mclapply(b, class))
[1] "Rle" "Rle" "Rle" "Rle" "Rle"
> sessionInfo()
R version 3.0.2 RC (2013-09-23 r63972)
Platform: x86_64-unknown-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=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

attached base packages:
[1] parallel  stats     graphics  grDevices utils     datasets  methods  
[8] base     

other attached packages:
[1] GenomicRanges_1.12.5 IRanges_1.18.3       BiocGenerics_0.6.0  

loaded via a namespace (and not attached):
[1] stats4_3.0.2
于 2013-09-27T18:01:37.790 に答える