I'll give Kathryn a copy of my (WinBUGS) solution to the bonus question later today or tomorrow. So she should be able to give you some more hints after that. In the meantime, I've posted a WinBUGS-friendly version of the data to the course web site (data.txt). Here are the variables you'll find in the data set:
- n.one -- The number of replicate populations produced in the one-generation experiment (100).
- n.one.min -- The mininum possible number of alleles in the population (because the sample contains at least that many copies, 34).
- n.time -- The number of generations for the time series (50).
- n.time.min -- The minimum possible number of alleles in any of the populations (because the sample contains at least that many copies, 52).
- one -- The replicate population data for the one-generation experiment.
- time -- The data for the time-series experiment.
Oh, one hint now. If the numbers you calculate directly for effective size are different from those you get from WinBUGS, don't worry. You should get a number fairly similar in the one-generation experiment, but the number in the time series experiment may differ by a factor of two.
Leave a comment