我正在使用rollapply函数来创建滚动logit估计。到目前为止,该函数正在正常工作,但我现在想要为数据的每次迭代提取该特定窗口的最后一个值,以使用它们来预测下一个下一个周期。我试图使用一个简单的应用函数,但它似乎不起作用。这是我的代码:
rollapply(zoo(DATASET),
width=300,
FUN = function(Z)
{
t = glm(formula= y ~ ( Sentiment-SentimentLag) + (Sentiment-SentimentLag10) + Volume , data = as.data.frame(Z),family=binomial(link='logit'));
new <- coeftest(t,vcov.=vcovHAC(t,sandwich = T))
pval <- new[,4]
#### this line does only return single values and not the vector I want
lastOb <- as.matrix(apply(Z,1,max))
},
by.column=FALSE, align="right" )
我做错了什么?