and an instance wide view for all of your items—also, 👀 finally got around to supporting @jupyter.org notebooks. #databs #rstats
anyways! a detailed post soon enough! super stoked 🙃🐇
Josiah
@josiah.rs
2.7K followers 509 following 764 posts
putting R in production @ricochet.rs doing spatial stats @ Esri blog: josiah.rs
![[services]
web = {
image = "nginx:latest",
environment = { DB_HOST = "postgres", DB_PORT = 5432 },
resources = {
limits = { memory = "512M", cpu = "0.5" }
}
}](https://cdn.bsky.app/img/feed_thumbnail/plain/did:plc:63atwz67mlrdquqnasarutmr/bafkreifmpl3o6ejyjbbjo4cbodktptz3bnszvb4o5w4n5gijtdtsytegou)
![use anyhow::bail;
use arrow_extendr::FromArrowRobj;
use extendr_api::prelude::*;
use geo_traits::{CoordTrait, PointTrait};
use geoarrow::array::{GeoArrowArrayAccessor, PointArray};
#[extendr]
fn geoarrow_centroid(x: Robj) -> anyhow::Result<Doubles> {
let res = PointArray::from_arrow_robj(&x)?;
let (sum_x, sum_y, n) = res.iter_values().filter_map(|xi| xi.ok()).fold(
(0.0_f64, 0.0_f64, 0usize),
|(sx, sy, n), v| {
let (x, y) = v.coord().unwrap().x_y();
(sx + x, sy + y, n + 1)
},
);
if n == 0 {
bail!("cannot compute centroid of empty array");
}
let res = [sum_x / n as f64, sum_y / n as f64]
.map(Rfloat::from)
.into_iter()
.collect::<Doubles>();
Ok(res)
}](https://cdn.bsky.app/img/feed_thumbnail/plain/did:plc:63atwz67mlrdquqnasarutmr/bafkreifvhwk7acara7a64jwglprbjgszw5vd7dmfox55qpyhhe4c5i3sry)
![library(wk)
library(geoarrow)
devtools::load_all()
#> ℹ Loading extendrtest
pnts <- wk::xy(rnorm(100, -180, 180), runif(100, -90, 90))
x <- geoarrow::as_geoarrow_array(pnts)
geoarrow_centroid(x)
#> [1] -182.711392 -2.806439
sf::st_as_sfc(pnts) |>
sf::st_combine() |>
sf::st_centroid()
#> Geometry set for 1 feature
#> Geometry type: POINT
#> Dimension: XY
#> Bounding box: xmin: -182.7114 ymin: -2.806439 xmax: -182.7114 ymax: -2.806439
#> CRS: NA
#> POINT (-182.7114 -2.806439)](https://cdn.bsky.app/img/feed_thumbnail/plain/did:plc:63atwz67mlrdquqnasarutmr/bafkreifb7wrdfujkxs356k5dj6peukdmtj2hekllcpbheckr4uhihc4r5a)
![ggplot() +
geom_sf(data = nodes[c(i, nb[[i]])]) +
geom_sf(data = nodes[i, ])](https://cdn.bsky.app/img/feed_thumbnail/plain/did:plc:63atwz67mlrdquqnasarutmr/bafkreifbvzushklwvmahtarwubkusbrgcng37bjbnqg2vdyume3wi2p62a)
![Screenshot of a failed github action workflow "Check for non-API calls."
Text output:
Run RESPONSE=$(curl -sf "https://dev.vm.ubuntu.ricochet.rs/app/non-api-checks/check?pkg=$PKG_NAME")
{
"last_updated": "2026-04-10 17:52:44 CEST",
"has_non_api_calls": true,
"symbols": [
"R_UnboundValue"
]
}
Error: non-API calls detected in b64
Error: Process completed with exit code 1.](https://cdn.bsky.app/img/feed_thumbnail/plain/did:plc:63atwz67mlrdquqnasarutmr/bafkreig2ldfi2cyrirfackg3uxrnvb5a76i7i5dpuapy3iphoxa5d2ouya)
![library(rlang)
macro_fun <- function(name) {
f <- new_function(alist(x=), quote(x^2))
nm <- sym(name)
env_bind(caller_env(), !!nm := f)
}
macro_fun("square2")
square2(5)
#> [1] 25](https://cdn.bsky.app/img/feed_thumbnail/plain/did:plc:63atwz67mlrdquqnasarutmr/bafkreigbivoy7rtc24uoetrt2ewvslfuzu3azgb7lrk6nc5unh4k6arugu)
