2
在R中,函数ar.yw如何估计方差?具体来说,数字“var.pred”来自哪里?它似乎不是来自通常的YW估计的方差,也不是平方残差的总和除以df(即使对于df应该是什么意见不一致,没有一个选项给出了等同于var.pred的答案) 。是的,我知道有比YW更好的方法;试图弄清楚R在做什么。ar.yw如何估计方差
set.seed(82346)
temp <- arima.sim(n=10, list(ar = 0.5), sd=1)
fit <- ar(temp, method = "yule-walker", demean = FALSE, aic=FALSE, order.max=1)
## R's estimate of the sigma squared
fit$var.pred
## YW estimate
sum(temp^2)/10 - fit$ar*sum(temp[2:10]*temp[1:9])/10
## YW if there was a mean
sum((temp-mean(temp))^2)/10 - fit$ar*sum((temp[2:10]-mean(temp))*(temp[1:9]-mean(temp)))/10
## estimate based on residuals, different possible df.
sum(na.omit(fit$resid^2))/10
sum(na.omit(fit$resid^2))/9
sum(na.omit(fit$resid^2))/8
sum(na.omit(fit$resid^2))/7
A guess ...这是预测包中的函数吗?你知道R包都是开源的吗? –
@ 42,不,这是“stats”包的一部分,即给我们“lm”和所有其他基本统计功能的相同包。因此,人们可能会认为,即使它是开源的,它会做一些合理的事情。 – chucky