Как создать файл матрицы t x nk для Gauss

Question

Я хотел бы использовать эконометрический тест данных панели, который реализован в статистическом программном обеспечении Гаусса. С пакетом R MASS::write.matrix Мне удалось создать файл ASCII и прочитать файл из Гаусса. Это хорошо работало для t x n матрицы. Но мне интересно, как экспортировать матрицу t x nk. Столбцы nk просто добавляются друг к другу?

Answer 1

Используя пример набора данных из пакета Plm, вот как вы можете изменить набор данных в АТХ формате пк матрицы:

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)];