Debugging my R script

@ramencult Very cool @ramencult excited to read it!

RT @ramencult: New preprint! We find a novel signal of positive selection in georeferenced genotype data, then use this signal to identify…

@vsbuffalo I just try and lean on the language itself very little. Handle anything even a little complex in python (more testable anyways). Use parallel on high cpu instances vs scatter where possible (usually cheaper).

@vsbuffalo We compile to run on DNAnexus with dxcompiler. Miniwdl is handy for small local tests.

@vsbuffalo It might be Stockholm syndrome but I now find WDL pretty good for complex analyses.

My feed is an even mix of anti-immigrant propaganda and recycled TikTok slop. RIP Twitter. I’m out.

Celebrating Hanukkah by eating 8 days worth of latkes in 20 minutes

RT @danrothenberg: ECMWF delivers an early Christmas present with what could arguably be the most important paper in numerical weather pred…

RT @fchollet: Today OpenAI announced o3, its next-gen reasoning model. We've worked with OpenAI to test it on ARC-AGI, and we believe it re…

@srikosuri What a bummer. I really didn’t appreciate how much schools and even daycares are still set up for one nonworking parent until I had kids.

@AlexTISYoung @vsbuffalo @jgschraiber labs with foundation $$ sometimes have these jobs too, but pay is still low and it’s a precarious gig (one job per city, temp funds, limited growth). To balance the downsides you could offer intellectual freedom and tenure, but that’s… a faculty job lol.

RT @spence_jeffrey_: What do GWAS and rare variant burden tests discover, and why? Do these studies find the most IMPORTANT genes? If not,…

@pastramimachine Yep Christmas morning tradition

RT @mbeisen: I don't think @DrJBhattacharya got everything right during the pandemic by any means, but if we're enumerating things incumben…

@pastramimachine Packed up the jnco cargo shorts for faculty life

@byersblake ADRB-1’s candidate gene story doesn’t hold up in large scale studies.

@pvtodorov @vsbuffalo Today I asked Claude for a command to filter reads with more than 15bp soft-clipping and it gave me a nice awk command filtering for >=10bp soft-clipping. Lots of little stuff like that will slip through I think.

@greebletheory @vsbuffalo Still learn bash/awk/sed/grep enough to do stuff. The fun part is coming up with good models linking bio/chem/stats and I think (hope?) that is where the human value added will come from. Also ggplot for exploratory analysis.