zx8754
@zx8754.bsky.social
Analysis - Data - Analysis - Data - Volleyball - Father - #rstats #STEM #ProstateCancer
November 4, 2025 at 7:34 AM
"my default assumption is nothing is going to work" - @clauswilke.com
October 23, 2025 at 2:04 PM
"my default assumption is nothing is going to work" - @clauswilke.com
We need to have a conversation about random seeds. Don't use 42.
blog.genesmindsmachines.com/p/if-your-ra...
blog.genesmindsmachines.com/p/if-your-ra...
If your random seed is 42 I will come to your office and set your computer on fire🔥
Figuratively. More likely you'll get a stern talking to.
blog.genesmindsmachines.com
October 23, 2025 at 8:33 AM
Reposted by zx8754
Men with BRCA1 and BRCA2 mutations should get annual #ProstateCancer screening.
The IMPACT trial finds BRCA1 carriers are 3 times more likely to have aggressive prostate cancers than non-carriers.
Results from the international trial shared at @myesmo.bsky.social
www.icr.ac.uk/about-us/icr...
The IMPACT trial finds BRCA1 carriers are 3 times more likely to have aggressive prostate cancers than non-carriers.
Results from the international trial shared at @myesmo.bsky.social
www.icr.ac.uk/about-us/icr...
October 17, 2025 at 8:05 AM
Men with BRCA1 and BRCA2 mutations should get annual #ProstateCancer screening.
The IMPACT trial finds BRCA1 carriers are 3 times more likely to have aggressive prostate cancers than non-carriers.
Results from the international trial shared at @myesmo.bsky.social
www.icr.ac.uk/about-us/icr...
The IMPACT trial finds BRCA1 carriers are 3 times more likely to have aggressive prostate cancers than non-carriers.
Results from the international trial shared at @myesmo.bsky.social
www.icr.ac.uk/about-us/icr...
![Barplots of population change ratios in "Most populous countries in Europe" 2000 vs 2025.
R code for the plot:
library(data.table)
library(ggplot2)
library(scales)
d <- fread("
Rank,Country,y2000,y2025
1,Russia,146000000,144000000
2,Germany,82800000,84080000
3,Britain,59510000,69550000
4,France,59330000,66660000
5,Italy,57640000,59150000
6,Ukraine,49430000,38980000
7,Spain,40000000,49320000
8,Poland,38650000,38140000")
d[, rate := y2025/y2000 ]
d[, Country := factor(Country, levels = d[ order(rate), Country]) ]
d[, lab := paste0(round(rate, 4), "\n",
label_number(accuracy = 0.1, scale_cut = cut_short_scale())(y2025 - y2000)) ]
d[, col := as.factor(sign(y2025-y2000)) ]
ggplot(d, aes(Country, rate, label = lab, colour = col)) +
geom_col(fill = "white", show.legend = FALSE, linewidth = 1.1) +
geom_hline(yintercept = 1, colour = "blue", linetype = "dashed") +
geom_text(vjust = 1.4, colour = "black") +
scale_x_discrete(name = NULL) +
scale_y_continuous(name = NULL, breaks = c(0, 1), limits = c(0, 1.3), expand = c(0, 0)) +
scale_color_manual(values = c("-1"="red", "1" = "darkgreen")) +
labs(title = "Most populous countries in Europe",
subtitle = "2000 vs 2025",
caption = "Source: @GEOMAPAS.GR | https://bsky.app/profile/simongerman600.bsky.social/post/3m2nz34kzyi27") +
theme_minimal() +
theme(panel.grid = element_blank())
ggsave("tmp.plot.jpeg", width = 8, height = 6)](https://cdn.bsky.app/img/feed_thumbnail/plain/did:plc:uudl4famc7jqkejakihsxg5a/bafkreifhd5umrvsvgajnglktemay2ghenu7f7icdlbsmfx2nh7gi3d4zrq@jpeg)
#rstats pipe help, where is this from? I am guessing it is "modify and re-assign".
%<+%
Google didn't help, Gemini says it is invalid R code or it is a regex. ChatGPT says it is from ggplot, searching GitHub didn't help. But found this function ggplot2::`%+replace%`(), are they related?
%<+%
Google didn't help, Gemini says it is invalid R code or it is a regex. ChatGPT says it is from ggplot, searching GitHub didn't help. But found this function ggplot2::`%+replace%`(), are they related?
September 1, 2025 at 7:48 AM
#rstats pipe help, where is this from? I am guessing it is "modify and re-assign".
%<+%
Google didn't help, Gemini says it is invalid R code or it is a regex. ChatGPT says it is from ggplot, searching GitHub didn't help. But found this function ggplot2::`%+replace%`(), are they related?
%<+%
Google didn't help, Gemini says it is invalid R code or it is a regex. ChatGPT says it is from ggplot, searching GitHub didn't help. But found this function ggplot2::`%+replace%`(), are they related?
As base::rbind also takes more than two dataframes as input, I expected at least the same from dplyr::bind_rows. And bind_rows handles rownames, factor levels, and mismatched column names better. #rstats
![#list of dfs
set.seed(1); listOfDfs <- lapply(1:3, \(i){ mtcars[ sample(1:nrow(mtcars), 4), ]})
#input as list
do.call(rbind, listOfDfs)
bind_rows(listOfDfs)
#input as dataframe
rbind(listOfDfs[[ 1 ]], listOfDfs[[ 2 ]], listOfDfs[[ 3 ]])
bind_rows(listOfDfs[[ 1 ]], listOfDfs[[ 2 ]], listOfDfs[[ 3 ]])](https://cdn.bsky.app/img/feed_thumbnail/plain/did:plc:uudl4famc7jqkejakihsxg5a/bafkreifzj3fqdrnu5y3fwibe3kehjmqhew6xtvzmj6hvi43fkckwuqeoeu@jpeg)
July 21, 2025 at 7:37 AM
As base::rbind also takes more than two dataframes as input, I expected at least the same from dplyr::bind_rows. And bind_rows handles rownames, factor levels, and mismatched column names better. #rstats
![set.seed(2025)
d1 <- data.frame(expand.grid(1:50, 1:50))
d1$hot <- factor("no", levels = c("no", "yes"))
d1[ sample(seq(nrow(d1)), 1), "hot" ] <- "yes"
d1$time <- factor("before", levels = c("before", "today"))
d2 <- data.frame(expand.grid(1:50, 1:50))
d2$hot <- factor("no", levels = c("no", "yes"))
d2[ sample(seq(nrow(d2)), 100), "hot" ] <- "yes"
d2$time <- factor("today", levels = levels(d1$time))
library(ggplot2)
ggplot(rbind(d2, d1), aes(Var1, Var2, colour = hot)) +
geom_point() +
facet_wrap(vars(time)) +
scale_colour_manual(values = c("no" = "grey90", "yes" = "red"), guide = "none") +
theme_void()](https://cdn.bsky.app/img/feed_thumbnail/plain/did:plc:uudl4famc7jqkejakihsxg5a/bafkreiciukvjwpu4httb775e6vp5mpunburjftcb5erxkqxc54qpwsrqt4@jpeg)
May 27, 2025 at 8:04 AM
![df <- read.table(text = "stu_id race1 race2 race3 race4 race5 race6
100 0 0 0 0 0 1
101 1 0 0 0 0 0
102 0 1 0 0 0 1", header = TRUE)
f <- function(...){
ix <- which(c(...) == 1)
if(length(ix) == 1) ix else 7L
}
#base
apply(df[, -1], 1, f)
# [1] 6 1 7
#to create new column
#df$race <- apply(df[, -1], 1, f)
#purrr
pmap_int(df[, -1], f)
# [1] 6 1 7
#to create new column
# df <- df %>%
# mutate(race = pmap_int(select(., -1), f))](https://cdn.bsky.app/img/feed_thumbnail/plain/did:plc:uudl4famc7jqkejakihsxg5a/bafkreibsidfwcbtnmx2ann6sd6qur7xz2tcjyidsefcnamew2zldnvhsxm@jpeg)
Movie you've watched more than six times using GIFs. ("Hard mode" no Star Wars, Star Trek or LOTR.)
April 26, 2025 at 10:18 PM
Movie you've watched more than six times using GIFs. ("Hard mode" no Star Wars, Star Trek or LOTR.)
![# tidy: count 1s on subset of columns
# base: try this
colSums(mtcars[, c("vs", "am", "gear", "carb") ] == 1)
# tidy: but I don't want quotes...
# base: ok, cols by index
colSums(mtcars[, 8:11 ] == 1)
# tidy: but I want to see column names!
# base: ok, let's use expressions
colSums(mtcars[, as.character(substitute(list(vs, am, gear, carb)))[ -1 ] ] == 1)
# tidy: oh no, too many nested brackets, my head hurts!
# base: wrap into a fucntion
count_ones_base <- function(data, ...) {
colSums(mtcars[, as.character(substitute(list(...)))[ -1 ] ] == 1)
}
count_ones_base(mtcars, vs, am, gear, carb)
# tidy: but I like pipes, weird function names like quos,
# and to be loud !! about it !!!
# base: ok.
library(dplyr)
library(tidyr)](https://cdn.bsky.app/img/feed_thumbnail/plain/did:plc:uudl4famc7jqkejakihsxg5a/bafkreicdna7pej7k4lgevbym5dhmx2nbletm6u7yjqxujur3dhmrzdk574@jpeg)
April 25, 2025 at 8:38 AM
Please stop sending me your datasets.
YouTube video by Darren Dahly
youtu.be
April 17, 2025 at 10:00 AM
#rstats #ggplot version of this interesting
approach
data and codes are at gist.github.com/zx8754/85d0c...
approach
data and codes are at gist.github.com/zx8754/85d0c...
April 8, 2025 at 9:38 AM
#rstats #ggplot version of this interesting
approach
data and codes are at gist.github.com/zx8754/85d0c...
approach
data and codes are at gist.github.com/zx8754/85d0c...