Dear R users

I am trying to exporting plots from R to an external folder, or to the
working directory, but the resolution of plots (pdf file) largely reduced.
Any way I can get same quality as my original plots?? e.g. I tested the
plotting part using one example and obtained pretty good (/readable)
quality for each plot in the (4*4) multiple graph. But when I did the loop
and tried to export those plots out using 'dev.copy', the quality was not
same. I do need this loop function as there are 400 datasets so I can not
handle them manually.

Also I am using 'xtable' to write Latex code for my summary tables of the
data (again there are 400 datasets). I saved those Latex code in a list and
use 'sink()' to save them outside the loop function, which works ok. But as
my table is wide, so I found that the Latex code produced by xtable can not
fit my purpose well, e.g. I have set size to 'tiny' but the table still too
wide. Is there an alternative package more handy to do the job? Also are
there a good way to write some hundred tables from R to Latex for producing
pdf documents in an easier fashion?? i.e. Just realised that Latex does not
like compiling large amount tables/figures in one go!

Any advices/ideas are greatly appreciated!


Best wishes
HJ




Below is my code...
=====================
SumTab<-function(Data=SortedDataInList,StartDate="30/1/12",EndDate="31/05/12",StartTime="22:50:00",EndTime="23:00:00"){
  Start= chron(StartDate,StartTime , format=c(dates="d/m/y", times="h:m:s"))
  End= chron(EndDate,EndTime , format=c(dates="d/m/y", times="h:m:s"))
  deltat <- times("00:10:00")
  TT <- seq(Start,End, by = times("00:10:00"))
  TT1 = substr(TT, 2, 18)


  Data1=Data
  for (i in 1:length(Data1)){
    SumTab1[[i]]= matrix(NA, nrow=5, ncol=ncol(SortedDataInList[[i]])-4)
    SortedDataInList1[[i]]= matrix(NA, nrow=length(TT1),
ncol=ncol(Data1[[i]]))

SortedDataInList1[[i]]=Data1[[i]][match(as.character(TT1),as.character(Data1[[i]][,1])),]


    SumTab1[[i]][1,]<-apply(SortedDataInList1[[i]][,4:16],2,min)
    SumTab1[[i]][2,]<-apply(SortedDataInList1[[i]][,4:16],2,mean)
    SumTab1[[i]][3,]<-apply(SortedDataInList1[[i]][,4:16],2,median)
    SumTab1[[i]][4,]<-apply(SortedDataInList1[[i]][,4:16],2,sd)
    SumTab1[[i]][5,]<-apply(SortedDataInList1[[i]][,4:16],2,max)



colnames(SumTab1[[i]])=c("vOL1","VOL2","VOL3","CUR1","CUR2","CUR3","THD1","THD2","THD3","RPD","RPR","RAPD","RAPR")
    rownames(SumTab1[[i]])=c("Min","Mean","Standard
Deviation","Median","Max")


SumLax[[i]]<-xtable(SumTab1[[i]],label=as.character(StationsInDir[i]),caption=as.character(StationsInDir[i]))


    par(mfrow=c(4,4),oma=c(4,0,2,0))

plot(1:length(SortedDataInList1[[i]][,1]),SortedDataInList1[[i]][,4],xlim=c(1,length(SortedDataInList1[[i]][,1])),ylab="Voltage
1 (v)",xlab="Date/Time ind.(10 min. int.)",type="p",cex=.001)

plot(1:length(SortedDataInList1[[i]][,1]),SortedDataInList1[[i]][,5],xlim=c(1,length(SortedDataInList1[[i]][,1])),ylab="Voltage
2 (v)",xlab="Date/Time ind.(10 min. int.)",type="p",cex=.001)

plot(1:length(SortedDataInList1[[i]][,1]),SortedDataInList1[[i]][,6],xlim=c(1,length(SortedDataInList1[[i]][,1])),ylab="Voltage
3 (v)",xlab="Date/Time ind.(10 min. int.)",type="p",cex=.001)

plot(1:length(SortedDataInList1[[i]][,1]),SortedDataInList1[[i]][,7],xlim=c(1,length(SortedDataInList1[[i]][,1])),ylab="Current
1 (A)",xlab="Date/Time ind.(10 min. int.)",type="p",cex=.001)

plot(1:length(SortedDataInList1[[i]][,1]),SortedDataInList1[[i]][,8],xlim=c(1,length(SortedDataInList1[[i]][,1])),ylab="Current
2 (A)",xlab="Date/Time ind.(10 min. int.)",type="p",cex=.001)

plot(1:length(SortedDataInList1[[i]][,1]),SortedDataInList1[[i]][,9],xlim=c(1,length(SortedDataInList1[[i]][,1])),ylab="Current
3 (A)",xlab="Date/Time ind.(10 min. int.)",type="p",cex=.001)

plot(1:length(SortedDataInList1[[i]][,1]),SortedDataInList1[[i]][,10],xlim=c(1,length(SortedDataInList1[[i]][,1])),ylab="Total
harmonic distortion 1 (%)",xlab="Date/Time ind.(10 min.
int.)",type="p",cex=.001)

plot(1:length(SortedDataInList1[[i]][,1]),SortedDataInList1[[i]][,11],xlim=c(1,length(SortedDataInList1[[i]][,1])),ylab="Total
harmonic distortion 2 (%)",xlab="Date/Time ind.(10 min.
int.)",type="p",cex=.001)

plot(1:length(SortedDataInList1[[i]][,1]),SortedDataInList1[[i]][,12],xlim=c(1,length(SortedDataInList1[[i]][,1])),ylab="Total
harmonic distortion 3 (%)",xlab="Date/Time ind.(10 min.
int.)",type="p",cex=.001)

plot(1:length(SortedDataInList1[[i]][,1]),SortedDataInList1[[i]][,13],xlim=c(1,length(SortedDataInList1[[i]][,1])),ylab="Real
power delivered  (mw)",xlab="Date/Time ind.(10 min.
int.)",type="p",cex=.001)

plot(1:length(SortedDataInList1[[i]][,1]),SortedDataInList1[[i]][,14],xlim=c(1,length(SortedDataInList1[[i]][,1])),ylab="Real
power received  (mw)",xlab="Date/Time ind.(10 min. int.)",type="p",cex=.001)

plot(1:length(SortedDataInList1[[i]][,1]),SortedDataInList1[[i]][,15],xlim=c(1,length(SortedDataInList1[[i]][,1])),ylab="Reactive
power delivered (MVAr)",xlab="Date/Time ind.(10 min.
int.)",type="p",cex=.001)

plot(1:length(SortedDataInList1[[i]][,1]),SortedDataInList1[[i]][,16],xlim=c(1,length(SortedDataInList1[[i]][,1])),ylab="Reactive
power received (MVAr)",xlab="Date/Time ind.(10 min.
int.)",type="p",cex=.001)
    title(main=StationsInDir[i],outer=TRUE,cex.main=2.0)
    mtext(side=1,Start, outer = TRUE, line=0,cex = 0.6)
    mtext(side=1,End, outer = TRUE, line=2,cex = 0.6)


    #pdf(paste(as.character(StationsInDir[i]),".pdf"))
    dev.copy(pdf,file=paste(as.character(StationsInDir[i]),".pdf"))
    dev.off()

  }
}

	[[alternative HTML version deleted]]

