Hi
Dataframe fdf contains
bin rate overlay
1 1 90 Assign First/cc _from_SN_53 RNC_20_to_SN_50 RNC_21_Success
Rate
2 2 93 Assign First/cc _from_SN_53 RNC_20_to_SN_50 RNC_21_Success
Rate
3 1 90 Assign First/cc _from_SN_50 RNC_21_to_SN_53 RNC_20_Success
Rate
4 2 94 Assign First/cc _from_SN_50 RNC_21_to_SN_53 RNC_20_Success
Rate
time
1 (04/01/08 16:02:30)
2 (04/01/08 16:07:30)
3 (04/01/08 16:02:30)
4 (04/01/08 16:07:30)
And then I write the following lines of code to plot a xygraph on a pdf
file
n <- length(unique(fdf$time))
x <-sort(unique(fdf$time))
y<-UTCsecs2chron(x)
xscales <- computeTimeScales(y)
yscales<-NULL
scales<-c(xscales,yscales)
ylab<-"session transfer rate"
xlab<-"time"
lgnd.txt<-levels(fdf$overlay)
celnet.trellis.device(device="pdf", file="ravi_st.pdf", width=10,height
= 10)
main<-"this is the first report"
formd="rate ~ as.numeric(time)"
print(
xyplot(formula(formd),
data = fdf, groups = overlay,
type = "b", as.table = TRUE, cex=0.20,
#subset = ok, commented to show breaks in graph
main = main, xlab = "Time", ylab = ylab,
scales = scales,
key = simpleKey(text = lgnd.txt, cex = 3/4,
points = FALSE, lines = TRUE),
page = function(n) annotate(opts$ann),
sub = "",
layout = c(1,1)
)
)
"computeTimeScales" <-
function(z, rot=0) ## z is a chron object
{browser()
## how many days do the data span?
n.days <- ceiling(diff(range(z, na.rm = TRUE)))
## compute x ticks, pretty labels
x <- as.numeric(z)
r <- range(x, na.rm = TRUE)
at.x <- seq(from = r[1], to = r[2], length = 7)
two_digits <- function(x){
x <- paste("0", x, sep="")
substring(x,nchar(x)-1)
}
## heuristics: use hh:mm when range of data falls within one day,
## otherwise use MM-DD hh:mm.
at.z <- chron(at.x)
hh <- paste(two_digits(hours(at.z)), two_digits(minutes(at.z)),
sep=":")
if(0){ ## hh:ss # this is commented as for more than one day
from/to time no date was printed
at.lbls <- hh
}
else { ## MM-DD\nhh:ss
m <- month.day.year(at.z)
dd <- paste(two_digits(m$month), two_digits(m$day), sep="-")
at.lbls <- paste(dd, hh, sep="\n")
}
list(x=list(at = at.x, rot=rot, labels = at.lbls))
}
"celnet.theme" <-
function()
{
celnet.theme <- canonical.theme("pdf", color = TRUE)
if(tolower(.Device) =="png")
celnet.theme$background$col <- "white"
else
canonical.theme("pdf", color = TRUE)
celnet.theme$background$col <- "transparent"
celnet.theme
}
"celnet.trellis.device" <-
function(device, file, width = NULL, height = NULL,...)
{
## analogous to trellis.device() -- this is just tailored to Celnet
if(is.null(width))
width <- 8
if(is.null(height))
height <- 6
if(tolower(device)=="png"){
if(missing(width)) width <- 72 * width
if(missing(height)) height <- 72 * height
}
if(tolower(device)=="x11")
trellis.device(device, width=width, height=height,
theme=celnet.theme())
if(tolower(device)=="ps" || tolower(device)=="postscript")
trellis.device("postscript", file=file, color = TRUE,
width=width, height=height, theme=celnet.theme())
else
trellis.device(device, file=file,
width=width, height=height, theme=celnet.theme())
}
"annotate" <-
function(str, ...)
## print a metadata message at the top-bottom of current trellis
display
{ else if(n==1){
x <- 0.95
hjust <- "right"
} else if(n==2){
x <- c(0.05, 0.95)
hjust <- c("left", "right")
} else if(n==3){
x <- c(0.05, 0.50, 0.95)
hjust <- c("left", "center", "right")
}
x <- unit(x, "npc")
y <- unit(1.5, "lines") ## at 0.5, 1.5, ... lines
for(i in seq(along = x)){
grid.text(label = str[i],
x = x[i], y = y, just = c(hjust[i], "bottom"), ...)
}
invisible(str)
}
n <- length(str) ## one line per string
if(n==0)
return(invisible(str))
}
But I did nt get any desired o/p on a file ravi_st.pdf
Can any one tell how can I get the desired o/p
[[alternative HTML version deleted]]
______________________________________________
R-help@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.
