Robert Cope

(Kangaroos on campus is normal but not quite like this)

A path with bushes on either side. A kangaroo is barely visible amongst the bushes.

🦘 certainly well camouflaged

Yeah I think if you had an institutional account and ran your whole class through it that'd be quite effective. I think I'd have some trouble getting IT approval, and most students manage to install Rstudio on their own computers so it'd be a lot of logistics just for the edge cases.

We send the cursed ones to use posit.cloud free accounts but I wish there was a nicer solution

Kanji Explorer available at https://jcarroll.shinyapps.io/kanjiexplorer/

I Vibe Coded an R Package
and it ... actually works?!?

https://jcarroll.com.au/2025/09/13/i-vibe-coded-an-r-package/

I wanted an #rstats package to exist so I spent $20 and got Claude Code to build it for me.

Plus a shiny app using it: https://jcarroll.shinyapps.io/kanjiexplorer/

my lil reproducibility talk from today / I really wanted to instill the PhD students some simple first practices and ways to step up your game from there github.com/tjmahr/2025-...

Wow. The discussions in the Annals of Statistics were... rough, in the 80s 😅

That being said, interesting paper and discussion. (Thanks to Yannis Yatracos who sent me this reference!)

Scientists uncover the secret lives of dugongs in Moreton Bay The world's largest health study of dugongs in Moreton Bay provides critical clues for helping conserve the enigmatic marine mammal, especially in other locations where numbers are falling.

Dugongs!! This is a great article. I was lucky to be part of this team in the distant past so it's great they are still doing cool science

Job alert: 2.75 -3 yr postdoc in Mathematical ecology/epidemiology or Decision Sci at the University of Queensland on the value of information for decisions with @hugepossum.bsky.social @katehelmstedt.bsky.social Katriona Shea & myself #MathSky 🌍 uq.wd3.myworkdayjobs.com/uqcareers/jo... please share

Happy to discover that R.A. Fisher already relied on what one may refer to as "resolution by cope"

Graphical abstract for "Vocal communication is seasonal in social groups of wild, free-living house mice." The abstract has, from top to bottom, a title, four middle image panels, and two bottom text panels. Image title: "Vocal communication in social groups of wild-free living house mice" Middle image panels from left to right: (1) An aerial snap shot of the region where the study site is located, an agricultural landscape in rural Switzerland. (2) An image of the study site, a small barn in the forest inhabited by mice. (3) An image of a radio frequency identification (RFID) box used to track mouse social interactions. A mouse is entering the box from the left while another sits outside. (4) A spectrogram showing example vocalizations - one low frequency squeak and one ultrasonic call - recorded from an RFID box. Bottom panels: Left: Data Collection  - 10 years of RFID-based tracking data (from 6,946 mice) - 15 months of acoustic monitoring (totaling 6,594 hours) - Machine learning for vocal detection and labeling (CNN) Right: Key Findings - Vocalization is seasonal (most in spring and summer) - Vocalization is associated with the presence of pups - Vocalization is correlated with social group dynamics

Very happy to share the latest from my postdoc‬!

10 yrs of mouse social networks + 1.25 yrs of acoustic data ➡️ insight into vocalization & sociality in a wild population of your favorite lab model 🐁

paper: bit.ly/4n93yyD
data: bit.ly/4lfFBEk
code: bit.ly/4kNnMwx

#bioacoustics #neuroskyence

1/8

Little Ume, our endangered tapir calf, has loved meeting so many of you over the past two weekends! 🩷 This growing baby needs lots of rest, so she’s been catching plenty of zzz’s behind the scenes and in her habitat. 😴

Please share. 2 tenure-track T&R math/stat jobs at the University of Queensland, Australia. Amazing weather, gorgeous campus & lifestyle. 17% employer retirement contribution, lots of leave. Strong global reputation. Supports diversity @stevenstrogatz.com uq.wd3.myworkdayjobs.com/en-US/uqcare...

Sorry, Modsim!

Should I go to the Australian Statistical Conference or Modism?!? They both look cool but are in direct competition...

YouTube video by JuliaHub Converting R data-wrangling code to julia

🔥🔥🔥

