Skip to content

Documentation on pattern #5

@hyanwong

Description

@hyanwong

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?

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions