我通过R中的BRugs
包为OpenBUGS中的空间建模创建了一些模拟疾病计数数据。
我有一个现有的.txt
文件,其中包含在BRugs
中运行模型所需的大部分信息。这是这个文本文件的样子:
list(N = 44,
O = c(4, 8, 3, 6, 2, 7, 7, 2, 10, 13, 5, 7, 9, 3, 3, 1, 2, 3, 19, 7, 5, 4,
9, 2, 4, 6, 14, 13, 6, 7, 9, 3, 4, 10, 10, 8, 23, 24, 18, 10, 6, 4, 17, 7),
E = c(7.20902956, 7.81436792, 3.46693788, 7.04393728, 7.37412184,
4.45749156, 6.3285374, 4.45749156, 6.3285374, 6.76878348, 9,...)
etc....
),
sumNumNeigh = 212
)
我想要完成的是将上面O
文件中的.txt
向量替换为我在R中生成的m
模拟疾病计数向量。这种方式,我可以使用m
个新的.txt文件与BRugs
一起运行我的模型。
例如,让我们说这个新的疾病计数向量被称为O_sim
:
[1] 6 4 6 9 9 6 11 6 7 6 11 9 7 10 9 10 9 3 10 4 6 8 8 4 6 9 8 6 9 8 7 8 5 10 8 6 4 4 5 5 8
[42] 10 8 11
我怎么能创建一个新的.txt文件,现在看起来像这样:
list(N = 44,
O = c(6,4,6,9,9,6,11,6,7,6,11,9,
7,10,9,10,9,3,10,4,6,8,8,4,6,
9,8,6,9,8,7,8,5,10,8,6,4,4,5,5,
8,10,8,11),
E = c(7.20902956, 7.81436792, 3.46693788, 7.04393728, 7.37412184,
etc...
),
sumNumNeigh = 212
)