2

我在使用该scatter3d函数时遇到了一些问题，特别是似乎不可能向函数规定轴的限制应该是什么。这是我的基本绘图功能，其中包含我正在使用的集合中的一些数据点：

library(rgl)
library(plot3D)
library(car) 

df <- data.frame(meanX = c(147.34694, 173.89244, 135.73004, 121.93766,
                           109.72152,  92.53709, 165.46588, 169.77744,
                           127.01796,  99.34347),
                 meanY = c(140.40816, 110.99128, 134.56023, 164.18703,
                           166.04051, 155.97329, 105.29377, 104.42683,
                           130.17066, 155.99696),
                 avgDist = c(40.788118, 12.957329, 14.24348, 39.10424,
                             34.694258, 25.532335, 21.491695,23.528944,
                             9.309201,  31.916879))

car::scatter3d(x = df$meanX, y = df$meanY, z = df$avgDist, surface = FALSE)

这绘制得很好，但设置了默认轴范围

xlim = c(90,200) 
ylim = c(100,200)
zlim = c(9,40)

问题是，这似乎完全无法改变。我知道RGL不允许裁剪，我不希望它在这里排除任何数据点，但是轴范围比每个轴的最大值高很多

max(df$meanX)
[1] 173.8924

max(df$meanY)
[1] 166.0405

max(df$avgDist)
[1] 40.78812

我想做的是将 X 和 Y 的轴限制从 70 设置为 185，据我所知，这应该可以通过以下代码实现：

car::scatter3d(x = df$meanX, y = df$meanY, z = df$avgDist, surface = FALSE,
          xlim = c(70,185), ylim = c(70,185))

但这只会产生相同的情节（没有错误或警告）。有人知道如何手动设置这些轴吗？

4

1 回答 1

1

在axis.scales = TRUE（现在是默认值）的情况下，用于设置每个轴的最小和最大标签和刻度的代码由非导出car:::nice函数选择。您可以调用axis.scales=FALSE，然后重做添加轴刻度和标签的工作。

with( df,  scatter3d( meanX , meanY, avgDist, axis.scales = FALSE, 
           xlim = c(70,185), ylim = c(70,185), zlim = c(9,40) ))

或者，您可以尝试重写代码以检查是否存在任何 xlim、ylim 或 zlim 参数，并将这些值替换为 car:::scatter3d 中这段代码设置的值：

car:::scatter3d
#------- code near the top of the console display------
if (axis.scales) {
    lab.min.x <- nice(minx)
    lab.max.x <- nice(maxx)
    lab.min.y <- nice(miny)
    lab.max.y <- nice(maxy)
    lab.min.z <- nice(minz)
    lab.max.z <- nice(maxz)
    minx <- min(lab.min.x, minx)
    maxx <- max(lab.max.x, maxx)
    miny <- min(lab.min.y, miny)
    maxy <- max(lab.max.y, maxy)
    minz <- min(lab.min.z, minz)
    maxz <- max(lab.max.z, maxz)
    min.x <- (lab.min.x - minx)/(maxx - minx)
    max.x <- (lab.max.x - minx)/(maxx - minx)
    min.y <- (lab.min.y - miny)/(maxy - miny)
    max.y <- (lab.max.y - miny)/(maxy - miny)
    min.z <- (lab.min.z - minz)/(maxz - minz)
    max.z <- (lab.max.z - minz)/(maxz - minz)
}

我认为只允许lim值在值的内部进行修改是可以的nice，但是这段代码出错了，所以需要逻辑来检测 lims 的缺失：

    lab.min.x <- max( nice(minx), xlim[1])
    lab.max.x <- min( nice(maxx), xlim[2])
    lab.min.y <- max( nice(miny), ylim[1])
    lab.max.y <- min( nice(maxy), ylim[2])
    lab.min.z <- max( nice(minz), zlim[1])
    lab.max.z <- min( nice(maxz), zlim[2])
于 2016-05-25T20:16:19.790 回答