問題タブ [genomicranges]

何千ものスコア列を持つ GRanges オブジェクト ( tomap ) と、関心領域がありメタデータがない ( roi ) オブジェクトがあります。tomapの各列の最大スコアをroiの対応する間隔にマップしようとしています。

また、スコア列の名前も保持したいと思います (私の実際のデータでは、これらは意味のある名前であり、score1、score2 などのように一般化できません)。特定の列に対して行うことはできますが、すべての列に一般化するのに苦労しています。

これが私がこれまでに得たものです：

ご覧のとおり、各スコア列を個別に指定するとこれが機能しますが、何千もの列に対してこれを行うにはどうすればよいでしょうか?

この簡単なコード:

やや不可解な例外で失敗します:

validObject(.Object) のエラー: 無効なクラス “GRanges” オブジェクト: 1: 'x@seqnames' は 'x' と平行ではありません 無効なクラス “GRanges” オブジェクト: 2: 'x@strand' は 'x' と平行ではありません

さらに、提供しようとするとseqlengths：

私は得る：

.normargSeqlengths(seqlengths, seqnames) のエラー: 指定された 'seqlengths' の長さは、シーケンスの数と等しくなければなりません

これは、プロセスで一部のデータが削除されたことを示しています。しかし、なぜそれが起こるのか理解できません。

ここで何が起こっているのかについての洞察に感謝します。

環境：