Wang minjie
@perlatex.bsky.social
R, Tidyverse, Tidymodels and Stan
Reposted by Wang minjie
![df = data.frame(x = -1:10, y = c(-3,2,1,2,5,3,6,7,3,6,10,11))
smooth_spline_x = \(formula, data) {
X = model.matrix(formula[-2], data)
# this is dumb but good enough for a hack
X %*% rep(1, ncol(X))
}
smooth_spline = \(formula, data, ..., weights = NULL) {
x = smooth_spline_x(formula, data)
y = eval(formula[[2]], data, environment(formula))
structure(
splinefun(x, y, ...),
class = c("smooth_spline", "function"),
formula = formula
)
}
predict.smooth_spline = \(object, newdata, ..., weights = NULL, se.fit = FALSE) {
fit = object(smooth_spline_x(attr(object, "formula"), newdata))
if (se.fit) {
list(fit = data.frame(fit = fit, lwr = 0, upr = 0), se.fit = NA_real_)
} else {
fit
}
}
df |>
ggplot(aes(x, y)) +
stat_smooth(formula = y ~ x, method = smooth_spline) +
geom_point()](https://cdn.bsky.app/img/feed_thumbnail/plain/did:plc:6fs5x6cofzyz4si2xhqiqj37/bafkreiacuzr4k36gdw3hfpomiayv5a6f3tegppea5czle5yxqaap5aqelm@jpeg)
Reposted by Wang minjie
Lightnews
Add to your homescreen for quick access.
Add to Home Screen
Works just like the mobile app.
Install this app on your iPhone: tap the share button and then Add to Home Screen.
1
Tap the Share button
in Safari's toolbar
2
Select "Add to Home Screen"
from the menu
3
Tap "Add"
The app will appear on your home screen
