[R] lattice --- different properties of lines corresponding to type=c("l", "a") respectively

[Deepayan Sarkar]

On Fri, Nov 27, 2009 at 1:13 AM, Kjetil Halvorsen
<kjetilbrinchmannhalvorsen at gmail.com> wrote:
> I think the subject says it all. I want to make a simple lattice plot,
> using xyplot with the
> argument     type=c("l","a").
> The problem then is that in the resulting plot it is
> difficult/impossible to see which plot corresponds to the average
> and which to the individual profiles. I triedthings like extra
> arguments  lwd=c(1,3)   or   col=c("blue","red")
> hoping this would be interpreteded parallely to the type=  argument,
> but no.  Like:
>
> xyplot(response ~ time|group, repa0, groups=~Participant, type=c("b",
> "a"), lwd=c(1, 3),
>       ylim=c(0, 10))

Are you looking for per-group average or overall average? All your
attempts will do per-group average, which is rather pointless here.
For a global average you need to call panel.average with all the data,
not through panel.superpose:

mypanel <- function(x, y, ..., groups, type){
    panel.xyplot(x, y, ..., groups = groups, type=type)
    panel.average(x, y, horizontal = FALSE,
                  col = "black", lwd = 3)
}

xyplot(response ~ time|group, repa0, groups=~Participant, type="l",
       ylim=c(0, 10), panel=mypanel)


-Deepayan

>
> and many other variants ...
>
> Then I vent for writing my own panel functions:
> mypanel <- function(x, y, ..., type){
>    panel.average(x, y, ...,  horizontal=FALSE)
>    panel.xyplot(x, y, ..., type=type)
>  }
> xyplot(response ~ time|group, repa0, groups=~Participant, type="l",
> ylim=c(0, 10),
>       panel=panel.superpose, panel.groups=mypanel)
>
> (which doesn't work)
> ???
> Kjetil Halvorsen
> To recreate the data:
>
> repa0 <-
> structure(list(group = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L,
> 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
> 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
> 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
> 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L,
> 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
> 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
> 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
> 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
> 2L, 2L, 2L, 2L, 2L), .Label = c("1", "2"), class = "factor"),
>    Participant = c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L,
>    11L, 12L, 13L, 14L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L,
>    11L, 12L, 13L, 14L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L,
>    11L, 12L, 13L, 14L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L,
>    11L, 12L, 13L, 14L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L,
>    11L, 12L, 13L, 14L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L,
>    11L, 12L, 13L, 14L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L,
>    11L, 12L, 13L, 14L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L,
>    11L, 12L, 13L, 14L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L,
>    11L, 12L, 13L, 14L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L,
>    11L, 12L, 13L, 14L), time = c(1L, 1L, 1L, 1L, 1L, 1L, 1L,
>    1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
>    2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
>    3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
>    4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
>    5L, 5L, 5L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
>    1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
>    2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
>    4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 5L,
>    5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L), response = c(2L,
>    10L, 8L, 4L, 0L, 4L, 10L, 2L, 3L, 4L, 1L, 2L, 3L, 3L, 5L,
>    2L, 10L, 3L, 0L, 3L, 10L, 1L, 7L, 3L, 2L, 10L, 8L, 5L, 8L,
>    0L, 10L, 2L, 0L, 2L, 10L, 10L, 3L, 2L, 10L, 4L, 6L, 0L, 5L,
>    0L, 6L, 2L, 0L, 2L, 10L, 5L, 3L, 2L, 6L, 4L, 5L, 0L, 5L,
>    0L, 5L, 0L, 0L, 0L, 10L, 0L, 0L, 2L, 5L, 0L, 5L, 0L, 3L,
>    8L, 6L, 1L, 3L, 0L, 5L, 4L, 10L, 0L, 3L, 3L, 6L, 5L, 2L,
>    0L, 2L, 0L, 4L, 0L, 5L, 4L, 5L, 1L, 0L, 2L, 0L, 5L, 3L, 1L,
>    8L, 2L, 4L, 0L, 0L, 0L, 10L, 4L, 3L, 0L, 0L, 10L, 3L, 3L,
>    6L, 1L, 4L, 0L, 0L, 0L, 10L, 3L, 3L, 2L, 2L, 10L, 2L, 3L,
>    5L, 0L, 1L, 0L, 0L, 0L, 10L, 2L, 3L, 2L, 0L, 6L)), .Names = c("group",
> "Participant", "time", "response"), row.names = c("1.1", "2.1",
> "3.1", "4.1", "5.1", "6.1", "7.1", "8.1", "9.1", "10.1", "11.1",
> "12.1", "13.1", "14.1", "1.2", "2.2", "3.2", "4.2", "5.2", "6.2",
> "7.2", "8.2", "9.2", "10.2", "11.2", "12.2", "13.2", "14.2",
> "1.3", "2.3", "3.3", "4.3", "5.3", "6.3", "7.3", "8.3", "9.3",
> "10.3", "11.3", "12.3", "13.3", "14.3", "1.4", "2.4", "3.4",
> "4.4", "5.4", "6.4", "7.4", "8.4", "9.4", "10.4", "11.4", "12.4",
> "13.4", "14.4", "1.5", "2.5", "3.5", "4.5", "5.5", "6.5", "7.5",
> "8.5", "9.5", "10.5", "11.5", "12.5", "13.5", "14.5", "1.11",
> "2.11", "3.11", "4.11", "5.11", "6.11", "7.11", "8.11", "9.11",
> "10.11", "11.11", "12.11", "13.11", "14.11", "1.21", "2.21",
> "3.21", "4.21", "5.21", "6.21", "7.21", "8.21", "9.21", "10.21",
> "11.21", "12.21", "13.21", "14.21", "1.31", "2.31", "3.31", "4.31",
> "5.31", "6.31", "7.31", "8.31", "9.31", "10.31", "11.31", "12.31",
> "13.31", "14.31", "1.41", "2.41", "3.41", "4.41", "5.41", "6.41",
> "7.41", "8.41", "9.41", "10.41", "11.41", "12.41", "13.41", "14.41",
> "1.51", "2.51", "3.51", "4.51", "5.51", "6.51", "7.51", "8.51",
> "9.51", "10.51", "11.51", "12.51", "13.51", "14.51"), class = "data.frame")
>
> ______________________________________________
> R-help at 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.
>