Folded Multidimensional Site Frequency Spectrum Format

240 views
Skip to first unread message

qspreng...@gmail.com

unread,
Feb 7, 2018, 2:22:57 PM2/7/18
to fastsimcoal

I have a question about the format for a folded multidimensional site frequency spectrum using the minor allele.  Our current input files (attached is a 2pop and 3pop SFS) have zeros in places where there are implausible events. Is it that all right? 

 

Also, when using a folded MSFS, in the command line would I use the  --foldedSFS option? Even for three or more populations? 

 

 e.g. ./fsc26 -t 3POPChrX.tpl -n 100000 -m -e 3POPChrX.est -M 0.001  -L 200 -q --multiSFS -c3 -B3 --foldedSFS    

Thanks!!
Quentin 
2POP_MSFS.obs
3POP_MSFS.obs

Laurent Excoffier

unread,
Apr 23, 2018, 5:30:43 PM4/23/18
to fastsimcoal
No for a folded sfs just use the -m option instead of the -d option. The --foldedSFS is to be used when the folded sfs has been computed by angsd for multiple pairs of populations 

Laurent Excoffier

unread,
Apr 23, 2018, 5:31:38 PM4/23/18
to fastsimcoal
also note that the -M 0.001 syntax is obsolete for fsc26 and can lead to problems, just use -M
Laurent
Reply all
Reply to author
Forward
0 new messages