Error estimating covariance matrix and individual admixture proportions

Hi, I am trying to perform a PCA with some strains whole-genome NGS data. When i run the command: python pcangsd.py -beagle name.beagle.gz -admix -o name_test -threads 10, I obtain "ValueError: k must be greater than 0." Can someone help me understand what I am doing wrong? Thanks in advance!

pca pcangsd

• 847 views • link

updated 6.2 years ago by Sej Modha

• written 6.3 years ago by carolina.santiago.t •

0 answers

No answers yet.

More posts like this

How is module-trait association calculated in WGCNA?
written 3.4 years ago by kivanvan •

I have read some papers using WGCNA to explore relationship between metabolomics and clinical outcomes. I wanted to use this method in my work as …
Calling star alleles with PGx-POP on my VCF file
written 4.0 years ago by Vanish007 •

Hi everyone, I am currently trying to have star alleles called on my phased vcf file. So far I have a vcf file with a …
FastaAlternateReferenceMaker is not outputting fasta sequence for multiple samples from a VCF
written 4.1 years ago by anikcropscience •

Hello, I have the coordinates for a gene of interest and I am trying to obtain a consensus sequence from a multi-sample `VCF` file where …
Ensembl vep. How to filter population frequency less than 1%?
written 4.8 years ago by Gonzalo •

Hi everyone, I have gotten many responses from the site even though I never asked a question, this is my first query. I am working …
somatic variant calling with Mutect2
written 5.0 years ago by Elisa •

Hi everybody, I am performing somatic variant calling on NGS data (NGS targeted sequencing performed with Illumina MiSeq) using Mutect2. I do not have the …
VEP is very slow. Fork doesnt seem to work
written 5.4 years ago by nhaus •

I am using VEP (v103) to annotate a small vcf file (~1000 variants). Nonetheless, it takes very long (>20 minutes) , which doesnt quite match …
How to change MirBAse ID with org.Hs.eg.db library
written 6.1 years ago by Van9 •

I recently started learning bioinformatics for my final thesis and I am doing MiRNA analysis in R. I am trying to obtain the ensembl id …
Cannot run HaplotypeCaller in GATK
written 6.4 years ago by ImBioMario •

Hello! I am trying to conduct a SNP calling with HaplotypeCaller but i got an error message that theres something wrong with my *.BAM files. …
Tophat Usage difficulty?
written 6.9 years ago by jaqx008

Hello All, I am having difficulty understanding and executing the tophat usage after looking up the help. its just not as simple as running other …
The problem to run kmergenie
written 11.4 years ago by ziyingliu1996 •

I am new to the RNA Seq assembly. I am trying to use velvet-Oases to do the work. To be able to determining the best …