我想对人口数据执行（立方）样条插值，以将年度数据"转换"为季度数据。我知道这样做有相当多的缺陷，但我需要这样做。
这是我的代码示例（使用通用输入数据）：
#--------------spline interpolation
x <- c(1973:2014)
population <- seq(500000, 600000, length.out=42)
list <- spline(x, population, n=4*length(x), method="fmm",
xmin=min(x), xmax=max(x), ties=mean)
x_spline <- list$x
pop_spline <- list$y
如何定义样条是"每季度"计算的，换句话说，在 1973.25、1973.5、1973.75、1974 等？很抱歉不是统计学专家：将年度数据"转换"为季度数据的最佳方法是什么："fmm"、"自然"、"周期"、"monoH.FC"或"hyman"？假设人口增长在一年中分布均匀。
为什么不使用splinefun：
func = splinefun(x=x, y=population, method="fmm", ties = mean)
然后定义所需的预测点：
func(seq(1973, 2014, 0.25))