Question: (Closed) filter based on variant annotation
gravatar for winstondunnmd
12 weeks ago by
winstondunnmd0 wrote:

Dear Bioconductors:

I am doing SKAT test of genes based on their position on the chromosome. I use the following TopmedPipeline function and codes. I hope to additionally filter based on variant annotation. I have identified missense_variants and recorded their in a vector. Is it possible to additionally filter based on a vector that contains of the missense_variants?

Many thanks!


aggVarList <- aggregateGRanges(aggByPos)
seqResetFilter(seqData, verbose=FALSE)
iterator <- SeqVarRangeIterator(seqData, variantRanges=aggVarList, verbose=FALSE)
assoc <- assocTestAggregate(iterator, nullmod, test="SKAT", AF.max=0.1, weight.beta=c(1,25))
ADD COMMENTlink written 12 weeks ago by winstondunnmd0

Hello winstondunnmd!

We believe that this post does not fit the main topic of this site.

TopMedPipeline isn't a Bioconductor package. You could try or

For this reason we have closed your question. This allows us to keep the site focused on the topics that the community can help with.

If you disagree please tell us why in a reply below, we'll be happy to talk about it.


ADD REPLYlink written 12 weeks ago by James W. MacDonald51k
Please log in to add an answer.
The thread is closed. No new answers may be added.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 333 users visited in the last hour