Error with smarteigenstrat.perl accepting pca output from smartpca.perl
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
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- 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