It worked for me..Here are the commands..
library(lattice)
testdata <- as.data.frame(t(structure(c(
1,2005,9.24,6.18,634,
2,2005,8.65,6.05,96,
8,2004,6.81,6.51,16,
9,2004,9.0,7.29,8,
10,2004,8.84,6.18,524,
11,2004,8.54,6.35,579,
12,2004,9.97,6.3,614,
12,2005,8.75,5.84,32,
), .Dim=c(5,8))))
colnames(testdata) <- c('month', 'year', 'mean','stdDev','miceCount')
testdata$month <- factor(testdata$month)
testdata$year <- factor(testdata$year)
testdata <- testdata[do.call("order", testdata), ]
png('lexstar_3241.png', width=600,
height=as.numeric(length(levels(testdata$year))*200), pointsize=8)
trellis.par.set(theme = col.whitebg())
with(testdata, print(barchart(as.numeric(mean) ~ month | year, data=testdata,
layout=c(1,length(levels(year))),
horizontal=FALSE,
scales=list(y=list(limits=c(1,max(as.numeric(mean))+max(as.numeric(stdDev))))),
main='Marble Burying - Level I',
xlab='Months',
ylab='Mean',
sd = as.numeric(as.character(stdDev)),
panel= function(x, y, ..., sd, subscripts) {
panel.barchart(x, y, ...);
sd <- sd[subscripts];
panel.segments(as.numeric(x), y - sd, as.numeric(x), y + sd,
col = 'red', lwd = 2);
}
)))
dev.off()
Thanks a lot Deepayan..
-Sandeep
-----Original Message-----
From: Sundar Dorai-Raj [mailto:[EMAIL PROTECTED]
Sent: Thursday, May 05, 2005 2:05 PM
To: Deepayan Sarkar
Cc: Ghosh, Sandeep; r-help@stat.math.ethz.ch
Subject: Re: [R] Re: Need some quick help with lattice - barchart
Deepayan Sarkar wrote on 5/5/2005 12:03 PM:
> On Thursday 05 May 2005 13:10, Ghosh, Sandeep wrote:
>
>>For the following code below, the x-axis ticks are 1,2,3,4,5,6,7 when I was
>>expection them to be 1,2,8,9,10,11,12. Please help me figure out where is
>>the mistake.
>
> [...]
>
>>colnames(testdata) <- c('month', 'year', 'mean','stdDev','miceCount')
>>testdata$month <- as.numeric(testdata$month)
>
> [...]
>
>>with(testdata, print(barchart(as.numeric(mean) ~ month | year,
>>data=testdata, layout=c(1,length(levels(year))),
>> horizontal=FALSE,
>
> [...]
>
> 'month' is numeric, so it's being coerced to be a shingle. Try using
> 'factor(month)' instead in the formula.
>
> Deepayan
>
Hi Deepayan,
That was my original thought too. But when I tried it I got:
Error in unit(x0, default.units, units.per.obs) :
'x' must be numeric
Changing horizontal to TRUE produces the plot, but I'm sure that's not
what Sandeep wants.
--sundar
______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
