Question

我正在运行（至少尝试）R中的gibbs采样器，并且在查看数据时遇到问题。

我正在使用R Studio

以下是我的代码：

###Step 2
###Draw from truncated normal
for(i in 1:length(cens.list[,1])){

  row.i <- as.numeric(cens.list[i,1])

  mu.i <- as.numeric(x.mat[row.i,] %*% beta)

  bio.set[row.i,var.name] <- rtnorm(1,mean=mu.i,sd=sqrt(var),lower=-Inf,upper=dl)
}

test <- bio.set[which(row.names(bio.set) %in% cens.list[,1]),]

bio.set包含我想要输入的变量和其他信息

cens.list包含要估算的删失值列表及其对应的行信息。

dl是检测限。由于低于此值，所估算的观察结果缺失

问题：当我尝试查看数据时（使用View()或edit()），估算的值都显示相同的数字（检测限）。无论我是通过测试数据框还是通过bio.set本身查看它们，都会发生这种情况。

但是，如果我通过键入

单独查看值

bio.set[995,var.name] #Where I set the row number to be one of the imputed rows.

它显示了正确绘制的值。

我可以继续假设数据存储正确并且只是一个查看器问题吗？或者有什么我做错了吗？

Answer 1

感谢您的回复。当我弄清楚我的问题是什么时，我正在努力创建一个可重现的例子。

基本上发生的事情是cens.list包含要估算的值的行号。创建它之后，由于某种原因删除了一些行，但前面的行号由R

存储

在我提供的代码中，您可以在循环中看到我引用了特定于数据帧当前状态的行号而不是row.names。

为了解决这个问题，我将使用id变量，因为我应该一直这样做。

在将来的问题中，我将确保我的代码是可重现的。

数据查看问题

1 个答案: