PDF output problem - R-help

Tue, Aug 27, 2002 9:09 AM #

Hello,
	I'm quite new to R, but I've already stepped into this problem:
I open a PDF device with:
pdf("Name-%d.pdf", width=10, height=10, onefile=FALSE)

And draw 4 histograms in a row, expecting the pdf device to
automatically number them from 1 to 4. What I get back is
only 2 images with names "Name-1.pdf" and "Name-2.pdf" that
contain the two last histograms I've plotted. Is it a known
problem? Or am I doing something wrong here?

I attach the script to this message in the hope somebody knows
the cause of this strange behaviour. Thanks in advance for
any help.

	Andre.

Andre Rabello DOS ANJOS, M.Sc.
Signal Processing, Data Analysis and Computing
Office: 32-2-A06, Tel: (+ 41 22) 767 5022
Fax: (+ 41 22) 767 8420
CERN - EP/HC Division
CH-1211 Geneve 23 - Suisse/Switzerland
http://cern.ch/rabello
-------------- next part --------------
# Hello emacs (ess), this is -*- R -*-
# $Id$
# Andre DOS ANJOS <Andre.dos.Anjos at cern.ch>

# This script will run R, producing histograms
# generated with data from nat.

# Calculates the standard deviation of a vector
std <- function (x) {
  return(sqrt(sum((x-mean(x))^2)/(length(x)-1)))
}

# This function will produce a fancy histogram in R. The histogram,
# besides having the normal look, will display the total number of
# entries, the average and standard deviation (aka, sigma). 
nathist <- function (d, breaks = "Sturges",
                     main = "NetLogger Histogram",
                     xlab = expression(paste("Time (", mu, "s)")),
                     ylab = "Counts") {
  r <- hist(d, breaks=breaks, main=main, xlab=xlab, ylab=ylab,
            border='black', col='yellow')
  y <- max(r$counts) * 0.9
  x <- r$breaks[length(r$breaks)-1] * 0.7
  text(x, y, paste("Entries: ", length(d)), pos=4, col='blue3')
  text(x, 0.97*y, paste("Average: ",mean(d)), pos=4, col='blue3')
  text(x, 0.94*y, paste("Deviation: ",std(d)), pos=4, col='blue3')
  rect(0.99*x, 0.92*y, 0.99*r$breaks[length(r$breaks)], y*1.03,
       border='black')
  ##text(r$mids, r$counts*1.005, r$counts, pos=3)
  ## return(r)
}

#Set the column names
colname <- c("ROW.NAMES", "PUT.L1RESULT", "GET.L1RESULT", "START.PROCESSING",
             "GAUDI.START.ALGO.PROCESSING", "ROSDATA.REQUESTED.1",
             "ROSDATA.RECEIVED.1", "ROSDATA.REQUESTED.2",
             "ROSDATA.RECEIVED.2", "GAUDI.END.ALGO.PROCESSING",
             "END.PROCESSING", "SEND.L2DECISION", "DELETE.L1MESG")

#Read-in the data table
my <- read.table("file.dat", col.names=colname, row.names=1)

# Data Transfer Overhead
datareq = my[,"ROSDATA.RECEIVED.1"] - my[,"ROSDATA.REQUESTED.1"] +
  my[,"ROSDATA.RECEIVED.2"] - my[,"ROSDATA.REQUESTED.2"]

pdf(file="nat-rose.pdf", heigh=10, width=10, onefile=TRUE)
nathist(datareq, main="L2PU/ROSE Data Transfer Overhead", breaks=25)

# Gaudi Timings
gaudi = my[,"GAUDI.END.ALGO.PROCESSING"] -
  my[,"GAUDI.START.ALGO.PROCESSING"] - datareq
pdf(file="nat-pesa.pdf", heigh=10, width=10, onefile=TRUE)
nathist(gaudi, main="PESA Timings", breaks=50)

# L2PU Overhead
l2pu = my[,"DELETE.L1MESG"] - my[,"GET.L1RESULT"] - datareq - gaudi
pdf(file="nat-l2pu.pdf", heigh=10, width=10, onefile=TRUE)
nathist(l2pu, main="L2PU Overhead", breaks=50)

# Queue Timing
queue = my[,"GET.L1RESULT"] - my[,"PUT.L1RESULT"]
pdf(file="nat-queue.pdf", heigh=10, width=10, onefile=TRUE)
nathist(queue, main="Queue Timings", breaks=50)

graphics.off()
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 232 bytes
Desc: not available
Url : https://stat.ethz.ch/pipermail/r-help/attachments/20020827/4d6cf11b/attachment.bin