I wonder if it is worth repeating the description of the "pattern" variable from the original PSMC paper. It says:
Given a maximum TMRCA in the 2N0 scale of Tmax, and a number of atomic time intervals n, let the boundaries of these intervals be ti = 0.1exp[i/n log(1 + 10Tmax)]−0.1, i = 0,…, n. To reduce the complexity of the search space, blocks of adjacent atomic intervals were combined to have the same population-size parameter via a user-specified pattern. On autosome and simulated data, Tmax = 15, n = 64 and the pattern is ‘14 + 252 + 14 + 16’, which means that the first population-size parameter spans the first four atomic time intervals, each of the next 25 parameters spans two intervals, the twenty-seventh parameter spans four intervals and the last parameter spans the last six time intervals. On X-chromosome data, Tmax = 15, n = 60 and the pattern is ‘16 + 24 + 13 + 132 + 13 + 24 + 1*6’.
I presume they wouldn't mind you simply quoting that in the documentation?
I wonder if it is worth repeating the description of the "pattern" variable from the original PSMC paper. It says:
I presume they wouldn't mind you simply quoting that in the documentation?