DReichLab / DReichLab/EIG

Error with smarteigenstrat.perl accepting pca output from smartpca.perl

Open
#75 4 comments 0 reactions 0 assignees View on GitHub

Nobody has claimed this yet.

Dominant language
C
Stars
207
Forks
61
PR merge metrics
No merged PRs in 30d

Description

kuretlab@Synuclein:~/Bayesian/EIG/bin$ $ ./smarteigenstrat.perl -output.geno -a output.snp -b output.ind -p output.evec -o output.chiseq -l output.log -k 10
smarteigenstrat -p output.chiseq.par >output.log
OOPS: output.evec contains too many entries

Hello,
I ran smartpca as instructed on my WGS dataset (349 individuals, 20 million SNPs). The smartpca.perl wrapper ran fine producing 10 principal components in an .evec file. However, when loading the pca file into smarteigenstrat.perl, the aforementioned error was returned. What I find interesting is that after outlier removal by smartpca, the .ind file had 349 individuals even when sigma was at default value and the log file noted the removal of outliers (around 30 individuals). Is there something I can do to fix this. Thanks in advance

Contributor guide

No contributing guide indexed for this repository

First steps

  1. Read the whole issue, then the project's contributing guide.
  2. Comment on the issue to say you are picking it up — it saves two people doing the same work.
  3. Fork the repository and make your change on a branch.
  4. Open a pull request that references the issue number.

Research direction

Start by reproducing the shown smarteigenstrat.perl command with output.evec, output.ind, and output.chiseq.par, then inspect how smartpca.perl produces output.evec and how smarteigenstrat.perl reads it. Compare the .evec entry count with the outlier-removal information in the log; done means the PCA output is accepted or the incompatibility and required input format are clearly documented.

Written by the indexing model from the issue text.

Assessment

Tech stack
perl
Domain
data
Issue type
Bug
Difficulty
4/5
Estimated time
3-5 days
Activity status
Stale
Clarity
Mostly clear
Newbie friendliness
35/100

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.