我需要从 Bioconductor 的 ShortRead 库中创建一个 ShortReadQ 类型的对象。
ShortReadQ 'signature(sread = "DNAStringSet", quality =
"QualityScore", id = "BStringSet")'
质量槽需要是从 QualityScore 继承的对象,我可以很容易地从我希望模拟的另一个 ShortReadQ 对象中确定它。
> class(quality(anotherObject))
[1] "SFastqQuality"
attr(,"package")
[1] "ShortRead"
在构造函数参数中使用该信息(“SFastqQuality”)的最佳方式是什么?
newObject<-ShortReadQ(sread=...,
quality=SFastqQuality(...),
id=...)