Adding points to a wireframe: 'x and units must have length >0' error

sessionInfo()

On Oct 24, 2011, at 4:12 PM, Megan Bartlett wrote:

Hi David,

Thanks for the suggestion - I changed pts1 to pts, but I still got
the same error as before. Do you know what else I'm doing wrong?

I do not. The code that works for me is:

wireframe(z ~ y*x, gridd, drape = TRUE, colorkey= TRUE,
? ? ? ? ?scales = list(arrows = FALSE),
? ? ? ? ?pts = pts1,
? ? ? ? ?panel.3d.wireframe =
? ? ? ? ?function(x, y, z,xlim, ylim, zlim, xlim.scaled, ylim.scaled,
zlim.scaled, pts1, ...)

? ? ? ? ? ? ? ? ? ?{
? ? ? ? ? ? ?panel.3dwire(x = x, y = y, z = z,
? ? ? ? ? ? ? ? ? ? ? ? ? xlim = c(-4, 0),
? ? ? ? ? ? ? ? ? ? ? ? ? ylim = c(0, 40),
? ? ? ? ? ? ? ? ? ? ? ? ? zlim = c(-160, 0),
? ? ? ? ? ? ? ? ? ? ? ? ? xlim.scaled = c(-0.5, 0.5),
? ? ? ? ? ? ? ? ? ? ? ? ? ylim.scaled = c(-0.5, 0.5),
? ? ? ? ? ? ? ? ? ? ? ? ? zlim.scaled = c(-0.5, 0.5),
? ? ? ? ? ? ? ? ? ? ? ? ? ...)
? ? ? ?xx <- xlim.scaled[1] + diff(xlim.scaled) * (pts$x - xlim[1]) /
diff(xlim)
? ? ? ?yy <- ylim.scaled[1] + diff(ylim.scaled) * (pts$y - ylim[1]) /
diff(ylim)
? ? ? ?zz <- zlim.scaled[1] + diff(zlim.scaled) * (pts$z - zlim[1]) /
diff(zlim)
? ? ? ? ? ? ?panel.3dscatter(x = xx,
? ? ? ? ? ? ? ? ? ? ? ? ? ? ?y = yy,
? ? ? ? ? ? ? ? ? ? ? ? ? ? ?z = zz,
? ? ? ? ? ? ? ? ? ? ? ? ? ? ?xlim = xlim,
? ? ? ? ? ? ? ? ? ? ? ? ? ? ?ylim = ylim,
? ? ? ? ? ? ? ? ? ? ? ? ? ? ?zlim = zlim,
? ? ? ? ? ? ? ? ? ? ? ? ? ? ?xlim.scaled = xlim.scaled,
? ? ? ? ? ? ? ? ? ? ? ? ? ? ?ylim.scaled = ylim.scaled,
? ? ? ? ? ? ? ? ? ? ? ? ? ? ?zlim.scaled = zlim.scaled,
? ? ? ? ? ? ? ? ? ? ? ? ? ? ?...)
? ? ? ? ?})

--
David.

Thanks,

Megan

PS. New code is:

wireframe(z ~ y*x, gridd, drape = TRUE, colorkey= TRUE,
? ? ? ? ? scales = list(arrows = FALSE),
? ? ? ? ? pts = pts1,
? ? ? ? ? panel.3d.wireframe =
? ? ? ? ? function(x, y, z,xlim, ylim, zlim, xlim.scaled,
ylim.scaled, zlim.scaled, pts, ...)

? ? ? ? ? ? ? ? ? ? {
? ? ? ? ? ? ? panel.3dwire(x = x, y = y, z = z,
? ? ? ? ? ? ? ? ? ? ? ? ? ?xlim = c(-4, 0),
? ? ? ? ? ? ? ? ? ? ? ? ? ?ylim = c(0, 40),
? ? ? ? ? ? ? ? ? ? ? ? ? ?zlim = c(-160, 0),
? ? ? ? ? ? ? ? ? ? ? ? ? ?xlim.scaled = c(-0.5, 0.5),
? ? ? ? ? ? ? ? ? ? ? ? ? ?ylim.scaled = c(-0.5, 0.5),
? ? ? ? ? ? ? ? ? ? ? ? ? ?zlim.scaled = c(-0.5, 0.5),
? ? ? ? ? ? ? ? ? ? ? ? ? ?...)
? ? ? ? xx <- xlim.scaled[1] + diff(xlim.scaled) * (pts$x -
xlim[1]) / diff(xlim)
? ? ? ? yy <- ylim.scaled[1] + diff(ylim.scaled) * (pts$y -
ylim[1]) / diff(ylim)
? ? ? ? zz <- zlim.scaled[1] + diff(zlim.scaled) * (pts$z -
zlim[1]) / diff(zlim)
? ? ? ? ? ? ? panel.3dscatter(x = xx,
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? y = yy,
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? z = zz,
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? xlim = xlim,
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ylim = ylim,
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? zlim = zlim,
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? xlim.scaled = xlim.scaled,
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ylim.scaled = ylim.scaled,
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? zlim.scaled = zlim.scaled,
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ...)
? ? ? ? ? })


