Dear useRs,
I'm trying to simply fill in the area under a curve using RGL. Here' the set up: x < c(0.75,75.75,150.75,225.75,300.75,375.75,450.75,525.75,600.75,675.75, 0.5,50.5,100.5,150.5,200.5,250.5,300.5,350.5,400.5,450.5, 0.25,25.25,50.25,75.25,100.25,125.25,150.25,175.25,200.25,225.25) y < c(0.05,4.91,9.78,14.64,19.51,24.38,29.24,34.11,38.97,43.84, 0.1,9.83,19.56,29.29,39.02,48.75,58.48,68.21,77.94,87.67, 0.15,14.74,29.34,43.93,58.53,73.13,87.72,102.32,116.91,131.51) z < c(0.05,0.55,0.7,0.78,0.83,0.87,0.9,0.92,0.93,0.94, 0,0.32,0.59,0.77,0.87,0.93,0.96,0.98,0.99,1, 0,0.39,0.66,0.82,0.9,0.95,0.97,0.99,0.99,1) dat < data.frame(x = x, y = y, z = z, ID = c(rep(c(1,2,3),each=10))) plot3d(dat, type = "n", ylab = "", xlab = "", zlab = "", axes = F, ylim = c(0,200)) lines3d(dat[1:10,]) lines3d(dat[11:20,]) lines3d(dat[21:30,]) axes3d(edge = c("x", "y+", "z"), nticks = 5, ylim = c(0,200)) bbox3d(color = c("black", "white"), lit = F, back = "line") Any ideas/tips on how to do this? thanks in advance, Patrick 
Well I've figures out a painful way of doing it with triangles 3d, as suggested in another post by Duncan Murdoch.
The code is disgusting, but it works. x < c(0,0.75,75.75,150.75,225.75,300.75,375.75,450.75,525.75,600.75,675.75, 675.75,0, 0,0.5,50.5,100.5,150.5,200.5,250.5,300.5,350.5,400.5,450.5, 450.5,0, 0,0.25,25.25,50.25,75.25,100.25,125.25,150.25,175.25,200.25,225.25,225.25,0) y < c(0,0.05,4.91,9.78,14.64,19.51,24.38,29.24,34.11,38.97,43.84, 43.84,0, 0,0.1,9.83,19.56,29.29,39.02,48.75,58.48,68.21,77.94,87.67, 87.67,0, 0,0.15,14.74,29.34,43.93,58.53,73.13,87.72,102.32,116.91,131.51, 131.51,0) z < c(0,0.05,0.55,0.7,0.78,0.83,0.87,0.9,0.92,0.93,0.94,0.0,0, 0,0.32,0.59,0.77,0.87,0.93,0.96,0.98,0.99,1,0,0,0, 0,0.39,0.66,0.82,0.9,0.95,0.97,0.99,0.99,1,0,0,0) dat < data.frame(x = x, y = y, z = z, ID = c(rep(c(1,2,3),each=13))) x1 < x[1:13] y1 < y[1:13] z1 < z[1:13] x2 < x[14:26] y2 < y[14:26] z2 < z[14:26] x3 < x[27:39] y3 < y[27:39] z3 < z[27:39] t1 < triangulate(as(cbind(x1,z1), "gpc.poly")) t2 < triangulate(as(cbind(x2,z2), "gpc.poly")) t3 < triangulate(as(cbind(x3,z3), "gpc.poly")) yfit1 < predict(lm(y1~ x1 + z1), newdata = data.frame(x1 = t1[,1], z1 = t1[,2])) yfit2 < predict(lm(y2~ x2 + z2), newdata = data.frame(x2 = t2[,1], z2 = t2[,2])) yfit3 < predict(lm(y3~ x3 + z3), newdata = data.frame(x3 = t3[,1], z3 = t3[,2])) plot3d(dat, type = "n", ylab = "", xlab = "", zlab = "", axes = F, ylim = c(0,200)) axes3d(edge = c("x", "y+", "z"), nticks = 5, ylim = c(0,200)) bbox3d(color = c("black", "white"), lit = F, back = "line") lines3d(x1, y1, z1) lines3d(x2, y2, z2) lines3d(x3, y3, z3) triangles3d(x= t1[,1], y= yfit1, z= t1[,2], lit = F, col = "white", front = "fill", lwd = 10) triangles3d(x= t2[,1], y= yfit2, z= t2[,2], lit = F, col = "white", front = "fill", lwd = 10) triangles3d(x= t3[,1], y= yfit3, z= t3[,2], lit = F, col = "white", front = "fill", lwd = 10)

Free forum by Nabble  Edit this page 