线框图（格子包）中的负长度向量是什么意思？

Question

我想使用 lattice 包在 R 中绘制线框。但是，我收到以下错误消息“不允许使用数据包 1 负长度向量的错误”。数据如下所示：

> result_mean
  experiment alpha beta   packet
1          0   1.0    1 3.000000
2          0   1.1    1 2.571429

创建数据的命令如下

  png(file=paste("foobar.png"),width=1280, height=1280);
    plot <- wireframe(result_mean$packet ~ result_mean$alpha * result_mean$beta,       
      data=result_mean, scales = list(arrows=FALSE, cex= .45, col = "black", font = 3),
      drape = TRUE, colorkey = TRUE, main = "Foo",
      col.regions = terrain.colors(100),
      screen = list(z = -60, x = -60),
      xlab="alpha", ylab="beta", zlab="mean \npackets");
  print(plot);
  dev.off();

我想知道这个错误消息是什么意思，是否有调试的好方法？

提前致谢！

Answer 1

调试点阵图形有点困难，因为 (a) 代码复杂且具有多层性，并且 (b) 错误以一种难以拦截的方式被捕获。但是，您至少可以通过某种方式诊断问题。

首先创建一个最小的例子。我怀疑你的问题是你的数据落在一行上，所以我创建了看起来像这样的数据：

d <- data.frame(x=c(1,1.1),
                y=c(1,1),
                z=c(2,3))

library(lattice)
wireframe(z~y*x,data=d)

现在确认全三维数据（定义平面的数据）工作正常：

d2 <- data.frame(expand.grid(x=c(1,1.1),
                             y=c(1,1.1)),
                 z=1:4)
wireframe(z~y*x,data=d2) 

所以问题是真的——你是否打算画一个由两个点组成的线框？如果是这样，你想在情节中出现什么？你可以稍微修改一下，将y 的值设置为有一点点不同——我试过了，但没有出现线框（但也没有错误）。

编辑：我做了更多的跟踪，使用各种debug() 咒语（并搜索lattice 包的源代码和R 本身的“负长度”）来推断如下：在一个名为 lattice:::panel.3dwire 的函数中，有一个对 C 函数 wireframePanelCalculations 的调用，您可以在 https://r-forge.r-project.org/scm/viewvc.php/pkg/src/threeDplot.c?view=markup&root=lattice 看到它

在这个函数中：

nh = (nx-1) * (ny-1) * ng; /* number of quadrilaterals */
sHeights = PROTECT(allocVector(REALSXP, nh));

在这种情况下，nx 为零，因此此代码要求 R 分配一个负长度向量，这就是错误的来源。

不过，在这种情况下，我认为诊断比显式调试更有用。