如何創建高斯矩陣文件

Question

我想使用在高斯統計軟件中實現的面板數據計量經濟學測試。使用R包MASS::write.matrix我已經設法生成一個ASCII文件並從高斯內讀取文件。這對t x n矩陣很有效。但我想知道如何導出一個t×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) 

這個文件然後可以從高斯與

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