Differences

This shows you the differences between two versions of the page.

--- ima2 [2012/05/06 15:01] – heidi
+++ ima2 [2013/05/03 10:34] (current) – heidi
@@ Line 21: / Line 21: @@
   * DNA
   * Microsatellites (STR)
@@ Line 36: / Line 37: @@
   * line 2 - number of populations (npops
   * line 3 - population names in order, seperated by one or more space. This order also corresponds to the order in which the populations are numbered in the population tree and the order in which the data occur for each locus.
-  * line 4 - the population string in modified Newick format. The string contains information on the topology of the tree for the sampled populations and information on the ordering of the internal nodes in time.  These
+  * line 4 - the population string in modified Newick format. The string contains information on the topology of the tree for the sampled populations and information on the ordering of the internal nodes in time.  These internal nodes correspond to ancestral populations. The ancestral populations are numbered beginning with npops for the most recent ancestral population and proceeds up to 2×(npops-1) for the ancestor of all the sampled populations. Sampled populations in the string are represented by their respective number. Ancestral populations are represented by a colon, i.e. ‘:’, followed by their ancestral population number.
-internal nodes correspond to ancestral populations. The ancestral populations are numbered beginning with npops for the most recent ancestral population and proceeds up to 2×(npops-1) for the ancestor of all the sampled populations. Sampled populations in the string are represented by their respective number. Ancestral populations are represented by a colon, i.e. ‘:’, followed by their ancestral population number.
     * If there is only a single population then the tree string is simply: 0.
     * If there are two populations then the tree string is: (0,1):2
@@ Line 49: / Line 49: @@
     - The mutation rate per year for the locus (not per base pair). This can be left blank, but is needed for estimating parameters on demographic scales. If there are multiple STRs in the locus then there can be multiple mutation rates on this line separated by spaces. If the locus is a HapSTR, then the first mutation rate given applies to the sequence portion of the locus with subsequent values corresponding to STR markers included in the locus.
     - If the mutation rate is given, it can be followed by a range of mutation rates that can be used (with ranges for other loci in the analysis) to set priors on the ratios of mutation rate scalars. The range is entered with an open parentheses, the lowest value, a comma, the highest value, and a closed parentheses (e.g. ‘(0.00001, 0.00004)’. The range must bracket the rate. For a locus with multiple mutation rates, and multiple ranges, each range follows its corresponding mutation rate immediately on line.
-  * line 7 - data for gene copy # 1 from population 0. The first 10 spaces are devoted to the sample name. The sequence or allele length (for SSM model) begins in column 11 of the file. The sequence for a given sample is given all on one line without gaps. For SSM or HapSTR data, the allele length assumes a step size of 1. This means that data from STRs that are multiples of lengths greater than 1 must be converted to counts of the number of base repeats (e.g. for a dinucloitide ‘CACACACACACA’ the length would be 6). If the data is for an SSM model locus and there are multiple STRs, then there will be one integer on each line for each STR, separated by a space. If the locus is HapSTR (joint IS and SSM) then the STR data is given on the line, beginning at column 11, followed by the sequence data. For SSM data, as for other types of data, only one gene copy is represented on each line of the data file. Diploid genotype data must be broken up and listed, with one data line for each gene copy.
+  * line 7 - data for gene copy # 1 from population 0. The first 10 spaces are devoted to the sample name. The sequence or allele length (for SSM model) begins in column 11 of the file. The sequence for a given sample is given all on one line without gaps. For SSM or HapSTR data, the allele length assumes a step size of 1. This means that data from STRs that are multiples of lengths greater than 1 must be converted to counts of the number of base repeats (e.g. for a dinucloitide ‘CACACACACACA’ the length would be 6). Any number less than 5 causes the program to stop with an error. If the data is for an SSM model locus and there are multiple STRs, then there will be one integer on each line for each STR, separated by a space. If the locus is HapSTR (joint IS and SSM) then the STR data is given on the line, beginning at column 11, followed by the sequence data. For SSM data, as for other types of data, only one gene copy is represented on each line of the data file. Diploid genotype data must be broken up and listed, with one data line for each gene copy.
   * lines 8 thru line  - the remainder of the data for locus 1. Each line contains the data for one sample.  The data for locus 1 for population 1 immediately follow those for population 0, and so on
   * Additional lines for additional loci.  Each locus begins with a line containing the information for that locus,  in the same format as for the first locus. The sample names and sample sizes  for additional loci and the inheritance scalars and mutation model for additional loci do not have to be the same as for locus 1 (generally they are not).