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