Results

You’ll find the results in three separate Excel workbooks:

Running the analyses

If you’re interested in how I ran the analyses, you’ll find the code in GWAS.r. Just put the code in the same directory as gypsymoth.csv and type source("GWAS.r") at the command prompt in an R console. It will churn for a day or a day and a half and put three CSV files in your working directory. The Excel workbooks come from opening the CSV files and adding a little formatting.

If I’d had a bit more time, I could have written a version of GWAS.r that ran the analyses much more quickly, but getting that right would have taken more time than simply running the crude version linked here.

Genomic prediction

If you want to see what a crude approach to genomic prediction looks like with the gypsymoth data, refer to genomic-prediction.nb.html (Rmd version)

LS0tCnRpdGxlOiAiR1dBUyBvbiBneXBzeSBtb3RocyByZXN1bHRzIgpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sKLS0tCgojIFJlc3VsdHMKCllvdSdsbCBmaW5kIHRoZSByZXN1bHRzIGluIHRocmVlIHNlcGFyYXRlIEV4Y2VsIHdvcmtib29rczoKCi0gW0dXQVMtTWFzcy1yZXN1bHRzLnhsc3hdKGh0dHA6Ly9kYXJ3aW4uZWViLnVjb25uLmVkdS9lZWIzNDgtcmVzb3VyY2VzL0dXQVMtTWFzcy1yZXN1bHRzLnhsc3gpCi0gW0dXQVMtUEQtcmVzdWx0cy54bHN4XShodHRwOi8vZGFyd2luLmVlYi51Y29ubi5lZHUvZWViMzQ4LXJlc291cmNlcy9HV0FTLVBELXJlc3VsdHMueGxzeCkKLSBbR1dBUy1URFQtcmVzdWx0cy54bHN4XShodHRwOi8vZGFyd2luLmVlYi51Y29ubi5lZHUvZWViMzQ4LXJlc291cmNlcy9HV0FTLVREVC1yZXN1bHRzLnhsc3gpCgojIFJ1bm5pbmcgdGhlIGFuYWx5c2VzCgpJZiB5b3UncmUgaW50ZXJlc3RlZCBpbiBob3cgSSByYW4gdGhlIGFuYWx5c2VzLCB5b3UnbGwgZmluZCB0aGUgY29kZSBpbiBbR1dBUy5yXShodHRwOi8vZGFyd2luLmVlYi51Y29ubi5lZHUvZWViMzQ4LXJlc291cmNlcy9HV0FTLnIpLiBKdXN0IHB1dCB0aGUgY29kZSBpbiB0aGUgc2FtZSBkaXJlY3RvcnkgYXMgYGd5cHN5bW90aC5jc3ZgIGFuZCB0eXBlIGBzb3VyY2UoIkdXQVMuciIpYCBhdCB0aGUgY29tbWFuZCBwcm9tcHQgaW4gYW4gYFJgIGNvbnNvbGUuIEl0IHdpbGwgY2h1cm4gZm9yIGEgZGF5IG9yIGEgZGF5IGFuZCBhIGhhbGYgYW5kIHB1dCB0aHJlZSBDU1YgZmlsZXMgaW4geW91ciB3b3JraW5nIGRpcmVjdG9yeS4gVGhlIEV4Y2VsIHdvcmtib29rcyBjb21lIGZyb20gb3BlbmluZyB0aGUgQ1NWIGZpbGVzIGFuZCBhZGRpbmcgYSBsaXR0bGUgZm9ybWF0dGluZy4KCklmIEknZCBoYWQgYSBiaXQgbW9yZSB0aW1lLCBJIGNvdWxkIGhhdmUgd3JpdHRlbiBhIHZlcnNpb24gb2YgYEdXQVMucmAgdGhhdCByYW4gdGhlIGFuYWx5c2VzIG11Y2ggbW9yZSBxdWlja2x5LCBidXQgZ2V0dGluZyB0aGF0IHJpZ2h0IHdvdWxkIGhhdmUgdGFrZW4gbW9yZSB0aW1lIHRoYW4gc2ltcGx5IHJ1bm5pbmcgdGhlIGNydWRlIHZlcnNpb24gbGlua2VkIGhlcmUuCgojIEdlbm9taWMgcHJlZGljdGlvbgoKSWYgeW91IHdhbnQgdG8gc2VlIHdoYXQgYSBjcnVkZSBhcHByb2FjaCB0byBnZW5vbWljIHByZWRpY3Rpb24gbG9va3MgbGlrZSB3aXRoIHRoZSBneXBzeW1vdGggZGF0YSwgcmVmZXIgdG8gW2dlbm9taWMtcHJlZGljdGlvbi5uYi5odG1sXShnZW5vbWljLXByZWRpY3Rpb24ubmIuaHRtbCkgKFtSbWQgdmVyc2lvbl0oZ2Vub21pYy1wcmVkaWN0aW9uLlJtZCkpCg==