[adegenet-forum] long processing with fasta2bin

Adrien Bessy adrienbessy at hotmail.com
Wed Dec 2 13:06:21 CET 2015




Hi there,
I would like to calculate Fst on 1292 genes between 5 populations.So I have 1292 fasta files containing around 30 individual sequences.I saw there is a pairwise.fst function that could be what I want. 
But before that I need to use the fasta2bin function in order to read alignments form my fasta files.The problem is that when I apply this function, for a fasta file containing 30 alignments, the processing seems to be very very long.
Here are my script and the results:--------------------------------------------------------------------------------library(adegenet)res <- fasta2DNAbin("/Users/abessy/Documents/SA_without_N/DSPR_without_N/140up_Chr3R_14124163_14125548.fas_without_N.fas")
> source("fasta.R")
Loading required package: ade4


   /// adegenet 2.0.1 is loaded ////////////


   > overview: '?adegenet'
   > tutorials/doc/questions: 'adegenetWeb()' 
   > bug reports/feature resquests: adegenetIssues()






 Converting FASTA alignment into a DNAbin object... 




 Finding the size of a single genome... 



-----------------------------------------------------------------------------------
I tried with the fasta2genlight function and the processing seems to be very very long too, even with using the parallel parameter:-----------------------------------------------------------------------------------library(adegenet)library(parallel)res1<- fasta2genlight("/Users/abessy/Documents/SA_without_N/DSPR_without_N/140up_Chr3R_14124163_14125548.fas_without_N.fas",parallel=TRUE)
> source("fasta.R")
Loading required package: ade4


   /// adegenet 2.0.1 is loaded ////////////


   > overview: '?adegenet'
   > tutorials/doc/questions: 'adegenetWeb()' 
   > bug reports/feature resquests: adegenetIssues()






 Converting FASTA alignment into a genlight object... 


---------------------------------------------------
Have you an idea why it is so long ?I have a iMac. 
I attached my fasta file.
Best wishes
Adrien


 		 	   		  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/adegenet-forum/attachments/20151202/d274b399/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 140up_Chr3R_14124163_14125548_fas_without_N_DSPR.fasta
Type: application/octet-stream
Size: 20924 bytes
Desc: not available
URL: <http://lists.r-forge.r-project.org/pipermail/adegenet-forum/attachments/20151202/d274b399/attachment-0001.obj>


More information about the adegenet-forum mailing list