R中的read.coda函数有问题

Question

我无法使用read.coda函数中的coda函数读取马尔可夫链蒙特卡罗（MCMC）输出，以便在R中进行贝叶斯诊断。我的问题是mcmc表示数据正在填充NA。这是函数返回的内容：

Markov Chain Monte Carlo (MCMC) output:
Start = 1 
End = 6 
Thinning interval = 1 
    Rho Sigma2
  1 0.204     NA
  2 0.205     NA
  3 0.206     NA
  4    NA  0.147
  5    NA  0.148
  6    NA  0.149

我的预期输出是

Markov Chain Monte Carlo (MCMC) output:
Start = 1 
End = 6 
Thinning interval = 1 
    Rho Sigma2
  1 0.204     0.147
  2 0.205     0.148
  3 0.206     0.149

格式化需要索引文件（.txt）告诉您每个变量的名称和行范围，链文件（.txt）包含变量值和迭代次数。这是一个基本的工作示例，它复制了我遇到的问题。我的index.file看起来像是：

Rho     1   3
Sigma2  4   6

其中第1-3行对应于Rho的值，对于Sigma2为4-6，而我的out.file如下：

1   0.204
2   0.205
3   0.206
4   0.147
5   0.148
6   0.149

有谁知道为什么我没有返回预期的输出？