You could consider doing this directly with lattice. See:
demo(intervals)
On 5/30/07, Nitin Jain <[EMAIL PROTECTED]> wrote:
> Hello,
>
> I would like to get the scales of y-axes dependent only on the data points in
> a particular panel. Have attached a test example below.
> When using 'relation="free"', it does not make the scales 'free', however
> when using 'relation="sliced"', I get a warning "Explicitly specified limits
> ignored in: limitsFromLimitlist(have.lim = have.ylim, lim = ylim, relation =
> y.relation," (although in this particular case, I get the desired result, but
> in my real data, I do not get the free y-scale for each panel). Can you
> please let me know what should be correct syntax?
>
> Thanks.
> -Nitin
>
> library(Hmisc)
>
> test1 <- data.frame(
> y = c(rnorm(33), rnorm(33, mean=10),
> rnorm(34, mean=100)),
> x = 1:100,
> f = factor(c(rep("a", 33), rep("b", 33), rep("c", 34))),
> g = factor(sample(LETTERS[1:2], size=100, replace=TRUE))
> )
>
>
> CI <- rnorm(100)
> lb <- test1$y - CI
> ub <- test1$y + CI
>
>
> xYplot(Cbind(y, lb,ub )~x|f,
> groups=g,
> scales = list(relation="free"), ## Changing it to sliced gives warning
> data=test1)
>
>
>
>
>
>
>
> ____________________________________________________________________________________
> Don't pick lemons.
>
> ______________________________________________
> [email protected] 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.
>
______________________________________________
[email protected] 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.