Hey guys, Not sure if I encountered a bug with the scatterplot3d function. Here's the calls I made: s3d1 <- scatterplot3d(TotLogDisttenp,TotDifftenp,TotMeasuredRSLtenp,pch=16,highlight.3d=TRUE,angle=40,type="h",main="MRSL ~ LogDist + Diff"); s3d1$plane3d(fitols); s3d1 <- scatterplot3d(TotLogDisttenp,TotDifftenp,TotMeasuredRSLtenp,pch=16,highlight.3d=TRUE,angle=130,type="h",main="MRSL ~ LogDist + Diff"); s3d1$plane3d(fitols); s3d1 <- scatterplot3d(TotLogDisttenp,TotDifftenp,TotMeasuredRSLtenp,pch=16,highlight.3d=TRUE,angle=210,type="h",main="MRSL ~ LogDist + Diff"); s3d1$plane3d(fitols); s3d1 <- scatterplot3d(TotLogDisttenp,TotDifftenp,TotMeasuredRSLtenp,pch=16,highlight.3d=TRUE,angle=310,type="h",main="MRSL ~ LogDist + Diff"); s3d1$plane3d(fitols); Essentially four plots showing the data from different angles. This includes the fit plane. The first two graphs make sense, but for the latter two, the fit plane is not making sense. Take a look at the attached png. Is it a bug? Thanks, Vivek -------------- next part -------------- A non-text attachment was scrubbed... Name: ols_3dplot_diff_up.png Type: image/png Size: 69163 bytes Desc: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20090803/f7edea95/attachment-0002.png>
Not sure you can call it a bug when the help page says for angles in range (180, 360) that some support functions may "not work properly". On Aug 3, 2009, at 2:12 PM, Vivek Ayer wrote:> Hey guys, > > Not sure if I encountered a bug with the scatterplot3d function. > Here's the calls I made: > > s3d1 <-scatterplot3d(TotLogDisttenp,TotDifftenp, TotMeasuredRSLtenp, > pch=16,highlight.3d=TRUE,angle=40, type="h",main="MRSL ~ LogDist + > Diff"); > s3d1$plane3d(fitols);> s3d1 <- scatterplot3d(TotLogDisttenp,TotDifftenp, > TotMeasuredRSLtenp,pch=16, highlight. > 3d=TRUE,angle=130,type="h",main="MRSL~ LogDist + Diff"); > s3d1$plane3d(fitols); > s3d1 <- scatterplot3d(TotLogDisttenp,TotDifftenp, > TotMeasuredRSLtenp,pch=16, highlight. > 3d=TRUE,angle=210,type="h",main="MRSL ~ LogDist + Diff"); > s3d1$plane3d(fitols);I suspect s3d1$plane3d(fitols) uses one of the futzed functions.> s3d1 <- scatterplot3d(TotLogDisttenp,TotDifftenp, > TotMeasuredRSLtenp,pch=16, highlight. > 3d=TRUE,angle=310,type="h",main="MRSL~ LogDist + Diff"); > s3d1$plane3d(fitols); > > Essentially four plots showing the data from different angles. This > includes the fit plane. The first two graphs make sense, but for the > latter two, the fit plane is not making sense. > > Take a look at the attached png. > > Is it a bug?David Winsemius, MD Heritage Laboratories West Hartford, CT
I finally found some time to look at it: Yes, it does not work properly for angles > 180 degrees. Will try to find a fix and make a new release soon. Best wishes, Uwe Ligges Vivek Ayer wrote:> Hey guys, > > Not sure if I encountered a bug with the scatterplot3d function. > Here's the calls I made: > > s3d1 <- scatterplot3d(TotLogDisttenp,TotDifftenp,TotMeasuredRSLtenp,pch=16,highlight.3d=TRUE,angle=40,type="h",main="MRSL > ~ LogDist + Diff"); > s3d1$plane3d(fitols); > s3d1 <- scatterplot3d(TotLogDisttenp,TotDifftenp,TotMeasuredRSLtenp,pch=16,highlight.3d=TRUE,angle=130,type="h",main="MRSL > ~ LogDist + Diff"); > s3d1$plane3d(fitols); > s3d1 <- scatterplot3d(TotLogDisttenp,TotDifftenp,TotMeasuredRSLtenp,pch=16,highlight.3d=TRUE,angle=210,type="h",main="MRSL > ~ LogDist + Diff"); > s3d1$plane3d(fitols); > s3d1 <- scatterplot3d(TotLogDisttenp,TotDifftenp,TotMeasuredRSLtenp,pch=16,highlight.3d=TRUE,angle=310,type="h",main="MRSL > ~ LogDist + Diff"); > s3d1$plane3d(fitols); > > Essentially four plots showing the data from different angles. This > includes the fit plane. The first two graphs make sense, but for the > latter two, the fit plane is not making sense. > > Take a look at the attached png. > > Is it a bug? > > Thanks, > Vivek > > > ------------------------------------------------------------------------ > > > ------------------------------------------------------------------------ > > ______________________________________________ > R-help at r-project.org mailing list > 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.