Check out this Tidier.jl (and general #julialang data) tutorial from @kevinbonham.com @juliahub.bsky.social!

#rstats #tidyverse

youtu.be/bw-N1lrDeHI?...

Sum-score effect sizes for multilevel Bayesian cumulative probit models | A. Solomon Kurz This is a follow-up to my earlier post, *Notes on the Bayesian cumulative probit*. This time, the topic we're addressing is: *After you fit a full multilevel Bayesian cumulative probit model of severa...

If I'm following correctly, you can do this as a Bayesian with brms like so: solomonkurz.netlify.app/blog/2022-07...

Lecturer: Mathematical Sciences (Academic Level B) Job in Mawson Lakes, Adelaide SA - SEEK Contribute to research in Mathematical Sciences focusing on Applied Mathematics including dynamics, modelling & computation

2-year lecturer position in our Applied Maths / Modelling / Computation group at Uni SA / Adelaide Uni - www.seek.com.au/job/83327354

Causal Inference With Observational Data and Unobserved Confounding Variables As ecology tackles progressively larger problems, we are moving beyond the scales at which randomised controlled experiments are feasible. Using observational data for causal inference raises the pro...

OK, #causalinference #ecology folk - my paper with the even excellent @lauradee.bsky.social on causal inference with observational data is now out in Ecology Letters! 🌍🧪🌊

onlinelibrary.wiley.com/doi/full/10....

Screenshot of the linked Quarto website, with input checkboxes to change different conditions for a regression model that predicts economic performance based on US political party, with a reported p-value

I’ve long used FiveThirtyEight’s interactive “Hack Your Way To Scientific Glory” to illustrate the idea of p-hacking when I teach statistics. But ABC/Disney killed the site earlier this month :(

So I made my own with #rstats and Observable and #QuartoPub ! stats.andrewheiss.com/hack-your-way/

I'm just saying this is syntatically correct #rstats code

The latest version of the {mapgl} #rstats package is now released, supporting scrollytelling-style story maps!

Build richly interactive map stories from your R data pipelines using @Mapbox or @maplibre as your backend.

Learn how in the new vignette: walker-data.com/mapg...

ANU-AAGI team featuring investigators: Eric Stone, Emi Tanaka, Alan Welsh and Francis Hui

Want to use your #statistics, ML/AI or #bioinformatics skills to solve problems in #agriculture?

We are seeking 3 postdocs (up to 3 years) to join at ANU Biological Data Science Institute as part of Analytics for 🇦🇺 Grains industry by 📅 19 Jan

We now also have a (WIP) website!

🌐 anu-aagi.github.io

YouTube video by TheCoatlessProfessor Run R Calculations on Remote Clusters with {future}

Scale R analyses across cluster systems without leaving RStudio! {future} makes it simple to login, schedule & parallelize your code.

🎥 Watch: www.youtube.com/watch?v=G8MC...
💻 Code: github.com/coatless-vid...

#rstats #HPC #datascience

Image shows sample palettes for UChicago, Star Trek and Tron Legacy

The {ggsci} #RStats 📦 “offers a collection of ggplot2 color palettes inspired by scientific journals, data visualization libraries, science fiction movies, and TV shows.” By Nan Xiao.
Use with:
scale_color_palname()
scale_fill_palname()
nanx.me/ggsci/

#RStats #ggplot2 #ggplot

If you're interested in trying out LLMs in #rstats but don't know where to begin, I've added a few two vignettes to elmer: elmer.tidyverse.org/articles/elm... and elmer.tidyverse.org/articles/pro...

The {qpdf} #RStats 📦can transform PDF files with functions like pdf_split(), pdf_combine(), and pdf_rotate_pages().
cran.r-project.org/web/packages...

The {pdftools} #R 📦can read content & metadata in PDFs. docs.ropensci.org/pdftools/
Both by @jeroenooms.bsky.social

GitHub - fmichonneau/2024-latinr-duckdb-arrow: Content of the Tutorial "Working with larger than memory data in R with Arrow and DuckDB" taught on 2024-11-19 Content of the Tutorial "Working with larger than memory data in R with Arrow and DuckDB" taught on 2024-11-19 - fmichonneau/2024-latinr-duckdb-arrow

I recently taught a tutorial on #DuckDB and #ApacheArrow with #Rstats at @latinrconf.bsky.social The slides and content are available on GitHub github.com/fmichonneau/...
The tutorial was recorded and I'll add the link to this thread when it becomes available.

Why we (usually) don't have to worry about multiple comparisons Applied researchers often find themselves making statistical inferences in settings that would seem to require multiple comparisons adjustments. We challenge the Type I error paradigm that underlies t...

Since it is inefficient to attempt to educate every reviewer individually, I am yeeting into your feed this clear paper from Gelman, Hill, and Yajima on how Bayesians can do even better than correcting for multiple comparisons. arxiv.org/abs/0907.2478

You can actually change R’s localization settings to get output in different languages! If you want to see what your computer is currently set to use, run Sys.getLocale(): > Sys.getlocale() ## [1] "en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8" There’s a bunch of output there—the first part (en_US.UTF-8) is the most important and tells you the language code. The code here follows a pattern and has three parts: - A language: en. This is the langauge, and typically uses a two-character abbreviation following the ISO 639 standard - A territory: US. This is the country or region for that language, used mainly to specify the currency. If it’s set to en_US, it’ll use US conventions (like “$” and “color”); if it’s set to en_GB it’ll use British conventions (like “£” and “colour”). It uses a two-character abbreviation following the ISO 3166 standard. - An encoding: UTF-8. This is how the text is actually represented and stored on the computer. This defaults to Unicode (UTF-8) here. You don’t generally need to worry about this.

For macOS and Linux (i.e. Posit Cloud), setting locale details is pretty straightforward and predictable because they both follow this pattern consistently: - en_GB: British English - fr_FR: French in France - fr_CH: French in Switzerland - de_CH: German in Switzerland - de_DE: German in Germany If you run locale -a in your terminal (not in your R console) on macOS or in Posit Cloud, you’ll get a list of all the different locales your computer can use. Here’s what I have on my computer: [1] "af_ZA" "am_ET" "be_BY" "bg_BG" "C" "ca_ES" "cs_CZ" "da_DK" "de_AT" "de_CH" "de_DE" "el_GR" "en_AU" "en_CA" [15] "en_GB" "en_IE" "en_NZ" "en_US" "es_ES" "et_EE" "eu_ES" "fi_FI" "fr_BE" "fr_CA" "fr_CH" "fr_FR" "he_IL" "hi_IN" [29] "hr_HR" "hu_HU" "hy_AM" "is_IS" "it_CH" "it_IT" "ja_JP" "kk_KZ" "ko_KR" "lt_LT" "nl_BE" "nl_NL" "no_NO" "pl_PL" [43] "POSIX" "pt_BR" "pt_PT" "ro_RO" "ru_RU" "sk_SK" "sl_SI" "sr_YU" "sv_SE" "tr_TR" "uk_UA" "zh_CN" "zh_HK" "zh_TW" For whatever reason, Windows doesn’t use this naming convention. It uses dashes or full words instead, like en-US or american or en-CA or canadian. You can see a list here, or google Windows language country strings (that’s actually RStudio’s official recommendation for finding Windows language codes)

Those functions have a locale argument, though, so it’s really easy to switch between languages: > wday(x, label = TRUE, abbr = FALSE, locale = "en_US") ## [1] Friday ## Levels: Sunday < Monday < Tuesday < Wednesday < Thursday < Friday < Saturday > wday(x, label = TRUE, abbr = FALSE, locale = "fr_FR") ## [1] Vendredi ## Levels: Dimanche < Lundi < Mardi < Mercredi < Jeudi < Vendredi < Samedi > wday(x, label = TRUE, abbr = FALSE, locale = "fr_BE") ## [1] Vendredi ## Levels: Dimanche < Lundi < Mardi < Mercredi < Jeudi < Vendredi < Samedi > wday(x, label = TRUE, abbr = FALSE, locale = "it_IT") ## [1] Venerdì ## Levels: Domenica < Lunedì < Martedì < Mercoledì < Giovedì < Venerdì < Sabato > wday(x, label = TRUE, abbr = FALSE, locale = "zh_CN") ## [1] 星期五 ## Levels: 星期日 < 星期一 < 星期二 < 星期三 < 星期四 < 星期五 < 星期六

You can also set the locale for your entire R session like this: > Sys.setlocale(locale = "de_DE") ## [1] "de_DE/de_DE/de_DE/C/de_DE/en_US.UTF-8" Now month() and wday() will use German by default without needing to set the locale argument: > month(x, label = TRUE, abbr = FALSE) ## [1] Juli ## 12 Levels: Januar < Februar < März < April < Mai < Juni < Juli < August < September < Oktober < ... < Dezember > wday(x, label = TRUE, abbr = FALSE) ## [1] Freitag ## Levels: Sonntag < Montag < Dienstag < Mittwoch < Donnerstag < Freitag < Samstag

Lots of the students in my dataviz class speak other languages, so I wrote up a little guide about how to change the locale in #rstats to get localized month/day names with {lubridate} and other R functions

datavizf24.classes.andrewheiss.com/news/2024-11...