Print plot to pdf, jpg or any other format when using scatter3d error

On Jan 3, 2011, at 8:17 PM, Jurica Seva wrote:

Hi,

I have been trying to output my graphs to a file (jpeg, pdf, ps, it
doesnt matter) but i cant seem to be able to get it to output. I tried
a few things but none of them worked and am lost as what to do now. I
am using the scatter3d function, and it prints out the graphs on tot
he screen without any problems, but when it comes to writing them to a
file i cant make it work. Is there any other way of producing
3dimensional graphs (they dont have to be rotatable/interactive after
the print out)?

The code is fairly simple and is listed down :

#libraries
library(RMySQL)
library(rgl)
library(scatterplot3d)
library(Rcmdr)

##############################################################################

#database connection
mycon <- dbConnect(MySQL(),
user='root',dbname='test',host='localhost',password='')
#distinct sessions
rsSessionsU01 <- dbSendQuery(mycon, "select distinct sessionID from
actiontimes where userID = 'ID01'")
sessionU01 <-fetch(rsSessionsU01)
sessionU01[2,]

#user01 data
mycon <- dbConnect(MySQL(),
user='root',dbname='test',host='localhost',password='')
rsUser01 <- dbSendQuery(mycon, "select
a.userID,a.sessionID,a.actionTaken,a.timelineMSEC,a.durationMSEC,b.X,b.Y,b.Rel__dist_,b.Total_dist_
from `actiontimes` as a , `ulogdata` as b where a.originalRECNO =
b.RECNO and a.userID='ID01'")
user01 <- fetch(rsUser01, n= -1)
user01[1,1]

#plot loop

for (i in 1:10){

userSubset<-subset(user01,sessionID ==
sessionU01[i,],select=c(timelineMSEC,X,Y))
userSubset
x<-as.numeric(userSubset$X)
y<-as.numeric(userSubset$Y)
scatter3d(

#??? I thought the function was scatterplot3d()

x,y,userSubset$timeline,xlim = c(0,1280), ylim = c(0,1024),
zlim=c(0,1800000),type="h",main=sessionU01[i,],sub=sessionU01[i,])
tmp6=rep(".ps")

Why name it ".ps" when you are using pdf.dev()?

tmp7=paste(sessionU01[i,],tmp6,sep="")
tmp7
rgl.postscript(tmp7,"ps",drawText=FALSE)

When you want to get output to the file device, you need to "surround"
the plotting commands. the pdf call goes at the beginning of the loop
and the dev.odd at the end. It seems dangerous to specify the number to
dev.off since there might be no or more than one device already open. If
you just use dev.off() you will close the last device that was opened,'

#pdf(file=tmp7) # Move before the scatterplot3d call

#dev.print(file=tmp7, device=pdf, width=600)

# you should also make sure you created a valid file string. We cannot
check since you have not offered a reproducible example.

#dev.off(2)
}

David Winsemius, MD
West Hartford, CT