annotate_assembly {ballgown} | R Documentation |
match assembled transcripts to annotated transcripts
annotate_assembly(assembled, annotated)
assembled |
|
annotated |
|
If gown
is a ballgown
object, assembled
can be
structure(gown)$trans
(or any subset). You can generate a
GRangesList
object
containing annotated transcripts from a gtf file using the
gffReadGR
function and
setting splitByTranscripts=TRUE
.
data frame, where each row contains assembledInd
and
annotatedInd
(indexes of overlapping transcripts in assembled
and annotated
), and the percent overlap between the two transcripts.
Alyssa Frazee
data(bg) gtfPath = system.file('extdata', 'annot.gtf.gz', package='ballgown') annot = gffReadGR(gtfPath, splitByTranscript=TRUE) info = annotate_assembly(assembled=structure(bg)$trans, annotated=annot)