Oh, look! A paper about the coolest piece of software I've ever built!
Oh, look! A paper about the coolest piece of software I've ever built!
Thank you for the kind words! I’m a big fan of your recent work on evaluating qpAdm. As I was reading it, I was like “slendr would've made things so much easier for these guys!”. :D
Cool stuff! Been referring to it a lot this year.
More exciting #rstats popgen announcements are on the horizon, so stay tuned!
[I have no idea how threads are supposed to work here. Let's see if I made some embarrassing blunder.]
The new slendr is already on CRAN, so go ahead and run `install.packages("slendr")` to get the latest goodies.
Finally, this update involved massive changes to the slendr internals. I'm actually a bit scared by the scale of those changes! Please report any bugs / issues at github.com/bodkan/slendr.
You can now take a normal slendr demographic model of arbitrary complexity and "inject" a customization SLiM snippet during a `compile_model()` step, overriding the default neutral mode of operation. This gives you infinite SLiM flexibility while keeping everything else in slendr exactly the same!
This thread is getting too long. If you want to know more, check out the WIP vignette with examples capturing a range of WF selection models you might normally use SLiM for: custom mutation types, genomic elements, scheduling specific mutations at given times etc. etc. www.slendr.net/articles/vig...
It was frustrating because slendr is built on top of two simulation engines: 1. msprime, 2. the almighty selection simulator SLiM, which was used exclusively for its handling of space, not for selection (?!). This is no longer true. Since v1.0, slendr models can include any WF-compatible selection.
This brings me to the main feature of the v1.0 release. Having personally come from an aDNA background and primarily studying history, slendr used to only simulate neutral scenarios, frustrating people who wanted to simulate selection happening along (potentially complex) demographic histories.
Simulating ancient individuals at defined points in time, getting various f-statistics, ancestry tracts (chromosome painting), VCF or EIGENSTRAT files, are all supported (powered by tskit!). Very useful for model testing and building intuition as it's really easy to get ground truth from any model.
To my big surprise, slendr became popular even for traditional nonspatial models. In few lines of R and no time, you can get from a model idea to meaningful statistics: a dealbreaker in classrooms and day-to-day exploratory work! See two toy models of f4 / divergence / 2D AFS. All without leaving R!
What is slendr? We originally developed it to simulate "realistic-looking" genomic data from geographically-explicit models. It would take too many posts to detail its awesome spatial features but check out the project website for tutorials and our paper with many cool examples. www.slendr.net
Hex sticker logo of the slendr R package
Let’s celebrate the Bluesky exodus with an exciting software update. 📢
slendr popgen simulation #rstats package has reached version 1.0! This magic number signifies that slendr can now do most of what I had hoped it would do.
TL;DR for those familiar with slendr: you can now simulate selection! 💣