Я немного погуглил и не могу найти документацию по этому вопросу. Я пытаюсь оценить допустимую модель обобщенных наименьших квадратов (FGLS) на данных временных рядов поперечного сечения в R. Например:
library(nlme)
foo <- gls(Y ~ factor(panel_ID) + X1 + X2, data = myData,
correlation=corARMA(p=1), method='ML', na.action=na.pass)
Когда я запускаю это (мой фрейм данных довольно большой, поэтому я не включаю его здесь), я получаю следующую ошибку:
# Error in array(c(X, y), c(N, ncol(X) + 1), list(row.names(dataMod), c(colnames(X), :
# length of 'dimnames' [1] not equal to array extent</code>
Кто-нибудь достаточно знаком с внутренней работой gls
или пакета nlme
в целом, чтобы сказать мне, что я делаю здесь неправильно? Или предложить другой способ сделать это (я также попробовал пакет plm
)?