Rafael Irizarry
@rafalab.bsky.social
2.9K followers 12 following 20 posts
Applied statistician. I tweet data-driven observations, data science educational materials, academic research updates, and the occasional joke.
Posts Media Videos Starter Packs
rafalab.bsky.social
We agree. I meant function in its most basic definition: every p-dimensional x in your dataset is mapped to a unique 2-dimensional f(x). I did not claim f is defined for all p-dimensional xs.

My point is that this f is difficult or impossible to describe. In contrast, we can write it down for PCA.
rafalab.bsky.social
I like it. Thank for the tip.
rafalab.bsky.social
I agree. I don't use them at all with genomics data, especially sparse noisy scRNA-Seq data.

It does appear to perform impressively well with high signal-to-noise ratio datasets, such as MNIST.
rafalab.bsky.social
It depends on the specific biological insight you want to highlight or communicate.
rafalab.bsky.social
The inertia here is strong. I have not been unable to convince collaborators to not use them on papers I am a co-author on...
I'll keep trying though.

At some point I might give up as I did with pie charts: simplystatistics.org/posts/2012-1...
Simply Statistics: I give up, I am embracing pie charts
simplystatistics.org
rafalab.bsky.social
As made clear in the blogpost I am not against UMAP either. But when I see a plot in a paper, I want to understand what I am being shown and why.

Other than to show different cell types have different expression patterns, which I already know, or to decorate,
why use UMAP to display in 2D?
rafalab.bsky.social
What am I supposed to learn from that plot? Cell types have different expression patterns. Those are markers for different cell types. So, this just confirms something obvious. Do all those non-linear shapes and tiny clusters represent anything biological?
rafalab.bsky.social
To be clear, as the post explains, UMAP can be useful for exploratory data analysis. My concern is their inclusion in papers as if they were results. What exactly is the reader supposed to learn? And how often are we misdirected by false clusters or artifactual shapes?
rafalab.bsky.social
Can you explain what the axes represent?

As mentioned in the post, UMAP can be useful for exploring data. But why are plots included in papers? What is the reader supposed to get out of them? The 2D distance between points can't be interpreted.

It seems the only reason is because they are pretty.
rafalab.bsky.social
There are plenty of alternatives. They don’t produce flashy art work but they do provide scientific insights.

If journals want art work no need to pretend we are analyzing data. Just paint pretty pictures.
rafalab.bsky.social
This is unfortunately true. I would say the main reasons are that the subject is hard and deep understanding is not incentivized enough.

But note understanding UMAP is much harder than understanding pvalues.
rafalab.bsky.social
🚨 We're hiring! 🚨

Join us at @HarvardBiostats and @dfcidatascience
We are filling THREE positions:

1️⃣ Tenure-track (any rank) in AI/ML
2️⃣ Assistant Professor in Single Cell Genomics
3️⃣ Lecturer & Director Training/Education
Need to apply separately for each role. Links 👇
rafalab.bsky.social
New method for identifying spatially variable genes in spatial transcriptomics.

We use spectral graph theory to extract 1D curves related to tissue morphology, creating a more relevant coordinate system. Improves sensitivity.

with @PhillipNicol et al.

www.biorxiv.org/content/10.1...
rafalab.bsky.social
Question for #rstats educators: Do you use a standard format for exercises?

I have 100s of exercises along with solutions in quarto documents (see image). But they are not formatted in a systematic way. I would like to easily create docs with and w/o solutions and export to LMS.
rafalab.bsky.social
Hoy inicia nuestro Curso en español de Ciencia de Datos: Fundamentos de R en HarvardX! Aprenderás desde los básicos hasta data wrangling, conceptos de programación e importación de datos. La versión gratis y en línea te da acceso a los videos, textos y ejercicios.
t.co/TGID4IwTSQ
rafalab.bsky.social
Starting to work on it 😅
rafalab.bsky.social
Data visualization is used in two ways in data analysis:

1) Communicating data insights.
2) Guiding analysis and spotting data issues.

While modern courses and books focus on the first, it's the second - exploratory data analysis - that is more important and harder to teach.
rafalab.bsky.social
Introducing the Data Science Postdoctoral Fellows Program at Harvard/DFCI!

- Join a research group in our department
- Co-mentoring opportunities with 2+ faculty
- Collaborate with DFCI investigators beyond our department
- Salary starts at $75K

Apply here: ds.dfci.harvard.edu/postdocs/