On Mon, Oct 24, 2011 at 12:04 PM, David Winsemius <dwinsemius at comcast.net

wrote:

On Oct 24, 2011, at 1:23 PM, Megan Bartlett wrote:

Hi,

I'm trying to follow the suggestions given by Deepayan Sarkar in this
message:

http://tolstoy.newcastle.edu.au/R/help/05/11/16135.html

to plot 3-D points on a wireframe plot. The problem is that I keep
getting a
partly formed plot- with the colored lattice visible but no axis
labels or
additional points- with the error message "error using packet 1, 'x'
and
'units' must have length >0". ?Does anyone know what I'm doing
wrong? The
code I've been using is:

# generate some fake data
data.frame(x = seq(-4, 0, 0.5), y = seq(0, 40, 5))-> df
expand.grid(x = df$x, y = df$y) -> gridd
(gridd$y* gridd$x) -> gridd$z

data.frame(x = runif(10, -4, 0), y = runif(10, 0, 40))-> pts1
pts1$z <- pts1$x*pts1$y

# plot

wireframe(z ~ y*x, gridd, drape = TRUE, colorkey= TRUE,
? ? ? ? scales = list(arrows = FALSE),
? ? ? ? pts = pts1,
? ? ? ? panel.3d.wireframe =
? ? ? ? function(x, y, z,xlim, ylim, zlim, xlim.scaled, ylim.scaled,
zlim.scaled, pts1, ...)

? ? ? ? ? ? ? ? ? {
? ? ? ? ? ? panel.3dwire(x = x, y = y, z = z,
? ? ? ? ? ? ? ? ? ? ? ? ?xlim = c(-4, 0),
? ? ? ? ? ? ? ? ? ? ? ? ?ylim = c(0, 40),
? ? ? ? ? ? ? ? ? ? ? ? ?zlim = c(-160, 0),
? ? ? ? ? ? ? ? ? ? ? ? ?xlim.scaled = c(-0.5, 0.5),
? ? ? ? ? ? ? ? ? ? ? ? ?ylim.scaled = c(-0.5, 0.5),
? ? ? ? ? ? ? ? ? ? ? ? ?zlim.scaled = c(-0.5, 0.5),
? ? ? ? ? ? ? ? ? ? ? ? ?...)
? ? ? xx <- xlim.scaled[1] + diff(xlim.scaled) * (pts1$x - xlim[1]) /
diff(xlim)
? ? ? yy <- ylim.scaled[1] + diff(ylim.scaled) * (pts1$y - ylim[1]) /
diff(ylim)
? ? ? zz <- zlim.scaled[1] + diff(zlim.scaled) * (pts1$z - zlim[1]) /
diff(zlim)
? ? ? ? ? ? panel.3dscatter(x = xx,

You assigned 'pts1' to 'pts', so you should have used 'pts' above:


? ? ? ? ? ? ? ? ? ? ? ? ? ? y = yy,
? ? ? ? ? ? ? ? ? ? ? ? ? ? z = zz,
? ? ? ? ? ? ? ? ? ? ? ? ? ? xlim = xlim,
? ? ? ? ? ? ? ? ? ? ? ? ? ? ylim = ylim,
? ? ? ? ? ? ? ? ? ? ? ? ? ? zlim = zlim,
? ? ? ? ? ? ? ? ? ? ? ? ? ? xlim.scaled = xlim.scaled,
? ? ? ? ? ? ? ? ? ? ? ? ? ? ylim.scaled = ylim.scaled,
? ? ? ? ? ? ? ? ? ? ? ? ? ? zlim.scaled = zlim.scaled,
? ? ? ? ? ? ? ? ? ? ? ? ? ? ...)
? ? ? ? })

Thanks so much!

-Megan Bartlett

? ? ? ?[[alternative HTML version deleted]]

David Winsemius, MD
West Hartford, CT


? ? ? ?[[alternative HTML version deleted]]