同事!我有面板数据:
Company year Beta NI Sales Export Hedge FL QR AT Foreign
1 1 2010 -2.2052800 293000 1881000 78.6816 0 23.5158 1.289 0.6554 3000
2 1 2011 -2.2536069 316000 2647000 81.4885 0 21.7945 1.1787 0.8282 22000
3 1 2012 0.3258693 363000 2987000 82.4908 0 24.5782 1.2428 0.813 -11000
4 1 2013 0.4006030 549000 4546000 79.4325 0 31.4168 0.6038 0.7905 71000
5 1 2014 -0.4508811 348000 5376000 79.2411 0 37.1451 0.6563 0.661 -64000
6 1 2015 0.1494696 355000 5038000 77.1735 0 …晚上好！现在我用外部回归器测试了 GARCH 模型：
library("tseries")
library("xts")
library("rugarch")
Rim <- read.csv("Rim.csv", header=T, sep=";")
Rcm <- read.csv("Rcm.csv", header=T, sep=";")
n=119
spec <- ugarchspec(variance.model = list(model = "sGARCH", garchOrder = c(1, 1), submodel = NULL, external.regressors = NULL, variance.targeting = FALSE), mean.model = list(armaOrder = c(0, 0), external.regressors = matrix(Rcm[,1])), distribution.model = "norm", start.pars = list(), fixed.pars = list())
my_lms=lapply(1:n,function(x) ugarchfit(spec=spec,data=Rim[,x],solver.control=list(trace=0)))
昨天我测试了这个模型，R 显示了这个模型的结果。今天我没有更改代码，但现在 R 给出了错误：
Warning messages:
1: In .sgarchfit(spec = spec, data = data, out.sample = out.sample, :
ugarchfit-->warning: solver failer to converge. …我有一个问题,如何构建这种类型的矩阵,节省更多的时间.任务是用户输入自然数(n = 1,2,3,4,5 ......)并且R必须打印这种类型的矩阵(在n = 4的矩阵类型下面).我观察到第一列和第二列与第6列和第7列对称,并且足以正确打印左侧的第一列和第二列.我还观察到第3列是第二列+ c(0,0,1,1,1,0,0),对应1 + 1 + 1 = 3 - 该列的数量.但是我不明白,对于一般情况,这个矩阵的算法是什么(对于n = 6的fe,dim是2n-1 x 2n-1(nrow x ncol).它是构造这个矩阵的最简单的变体吗？按列或函数外部列允许简化该任务？
1 1 1 1 1 1 1
1 2 2 2 2 2 1
1 2 3 3 3 2 1
1 2 3 4 3 2 1
1 2 3 3 3 2 1
1 2 2 2 2 2 1
1 1 1 1 1 1 1