Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

0 % of model's snps used for trans QTL #189

Open
elnetwas opened this issue Jan 19, 2024 · 0 comments
Open

0 % of model's snps used for trans QTL #189

elnetwas opened this issue Jan 19, 2024 · 0 comments

Comments

@elnetwas
Copy link

Hello MetaXcan,
I am attempting to use predict.py and receiving the INFO - 0 % of models' snps used message. My model in db format is similar to that of elastic net, obtained using the tensorQTL tool, where I perform a trans-mapping. I format these QTLs appropriately for elastic net with 57,997 variants, resulting in the following format:
image

Considering that the gene variable is predetermined and always set as ENSG0000269981, as it pertains to a single phenotype, and my imputed genotype (which is the same one used in tensorQTL to obtain the above-mentioned variants) is in the following format:
image

I apply the following script using Predict.py:
python3 $METAXCAN/Predict.py --model_db_path trans_associations_mucosa.db --model_db_snp_key varID --vcf_genotypes /mnt/lustre/scratch/nlsas/home/otras/fmx/fgs/carmen/telomeros/GENOTIPO/resultados/GENOTIPOS_IMPUTADOS_FILTRADO_FASTQTL.vcf.gz --vcf_mode genotyped --prediction_output trans_associations_mucosa_umbral_dos_predict.txt --prediction_summary_output trans_associations_mucosa_umbral_dos_summary.txt --verbosity 9 --throw --on_the_fly_mapping METADATA "{}{}{}_{}_b37"
However, I consistently get "0 % of models' snps used." I have several questions: Can predict.py be applied to trans-QTLs? Or is it possible that the low number of variants is causing the result of "0% of models' snps used"?

Thank you for your time!!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant