代码创建的矢量除最后一项外仅具有NA值。
我尝试了下面显示的代码。
data <- subset(data, YYYY > 1975)
years <- unique(data$YYYY)
q95 <- rep(1,124)
for(m in length(years)){
data_sub <- subset(data, YYYY==years[m])
q95[m] <- quantile(data_sub$RGES,probs=1-0.95)
}
我期望一个向量,该向量在每个循环中分配的值从1到124。