Reposted by Eric Book
Working with spatial data in R? The terra cheat sheet is a quick reference for raster/vector operations, projections, plotting, and more.
Check it at rspatial.org/_static/terra_cheatshe…
#rstats #GISchat #rspatial
Eric Book
@erc-bk.fosstodon.org.ap.brid.gy
11 followers 8 following 89 posts
Pathologically curious about most things, but this space is mostly focused on data science and engineering. #tfr 🌉 bridged from ⁂ fosstodon.org/@erc_bk, follow @ap.brid.gy to interact
![time_slices <- 300 # 300 days
block_length <- 14 # 2 weeks
n_samples <- 200 # 200 resamples
block_idxs <- make_block_indices(
n_time = dim(r5to10)[2],
time_slices = time_slices,
block_length = block_length,
n_samples = n_samples
)
plan(future.mirai::mirai_multisession, workers = 10L)
tib_vario_params <- purrr::map(
asplit(block_idxs, 1),
\(x) {
estim_params_vgm(x, dat = r5to10)
}
) |>
futurize() |>
list_rbind()
mirai::daemons(0)](https://cdn.bsky.app/img/feed_thumbnail/plain/did:plc:k7ohhh6vlfkiw2f7a4mql3fh/bafkreif3qtaystniwxrmb27drggzdrqzkoa2ujbvp624d3balogflyrlim)
![Running dbreg::dbreg() on a full year of NYC taxi data... and it takes less than 2 seconds.
dbreg(
tip_amount ~ fare_amount + passenger_count | month + vendor_name,
path = "read_parquet('nyc-taxi/**/*.parquet')", ## path to hive-partitioned dataset
vcov = "hc1"
)
#> [dbreg] Estimating compression ratio...
#> [dbreg] Data has 178,544,324 rows and 24 unique FE groups.
#> [dbreg] Using strategy: compress
#> [dbreg] Executing compress strategy SQL
#>
#> Compressed OLS estimation, Dep. Var.: tip_amount
#> Observations.: 178,544,324 (original) | 70,782 (compressed)
#> Standard Errors: Heteroskedasticity-robust
#> Estimate Std. Error t value Pr(>|t|)
#> fare_amount 0.106744 0.000068 1564.742 < 2.2e-16 ***
#> passenger_count -0.029086 0.000106 -273.866 < 2.2e-16 ***
#> ---
#> Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1](https://cdn.bsky.app/img/feed_thumbnail/plain/did:plc:qww63z2t5rrhniusz5x7gk23/bafkreiajmnp6jiaycebutwid2lxx57hrz3zc3qp2p5ijnoqiyqag34xpe4)