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.
LS0tCnRpdGxlOiAiR1dBUyBvbiBneXBzeSBtb3RocyByZXN1bHRzIgpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sKLS0tCgojIFJlc3VsdHMKCllvdSdsbCBmaW5kIHRoZSByZXN1bHRzIGluIHRocmVlIHNlcGFyYXRlIEV4Y2VsIHdvcmtib29rczoKCi0gW0dXQVMtTWFzcy1yZXN1bHRzLnhsc3hdKGh0dHA6Ly9kYXJ3aW4uZWViLnVjb25uLmVkdS9lZWIzNDgtcmVzb3VyY2VzL0dXQVMtTWFzcy1yZXN1bHRzLnhsc3gpCi0gW0dXQVMtUEQtcmVzdWx0cy54bHN4XShodHRwOi8vZGFyd2luLmVlYi51Y29ubi5lZHUvZWViMzQ4LXJlc291cmNlcy9HV0FTLVBELXJlc3VsdHMueGxzeCkKLSBbR1dBUy1URFQtcmVzdWx0cy54bHN4XShodHRwOi8vZGFyd2luLmVlYi51Y29ubi5lZHUvZWViMzQ4LXJlc291cmNlcy9HV0FTLVREVC1yZXN1bHRzLnhsc3gpCgojIFJ1bm5pbmcgdGhlIGFuYWx5c2VzCgpJZiB5b3UncmUgaW50ZXJlc3RlZCBpbiBob3cgSSByYW4gdGhlIGFuYWx5c2VzLCB5b3UnbGwgZmluZCB0aGUgY29kZSBpbiBbR1dBUy5yXShodHRwOi8vZGFyd2luLmVlYi51Y29ubi5lZHUvZWViMzQ4LXJlc291cmNlcy9HV0FTLnIpLiBKdXN0IHB1dCB0aGUgY29kZSBpbiB0aGUgc2FtZSBkaXJlY3RvcnkgYXMgYGd5cHN5bW90aC5jc3ZgIGFuZCB0eXBlIGBzb3VyY2UoIkdXQVMuciIpYCBhdCB0aGUgY29tbWFuZCBwcm9tcHQgaW4gYW4gYFJgIGNvbnNvbGUuIEl0IHdpbGwgY2h1cm4gZm9yIGEgZGF5IG9yIGEgZGF5IGFuZCBhIGhhbGYgYW5kIHB1dCB0aHJlZSBDU1YgZmlsZXMgaW4geW91ciB3b3JraW5nIGRpcmVjdG9yeS4gVGhlIEV4Y2VsIHdvcmtib29rcyBjb21lIGZyb20gb3BlbmluZyB0aGUgQ1NWIGZpbGVzIGFuZCBhZGRpbmcgYSBsaXR0bGUgZm9ybWF0dGluZy4KCklmIEknZCBoYWQgYSBiaXQgbW9yZSB0aW1lLCBJIGNvdWxkIGhhdmUgd3JpdHRlbiBhIHZlcnNpb24gb2YgYEdXQVMucmAgdGhhdCByYW4gdGhlIGFuYWx5c2VzIG11Y2ggbW9yZSBxdWlja2x5LCBidXQgZ2V0dGluZyB0aGF0IHJpZ2h0IHdvdWxkIGhhdmUgdGFrZW4gbW9yZSB0aW1lIHRoYW4gc2ltcGx5IHJ1bm5pbmcgdGhlIGNydWRlIHZlcnNpb24gbGlua2VkIGhlcmUuCgojIEdlbm9taWMgcHJlZGljdGlvbgoKSWYgeW91IHdhbnQgdG8gc2VlIHdoYXQgYSBjcnVkZSBhcHByb2FjaCB0byBnZW5vbWljIHByZWRpY3Rpb24gbG9va3MgbGlrZSB3aXRoIHRoZSBneXBzeW1vdGggZGF0YSwgcmVmZXIgdG8gW2dlbm9taWMtcHJlZGljdGlvbi5uYi5odG1sXShnZW5vbWljLXByZWRpY3Rpb24ubmIuaHRtbCkgKFtSbWQgdmVyc2lvbl0oZ2Vub21pYy1wcmVkaWN0aW9uLlJtZCkpCg==