我需要从BioConductor的Shortread库中创建一个类型的Shortreadq对象。

ShortReadQ 'signature(sread = "DNAStringSet", quality =
          "QualityScore", id = "BStringSet")'

质量插槽需要是一个从质量分离的对象,我可以轻松地从我希望模仿的另一个Shortreadq对象中确定。

> class(quality(anotherObject))
[1] "SFastqQuality"
attr(,"package")
[1] "ShortRead"

在Cromentor参数中,使用该信息(“ SfastQQuelaligal”)的最佳方法是什么?

newObject<-ShortReadQ(sread=...,
             quality=SFastqQuality(...), 
             id=...)
有帮助吗?

解决方案 3

感谢您的回复。他们带领我找到了有效的解决方案

newObject<-ShortReadQ(sread=...,
             quality=new(Class=class(quality(anotherObject)),theFirstParameter=...), 
             id=...)

其他提示

这可以做你想要的吗?

quality = new(class(old.quality.obj)[[1]]))

您可能需要获得get函数:

a <- get(class(object))
a(...)
许可以下: CC-BY-SA归因
不隶属于 StackOverflow
scroll top