マルチパネルプロットのラティスパネル.abline

Question

以下のコードを使用して、複数のパネルを含む格子xyplotを作成しようとしています。問題は、すべてのパネルにすべての中央線を追加していることです（つまり、アブラインはSubjectでグループ化されていません）。サブジェクトごとのメジアンラインをサブセット化する方法はありますか？

以下は、PKPDdatasetsパッケージのTheophデータセットを使用した再現可能な例です。

# devtools::install_github("dpastoor/PKPDdatasets") 
library(PKPDdatasets) 
library(lattice) 
library(latticeExtra) 
library(dplyr) 

Theoph_Data <- Theoph 

Theoph_Data1 <- Theoph_Data %>% group_by(Subject) %>% mutate(Median= median(Time)) 

print(Theoph_Data1) 

Theoph_PK <- xyplot(conc ~ Time| Subject, data=Theoph_Data1, 
        panel = function(x,y,...) { 
        panel.xyplot(x,y,...) 
        panel.abline(v=Theoph_Data1$Median) 
        }) 
Theoph_PK 

Answer 1

あなたの例では、最高ではないですが、すべての個人が同じ期間の中央値を持っているので、基本的な考え方はこれです：

Theoph_Data <- Theoph 
xyplot(conc ~ Time| Subject, data=Theoph_Data1, 
    panel = function(x,y,subscripts,...) { 
     panel.xyplot(x,y, subscripts=subscripts,...) 
     panel.abline(v=median(Theoph_Data1$Time[subscripts])) 
    } 
) 

あなたのパネル機能でsubscripts=パラメータを要求することができます。これにより、各パネルにプロットされた観測値のインデックスが得られます。それらのサブスクリプトを使用して元のデータをサブセット化することができます。ここでは関数内の中央値を計算します（単一の値を持つ単一の線をプロットすることを確実にするため）。