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

Question

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

Answer 1

Используя пример набора данных из пакета plm, вы можете изменить форму набора данных в матричном формате a t x 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)];