I’ve just uploaded Project #2. As usual, you can find it either from the Lecture Schedule page or from the direct link below. This project is different from any of the lab exercises you’ve done so far. There isn’t data to analyze, and there aren’t simulations to do. Instead, there’s a paper by Robin Waples and David Teel to read: Conservation Genetics of Pacific Salmon I. Temporal Changes in Allele Frequency. After reading the paper, I have five questions for you to answer. When I grade this project, I will be evaluating how well you use what you’ve learned about genetic drift and natural selection to answer the questions. None of the answers need to be more than a couple of paragraphs. Feel free to submit them in whatever form you find convenient, Word document, R notebook, PDF, or Pages are formats I know I can handle easily. If you send it in a form I don’t recognize, I’ll be in touch.
I’ve just posted Lab 7. You’ll find it in the usual places, i.e., from the Lab Schedule page or from the direct link below. As you’ll see, this lab exercise will allow you to explore the coalescent by comparing the time to coalescence of all alleles in a sample as a function of the rate of migration, the number of populations exchanging genes, and the migration model, i.e., either the finite island model you explored in last week’s exercise or the one-dimensional stepping stone. As noted in the exercise itself, the run_simulation() function will run 1000 samples by default. I recommend picking a smaller number first, say 50 or 100, to get a sense of how long a full run will take before you start the simulation. It shouldn’t be surprising that simulations take longer the more populations that you specify. There’s no need to try simulating more than 100 populations. You’ll see any patterns associated with differing numbers of populations by the time you get to that many.
I’ll have grades returned for Week 5’s lab by the end of the day tomorrow. This week’s lab is a bit more involved for a couple of reasons.
- Each set of simulation conditions will take you longer to run. Each set takes about half an hour on my MacBook. You’ll want to get an early start on the exercise so you have a chance to explore at least five different sets of conditions.
- The questions ask you to stretch a little and interpret the results of the simulations relative to theoretical expectations.
As usual, you can get directly to the exercise from the link below, or you can get to it through the Lab Schedule page.
I just posted the lab exercise for week 5. You can find it from the Lab Schedule page, or you can click on the link below. As you’ll see, this week’s exercise is a bit different. Rather than analyzing data, you’ll be running a small simulation to explore some properties of drift. Provided that you show results from at least 10 different sets of parameters, you’ll get 10 out of 10 points for this exercise. The next two weeks will be similar, except that I’m liable to ask you to stretch a bit and explain some aspect of the results you see in weeks 6 and 7.
It didn’t take me as long to add notes on fertility selection as I thought it might. I found some old notes that I could borrow from. You can find links to the new versions on the lecture detail page or you can click on the links below to go directly to them. The fertility selection part of the notes comes at the end.
Project 1 is now available. You can find it from the Lab Schedule, or you can download it directly from this link. You’ll be analyzing the same data as you did last week, but it’s in a different format. If you continue to use population genetics in your research after this course, one of the things you’ll discover is that there are a variety of different formats used by different packages. In addition to using R for analyses, you’ll find that you need to learn how to wrangle data in R (or in Python), but I’ll do my best to protect you from that this semester.
I included a few details about the Bayesian model that Hickory uses that weren’t in the notes I posted earlier – the among population allele frequency distribution and extending the model to includ locus- and population-specific effects. I’ve edited the notes to include those details, and you’ll find new HTML and PDF versions on the website. They’ll appear in the consolidated “book” version of the notes in mid-late December.
In addition to the additional reading on F-statistics that I already provided in the lecture detail associated with today’s lecture, I just added a link to an R notebook that delves a bit more into how Bayesian inference is typically implemented and that illustrates a bit better how the likelihood, prior, and posterior are related to one another. We’ll spend some time exploring these ideas at the start of today’s lecture
I just posted the notes and data for Lab #3. You’ll find them in the Lab Schedule as usual. As I’ll explain in lecture, we’re going to use LEA instead of STRUCTURE for this exercise. STRUCTURE has some significant advantages over LEA, but LEA runs much faster. More importantly, I haven’t been able to get STRUCTURE to run on my MacBook since I upgraded it to Big Sur. If you have a Mac and you’d like to try getting STRUCTURE to run, I’d be delighted if you figured out how. But even if you do figure it out, I recommend that you use LEA for this project. Analyzing the data with STRUCTURE would take a very long time.
Next week Project #1 will build on the analysis you do this week with LEA to compare it with a different method (DAPC) and to address some of the questions that the authors of the paper from which these data were extracted.
Many of you submitted lab assignments in a format that made it clear that you already know about R notebooks. We’ll spend a little time discussing them in class this morning, including the best way for you to send me an R notebook from a lab assignment. In the meantime, here’s a link to a quick introduction to R notebooks (and another to a blog post I made several years ago after I belatedly discovered how useful they are). I’ll also add the quick introduction link to the lab resources for this week’s lab.