r - Gauss の atx nk マトリックス ファイルを作成するにはどうすればよいですか

Question

ガウス統計ソフトウェアに実装されているパネル データの計量経済学的検定を使用したいと考えています。R パッケージMASS::write.matrixを使用して、ASCII ファイルを生成し、Gauss 内からファイルを読み取ることができました。これは、txn マトリックスでうまく機能しました。しかし、atx nk マトリックスをエクスポートする方法が気になります。nk 列は単純に互いに追加されますか?

Answer 1

plm パッケージのサンプル データセットを使用して、データセットを atx nk マトリックス形式で再形成する方法を次に示します。

library(dplyr)
library(tidyr)
library(plm) # For the example dataset
data("Produc", package = "plm")
spreadvariable <- c("pcap", "pc", "emp", "unemp")
gaussmatrixfile <- file.path(tempdir(),"gaussmatrix.prn")
gaussmatrixfile
Produc %>%
    select_("year", "state", .dots = spreadvariable) %>% 
    gather(variable, value, -year,-state) %>%
    unite(state_variable, state, variable) %>%
    spread(state_variable, value) %>%
    MASS::write.matrix(gaussmatrixfile)

このファイルは、Gauss から次のコマンドで読み取ることができます。

load datax[t,n*k+1]     = gaussmatrix.prn;        /* t x nk matrix */
/* Remove the first column with years*/
x = datax[.,2:cols(datax)];