为什么我在Eviews和R（rugarch）中估计GARCH-M模型会得到非常不同的结果？

Question

我不知道是否允许在StackExchange上的不同网站上发布相同的问题。我已将它贴在Cross Vlidated上并且没有答案，所以我想在这里试一试。

这是原始question。如果有人能解决我的问题，我真的很感激。

我在R中的代码如下：

#get data
re <- read.table("return.csv",sep=",",header=TRUE)
re[,1] <- as.character(re[,1])
re[,1] <- as.Date(re[,1])

#create dummy variable
re$T <- re[,1]>as.Date("2010-03-31")
re$T <- as.numeric(re$T)
T <- as.matrix(re$T)

library("xts")
xts <- as.xts(re[,-1],order.by=re[,1])

library("rugarch")
#GARCH specification
garchspec <- ugarchspec(variance.model = list(model = "sGARCH", garchOrder = c(1, 1), 
                        submodel = NULL, external.regressors = T, 
                        variance.targeting = FALSE), 
                        mean.model = list(armaOrder = c(0, 0), include.mean = TRUE,
                        archm = TRUE, archpow = 1, arfima = FALSE, 
                        external.regressors = NULL, archex = FALSE), 
                        distribution.model = "ged",
                        start.pars = list(), fixed.pars = list())

#fitting
fit <- ugarchfit(spec=garchspec, data=xts[,1], out.sample = 0,solver="solnp",
                 solver.control = list(trace=0), fit.control = 
                 list(stationarity = 1, fixed.se = 0, scale = 0, rec.init = 0.7))
show(fit)

它有什么错误吗？