Question: Error running smartpca from eigensoft
0
gravatar for scho249
17 months ago by
scho24920
scho24920 wrote:

Hi,

I have encountered an error trying to run smartpca from eigensoft.

I have converted ped file of my SNPs to eigensoft file format using fcgene and entered the following command line

to run smartpca:

 perl ./bin/smartpca.perl -i out_eigensoft.ped -a out_eigensoft.pedsnp -b out_eigensoft.pedind -m 0 -o out.pca -p out.plot -l out.log -e out.eval 

 

When I ran it, it gives an error claiming that:

OOPS couldn't open file out.pca.evec for reading at ./bin/smartpca.perl line 65.

 

All the files are in the bin folder, so I believe this is not the path problem...

I have been stuck for days now. Could someone help me please?

 

Shane 

 

snp software error • 1.7k views
ADD COMMENTlink written 17 months ago by scho24920
1

Could you paste the whole content of "out.log" to pastebin.com or here, if it's not too big? smartpca.perl runs several steps internally and I believe one of them dies, an early step creates out.pca.evec and a later step depends on out.pca.evec and doesn't find it

ADD REPLYlink written 17 months ago by Philipp Bayer4.3k

Hi Philipp,

Thanks and sorry for the late reply.

There isn't any output in the "out.log" file but only in "out.pca.par" file and this is the content it output:

genotypename: /gpfs1m/projects/landcare00051/eigensoft/dart/all/out_eigensoft.ped
snpname: /gpfs1m/projects/landcare00051/eigensoft/dart/all/out_eigensoft.pedsnp
indivname: /gpfs1m/projects/landcare00051/eigensoft/dart/all/out_eigensoft.pedind
evecoutname: /gpfs1m/projects/landcare00051/eigensoft/dart/all/pca/out.pca.evec
evaloutname: /gpfs1m/projects/landcare00051/eigensoft/dart/all/pca/out.eval
altnormstyle: NO
numoutevec: 10
numoutlieriter: 0
numoutlierevec: 10
outliersigmathresh: 6
qtmode: 0

 

I have set the PATH variable as well, so I was thinking it shouldn't have any problem locating the binaries....

 

ADD REPLYlink written 17 months ago by scho24920
1

Have you tried using smartPCA.perl with the example data provided in the archive (run "perl example.perl" in the EIGENSTRAT folder)? That way you'll at least know that PATH is correctly set up

If it works, you may run the steps that smartpca.perl runs manually, these are the steps it runs:

smartpca -p your_par_file

ploteig -i your_evec_file -c 1:2 -p your_populations -x -y -o your_plot_file

evec2pca.perl your_principal_components_number(10) your_evec_file your_individuals_file your_pca_file
ADD REPLYlink modified 17 months ago • written 17 months ago by Philipp Bayer4.3k
1

Hi Philipp,

Thank you so much for your advice. The example wasn't running as well and it turns out that the binary file 'smartpca' was not executable. So I used 'sudo chmod ug+x' and now it works fine.

Again thank you so much for your help!

Shane

ADD REPLYlink written 16 months ago by scho24920
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 666 users visited in the last hour