如何在3D曲面上添加置信区间？

Question

我有一个矩阵40x40，该矩阵是通过使用库akima进行插值获得的值创建的，以创建3D曲面。

我使用蒙特卡洛模拟从预测值中估计出95％的CI，现在我想将0年的值添加到3D图形中。

我做错了什么，我不知道如何绘制垂直线来表示配置项。 我的台词看起来像这样：

我想在此图像上使用CI：

这是我的数据保管箱链接，因为它比此处允许发布的空间长：https://www.dropbox.com/s/c6iyd2r00k5jbws/data.rtf?dl=0

和我的代码：

    persp(xyz,theta = 45, phi = 25,border="grey40", ticktype = "detailed", zlim=c(0,.8))->res2

   y.bin <- rep(1,25)
    x.bin <- seq(-10,10,length.out = 25)

points (trans3d(x.bin, y.bin, z = y0, pmat = res2), col = 1, lwd=2)
lines (trans3d(x.bin, y.bin, z = LCI, pmat = res2), col = 1, lwd=2)
lines (trans3d(x.bin, y.bin, z = UCI, pmat = res2), col = 1, lwd=2)

Answer 1

问题在于，上下置信区间被绘制为一条直线。如果您以一定的间隔在这些点上循环，然后在该点的上限值和下限值之间绘制线，则该图看起来更接近您想要的值。请注意，示例y0中的点值并不是很多3d间隔。

# data from link is imported
persp(xyz,theta = 45, phi = 25,border="grey40", ticktype = "detailed", zlim=c(0,.8))->res2

y.bin <- rep(1,25)
x.bin <- seq(-10,10,length.out = 25)

# y0 points
points (trans3d(x.bin, y.bin, z = y0, pmat = res2), col = 1, lwd=2)

# lines between upper and lower CIs for each location
for(i in 1:length(LCI)){
  lines (trans3d(rep(x.bin[i],2), rep(y.bin[i],2), z = c(LCI[i],UCI[i]), pmat = res2), col = 1, lwd=2)
}