Rare variant meta-analyzing tests · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

Rare variant meta-analyzing tests

Hi all,

I have 5 different (dosage VCF files, with both continuous and binary phenotypes) datasets and I'm trying to meta-analyze rare variants in each dataset. Ideally I would like to group the variants per gene and then meta-analyze the p-values per gene from all 5 datasets, but I haven't found any software that easily does that. I played around a bit with rare metal and rvtests, but so far I keep on running into errors and weird results. Does anyone has any tips, ideas or tricks how to best do this? Any help is really appreciated!

meta analyze rare variants

• 1,690 views • link

updated 7.2 years ago by swetansu •

• written 8.5 years ago by Floris Brenk

1 answer

Have you tried SKAT/SKATO? R package SKAT has the details!! Worth exploring custom weighing scheme besides default MAF based weighing..

• 0 views • link

written 7.2 years ago by swetansu •

Log in to answer this question.

More posts like this

Extracting LOF variants and incidence from GnomAD for a list of genes
written 2.6 years ago by hello_bioinf •

Hi all, I have a list of around 200 genes which I would like to observe in GnomAD v3.2.1 I want to see the general …
rvtests : all genes are "Gene xxx has 0 variants, skipping"
written 4.6 years ago by Pierre Lindenbaum

Hi, this is the first time I'm using rvtests https://github.com/zhanxw/rvtests (version: 20190205, git: c86e589efef15382603300dc7f4c3394c82d69b8) the log file (including the command line) is: ``` Effective Options …
SKATO hangs after going through several genes
written 4.7 years ago by @anna •

Hi there I'm doing skato analysis with rvtest/2.1.0. My script: rvtest --inVcf file.vcf.gz --pheno phenofile --out newfile --geneFile refFlat_hg19.txt --kernel skato it runs for a …
How to conduct germline sharing WES analysis of three probands in a large family with the same tumo…
written 6.4 years ago by Maverick77 •

I have a question regarding sharing. Is there a 'sharing of the same putitative genes across 3 probands' type of tutorial or pipeline you could …
Add gene annotation/information to CNV plot in R
written 8.7 years ago by Floris Brenk

Hi all, Im have a few interesting (genotyping) CNVs and now I want to plot them using the B allele frequency and LogR ratio. There …
Merge several association signals in a region and use as a phenotype
written 9.1 years ago by Floris Brenk

Dear all, I got some case/control association signals in one locus and now want to perform a bit more advanced eQTL analysis in this region. …
RNA-seq de novo transcript assembly from one gene
written 10.2 years ago by Floris Brenk

Hi all, I was wondering whether anyone has experience with de novo transcript assembly of RNA-seq (100 bp PE Illumina reads) of only one gene. …
Problem with visualization of Haloplex exomes in IGV, can not find variant back bam file
written 11.1 years ago by Floris Brenk

Hi all, I'm analysing a few Haloplex exomes and after our standard GATK pipeline (not removing duplicates, due to haloplex exomes) I have a number …
Is it possible to perform meta-analysis of gene-based LMM tests for a binary trait?
written 11.4 years ago by Nick •

Hi everyone, We have exome-chip data (disease cases and healthy controls) from three different studies. Currently, for each study I have performed group-based tests (CMC …
Cumulative Effect Of Eqtl Snps Testing
written 12.5 years ago by Floris Brenk

<p>Hi all,</p> <p>I recently performed a basic eQTL analysis and now want to do some basic followup analysis. So first of all does anyone has …