Thank you very much, Sarah, for your help. The colvar argument was certainly what I needed. I will follow your suggestion. Cheers, Jorge.-
On Thu, Jul 14, 2016 at 1:29 PM, Sarah Goslee <sarah.gos...@gmail.com> wrote: > I assume you want some variant of: > with(d0, scatter3D(X, Z, Y, bty = "b2", col = cols, colvar = c(1, 2, 3), > pch = 20, cex = 4, ticktype = "detailed", colkey = FALSE, phi = > 20, theta = -140, zlim = c(-14, 14), xlim = c(-14, 14), ylim = c(-14, 14), > xlab = "X", ylab = "Z", zlab = "Y", CI = CI2)) > mtext('Experiment 2', line = -0.5, side = 3, cex = 1.3) > > You should reread the help for scatter3D, paying close attention to > col and colvar arguments. The default value of colvar is z, and that's > what scatter3D() was using. Setting colvar to your chosen levels gives > you the result you expect. > > Thanks for the complete reproducible example. I wouldn't have even > looked at your problem without something to test. > > Sarah > > On Thu, Jul 14, 2016 at 1:55 PM, Jorge I Velez <jorgeivanve...@gmail.com> > wrote: > > Dear R-help, > > > > I am using the plot3D package to produce 3D spheres along with 95% CIs > > distinguishing each sphere with a predefined colour (see the reproducible > > example at the end). > > > > I have been successful in producing a similar plot using a different data > > set (kindly see > > https://www.dropbox.com/s/snpgiqgqaiqgpiv/exp1combined.pdf?dl=0), but in > > this case I cannot arrange the colours as desired. Specifically, the > > colours obtaind do not correspond to those in the "cols" object below > (kindly > > see https://www.dropbox.com/s/nhljsuare84g811/test_exp2.pdf?dl=0). > > > > Any help/tips would be greatly appreciated. > > > > Thank you very much in advance. > > > > Best regards, > > Jorge Velez.- > > > > > > ## package needed > > if(!require(plot3D)) install.packages("plot3D") > > require(plot3D) > > > > ## data to be plotted > > d0 <- structure(list(X = c(2.5, -2, 1), Z = c(-3.5, 4, -1), Y = c(8, > > -8.5, -1)), .Names = c("X", "Z", "Y"), row.names = c("high", > > "low", "medium"), class = "data.frame") > > d0 > > > > ## confidence intervals to be added > > CI2 <- structure(list(z = structure(c(2, 3, 4, 2, 3, 4), .Dim = c(3L, > > 2L), .Dimnames = list(c("high", "low", "medium"), c("97.5%", > > "97.5%"))), y = structure(c(4, 5, 5, 4, 5, 5), .Dim = c(3L, 2L > > ), .Dimnames = list(c("high", "low", "medium"), c("97.5%", "97.5%" > > ))), x = structure(c(3, 7, 5, 3, 7, 5), .Dim = c(3L, 2L), .Dimnames = > list( > > c("high", "low", "medium"), c("97.5%", "97.5%"))), alen = 0, > > lwd = 2), .Names = c("z", "y", "x", "alen", "lwd")) > > > > ## colours I would like to have > > cols <- c("#0080ff", "#ff00ff", "darkgreen") > > > > # this produces the 3D plot, but the colours are not properly assigned > > with(d0, scatter3D(X, Z, Y, bty = "b2", col = cols, > > pch = 20, cex = 4, ticktype = "detailed", colkey = FALSE, phi > = > > 20, theta = -140, zlim = c(-14, 14), xlim = c(-14, 14), ylim = c(-14, > 14), > > xlab = "X", ylab = "Z", zlab = "Y", CI = CI2)) > > mtext('Experiment 2', line = -0.5, side = 3, cex = 1.3) > > > [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.