用R将年度数据分解(插值)为季度数据
用R将年度数据分解(插值)为季度数据
提问于 2016-08-17 13:18:43
我的数据采用以下形式：
df <- data.frame(year=c(1992:2015), share=c(31.9, 36.8, 38.2, 39.9, 36.3, 36.5, 35.6, 35.2, 34.8, 33.2, 33.5, 34.6,
36.3, 36.2, 38.1, 37.2, 35.9, 33.2, 36.9, 36.0, 33.9, 33.7, 34.3, 35.1))该数据由平均年值组成。假设变化是线性的，我想将数据转换为平均季度值。
我试过使用library("tempdisagg")，但我无法让它工作。我还在stats中尝试了一些使用library("splines")的其他方法，但都没有效果。也许是因为我处理的日期和时间序列格式不正确。
为了澄清，预期的输出如下所示：
y_q share
1992q1 values
1992q2 values
1992q3 values
1992q4 values
1993q1 values
1993q2 values
1993q3 values
1993q4 values
1994q1 values
1994q2 values
1994q3 values
1994q4 values任何帮助都将不胜感激。
回答 2
Stack Overflow用户
发布于 2016-08-17 13:50:35
“近似”函数可以对输出所需的值执行线性插值。
#create sequence of dates
quarter<-seq(as.Date("1992-01-01"), as.Date("2015-01-01"), by="quarter")
#create the linear interploated values
estshare<-approx(df$share, n=length(quarter))
newdf<-data.frame(quaters<-quarter, share<-estshare$y)在这种情况下，我假设股价是从一年的第一天开始。如果它是平均，您可能想改变日期序列的开始日期和结束日期，“季度”从一年的第一次到年中。
Stack Overflow用户
发布于 2016-08-18 08:21:05
通过线性年度插值，您可以在基本包中进行以下操作：
annual <- data.frame(year=c(1992:2015), share=c(31.9, 36.8, 38.2, 39.9, 36.3, 36.5, 35.6, 35.2, 34.8, 33.2, 33.5, 34.6, 36.3, 36.2, 38.1, 37.2, 35.9, 33.2, 36.9, 36.0, 33.9, 33.7, 34.3, 35.1))
# Annual delta for interpolation
annual$delta<-c(NA, diff(annual$share,1))
# Quarterly table
ref<-data.frame(quarter=paste0("Q", 1:4), nb=1:4)
quart<-merge(annual, ref)
quart<-quart[order(quart$year, quart$quarter),]
# quarterly value calculation with evolution (loop)
quart$quarterly<-NA
quart$quarterly[1:4]<-quart$share[1:4]/4
for (i in (2:dim(annual)[1])) {
quart$quarterly[quart$year==annual$year[i]] <- sum(quart$quarterly[quart$year==annual$year[i-1]])/4+ (quart$delta[quart$year==annual$year[i]] * quart$nb[quart$year==annual$year[i]])/10
}
# /10 : /(1+2+3+4)
# Check :
summary(annual$share == aggregate(quarterly ~ year, data=quart, FUN=sum)[,2])
plot(quart$quarterly, typ="l")它可能很难看，但我更喜欢使用基本的函数来理解底层代码，这样我就可以根据其他情况调整代码。
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/38997849
复制相关文章
相似问题
腾讯云开发者
