How do I perform conditional annotation of lattice

Not sure if this adds much to Ken Knoblauch previous suggestion. But:



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)

require(loa)

panel.temp = function(x=x, y=y, z=z, ...){
                         panel.xyplot(x, y, z =1, ...)
                         panel.abline(v=x[!is.na(z)], col="red")
}
loaPlot(EventA ~ Day*Activity|Subject, panel=panel.temp, type="b")



The loa package is still a little messy, not yet on CRAN and the help is only half written but you are welcome to have a play:

https://sites.google.com/site/karlropkins/rpackages/loa



The basic plot is z ~ x * y, where z is information to be displayed at (x,y) and the package includes a very inelegant alternative to subscripting.



Regards,



Karl Ropkins

ITS Leeds

------------------------------





Date: Tue, 7 Jan 2014 15:50:39 +0000
From: Ken Knoblauch <ken.knoblauch at inserm.fr<mailto:ken.knoblauch at inserm.fr>>
To: <r-<mailto:r-help at stat.math.ethz.ch>help at stat.math.ethz.ch>
Subject: Re: [R] How do I perform conditional annotation of lattice
        panel   plots?
Message-ID: <loom.20140107T164813-513 at post.gmane.org<mailto:loom.20140107T164813-513 at post.gmane.org>>
Content-Type: text/plain; charset="utf-8"

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