Generate a state transition template matrix based on numbers of alleles at QTLs.
generateRateTransTemplate(numsAllelesAtSites, nameRows = FALSE)
numsAllelesAtSites | integer vector specifying the number of alleles at each site, the length of the vector defining the number of QTLs. |
---|---|
nameRows | logical indicating whether the returned matrix should have named rows (for explicative purpose only, default FALSE). The column names are different for each row, i.e. the row-names excluding the genotype corresponding to that specific row. |
each element of the matrix is a multiplier for the corresponding genotype transition rate (see function simulateEpidemic in the patherit package). This multiplier is nonzero only for transitions that need exactly one mutation (haming distance = 1) in order to convert one genotype into another. The