event_times が vector であり、 がこれらのインデックスと順番に並ぶevent_times <- c(5,3,6,8)
ことが予想される場合、グループを event_times へのインデックスとして使用できます。これは、パネル引数の x および y のインデックスになります。groups
xyplot(Y~time,groups=subject_identifier,data=cdat,
panel=function(x,y,groups,...){
panel.xyplot(x,y,groups=groups,...)
panel.segments(x0=x[(grp-1)*10+event_times[grp]],
x1=x[(grp-1)*10+event_times[grp]],
y0=y[(grp-1)*10+event_times[grp]]-.2,
y1=y[(grp-1)*10+event_times[grp]]+.2, groups=groups,...)},
type=c('l'))
したがって、最初のグループは time=5 でティックされ、2 番目のグループは time=3 でティックされます。(group-1)*10
を使用したよりクリーンなアプローチがあるかもしれませんが、相殺する必要がありますpanel.superpose
。もちろん、ティックの他のオプションも可能ですが、問題のプレゼンテーションはかなり曖昧です。
