博文

目前显示的是 五月, 2025的博文

inlabru 时空分析

library ( INLA ) # 用于之后的建模 library ( MASS ) # 用于真实的多元正态分布模拟 library ( tidyverse ) library ( inlabru ) library ( fdmr ) library ( sf ) library ( inlatools )   covid19_data <- fdmr :: load_tutorial_data ( dataset = "priors" , filename = "covid19_dataBris.rds" ) covid19_sf <- covid19_data |> st_as_sf ( coords = c ( "LONG" , "LAT" )) coordinates <- st_coordinates ( covid19_sf ) locs <- unique ( coordinates ) # 计算初始参数 initial_range <- diff ( range ( locs [, 1 ])) / 3 max_edge <- initial_range / 2 # 2. 构建改进的网格,mesh从地图构建,此处不对。 mesh <- fm_mesh_2d (   loc = locs ,   max.edge = c ( 1 , 2 ) * max_edge ,   offset = c ( initial_range , initial_range ),   cutoff = max_edge / 5 ) # 3. 定义SPDE模型 spde <- inla.spde2.pcmatern (   mesh = mesh ,   prior.range = c ( initial_range , 0.5 ),   prior.sigma = c ( 1 , 0.01 ) ) # 4. 设置时间趋势先验, rw1_prior <- list ( theta = list ( prior = "pc.prec" , param = c ( 1 , ...

计算工作日

 library(tidyverse) # 定义节假日和补班日 holidays_2024 <- as.Date(c(   "2024-01-01",  # 元旦   "2024-02-10", "2024-02-11", "2024-02-12", "2024-02-13", "2024-02-14", "2024-02-15", "2024-02-16",  # 春节   "2024-04-04", "2024-04-05", "2024-04-06",  # 清明节   "2024-05-01", "2024-05-02", "2024-05-03",  # 劳动节   "2024-06-08", "2024-06-09", "2024-06-10",  # 端午节   "2024-09-15", "2024-09-16", "2024-09-17",  # 中秋节   "2024-10-01", "2024-10-02", "2024-10-03", "2024-10-04", "2024-10-05", "2024-10-06", "2024-10-07"  # 国庆节 )) workdays_2024 <- as.Date(c(   "2024-02-04", "2024-02-18",  # 春节补班   "2024-04-07",  # 清明补班   "2024-04-28", "2024-05-11",  # 劳动节补班   "2024-06-02",  # 端午补班   "20...