私はかなりグーグルで検索しましたが、これに関するドキュメントが見つかりません。R の断面時系列データで実行可能な一般化最小二乗 (FGLS) モデルを推定しようとしています。たとえば、次のようになります。
library(nlme)
foo <- gls(Y ~ factor(panel_ID) + X1 + X2, data = myData,
correlation=corARMA(p=1), method='ML', na.action=na.pass)
これを実行すると (データ フレームが非常に大きいため、ここには含めません)、次のエラーが発生します。
array(c(X, y), c(N, ncol(X) + 1), list(row.names(dataMod), c(colnames(X), : 'dimnames' の長さ [1] が等しくない) のエラー配列の範囲まで
ここで私が間違っていることを教えてくれる、パッケージの内部の仕組みglsや一般的なパッケージに精通している人はいますか? または、これについて別の方法を提案しますか（パッケージnlmeも試しました）？plm