R help - Jan 2014 - How do I perform conditional annotation of lattice panel plots?

My apologies for asking this question that may have been asked before. I am
trying to plot activity dependent on time conditioned by the subject. Code for
sample data below.

So I have something like this

xyplot(Activity~Time|Subject).
This works fine, but now I want to show where on these activity curves Event A
occurs. This is to explore the relationship between A and activity.

I tried this:

xyplot(EventA+Activity~Time|Subject)

but then the numerical ordering gets out of order.
Ideally Event A could be denoted by a vertical line cutting through the curve at
the day Event A occurred, but some other way to denote it would be great!

Here is my sample data as R code:

Subject<-c(1,1,1,1,1,1,2,2,2,2,2,2,3,3,3,3,3,3,4,4,4,4,4,4)
Day<-c(1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,4,5,6)
Activity<-c(2,3,4,3,7,4,5,8,2,8,4,6,2,5,3,8,9,5,6,3,4,5,6,7)
EventA<-c("Yes",NA,"Yes",NA,NA,NA,"Yes",NA,NA,NA,NA,NA,
          "Yes",NA,NA,NA,NA,NA,NA,NA,"Yes",NA,NA,NA)
Data<-data.frame(cbind(Subject,Day,Activity,EventA))

  Subject Day Activity EventA
1        1   1        2    Yes
2        1   2        3   <NA>
3        1   3        4    Yes
4        1   4        3   <NA>
5        1   5        7   <NA>
6        1   6        4   <NA>
7        2   1        5    Yes
8        2   2        8   <NA>
9        2   3        2   <NA>
10       2   4        8   <NA>
11       2   5        4   <NA>
12       2   6        6   <NA>
13       3   1        2    Yes
14       3   2        5   <NA>
15       3   3        3   <NA>
16       3   4        8   <NA>
17       3   5        9   <NA>
18       3   6        5   <NA>
19       4   1        6   <NA>
20       4   2        3   <NA>
21       4   3        4    Yes
22       4   4        5   <NA>
23       4   5        6   <NA>
24       4   6        7   <NA>
Any ideas on how to do this would be appreciated!

Thanks very much in advance!


This email message may contain legally privileged and/or confidential
information.  If you are not the intended recipient(s), or the employee or agent
responsible for the delivery of this message to the intended recipient(s), you
are hereby notified that any disclosure, copying, distribution, or use of this
email message is prohibited.  If you have received this message in error, please
notify the sender immediately by e-mail and delete this email message from your
computer. Thank you.
	[[alternative HTML version deleted]]

Chen, George <George.Chen <at> roswellpark.org>
writes:
> My apologies for asking this question that may have 
been asked before. I am trying to plot activity
> dependent on time conditioned by the subject. 
Code for sample data below.
> So I have something like this
> xyplot(Activity~Time|Subject).
> This works fine, but now I want to show 
where on these activity curves Event A occurs. 
This is to explore the
> relationship between A and activity.
> 
> I tried this:
> 
> xyplot(EventA+Activity~Time|Subject)
> 
> but then the numerical ordering gets out of order.
> Ideally Event A could be denoted by a vertical line 
cutting through the curve at the day Event A occurred,
but
> some other way to denote it would be great!
> 
-----snip -------------
> Any ideas on how to do this would be appreciated!
> 
> Thanks very much in advance!
> 
> This email message may contain legally privileged 
and/or confidential information.  If you are not the
> intended recipient(s), or the employee or agent 
responsible for the delivery of this message to the
> intended recipient(s), you are hereby notified 
that any disclosure, copying, distribution, or use of
> this email message is prohibited.  If you have 
received this message in error, please notify the sender
> immediately by e-mail and delete this email 
message from your computer. Thank you.
> 	[[alternative HTML version deleted]]
You used Time in your formula but Day in 
your data so I changed the formula to Day
assuming that's what you meant

Perhaps, not the most elegant solution, but 
how about something like this

library(lattice)
Subject<-c(1,1,1,1,1,1,2,2,2,2,2,2,3,3,3,3,3,3,4,4,4,4,4,4)
Day<-c(1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,4,5,6)
Activity<-c(2,3,4,3,7,4,5,8,2,8,4,6,2,5,3,8,9,5,6,3,4,5,6,7)
EventA<-c("Yes",NA,"Yes",NA,NA,NA,"Yes",NA,NA,NA,NA,NA,
          "Yes",NA,NA,NA,NA,NA,NA,NA,"Yes",NA,NA,NA)
Data<-data.frame(cbind(Subject,Day,Activity,EventA))

xyplot(Activity~ Day | Subject)
xyplot(Activity ~ Day | Subject, data = Data,
	subscripts = TRUE,
	panel = function(x, y, subscripts, ...){
		panel.xyplot(x, y)
		wh <- Data[subscripts, ]
		panel.abline(v = wh$Day[!is.na(wh$EventA)])
	})


-- 
Kenneth Knoblauch
Inserm U846
Stem-cell and Brain Research Institute
Department of Integrative Neurosciences
18 avenue du Doyen L?pine
69500 Bron
France
tel: +33 (0)4 72 91 34 77
fax: +33 (0)4 72 91 34 61
portable: +33 (0)6 84 10 64 10
http://www.sbri.fr/members/kenneth-knoblauch.html