r语言 - 如何为高斯创建 t x nk 矩阵文件

Question

我想使用在高斯统计软件中实现的面板数据计量经济学测试。使用 R 包MASS::write.matrix我设法生成了一个 ASCII 文件并从高斯中读取该文件。这适用于 t x n 矩阵。但是我想知道如何导出 t x nk 矩阵。nk 列会简单地相互追加吗？

Answer 1

使用 plm 包中的示例数据集，下面介绍了如何以 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)];