Как добавить среднюю строку к xyplot в R?

Question

Привет Я пытаюсь добавить среднюю линию к решетке с помощью графика XY следующего кода ... который не работает ...

library(lattice) 
temp <- data.frame(c(1,1,1,1,1,1,2,2), 
     c(1,1,2,2,1,1,1,1), 
     c(1,1,1,2,1,2,1,2), 
     c(0,0,0,150,150,0,0,0), 
     c(3,10,4,15,14,2,4,5)) 
names(temp) <- c("instNo","NoOfAC","NoReleaseLevels","trExt","maxDelay") 

tempavg <- aggregate(temp[5], by=list(temp$NoOfAC,temp$NoReleaseLevels 
            ,temp$trExt),FUN = mean) 

names(tempavg) <- c("NoOfAC","NoReleaseLevels","trExt","meanDelay") 
xyplot(maxDelay ~ NoOfAC |trExt + NoReleaseLevels, 
      data = temp, 
      panel = function(x, y, tempavg) { 
       panel.xyplot(x, y) 
       panel.abline(tempavg$meanDelay ~ y |trExt + NoReleaseLevels) 
      }, 
xlab = "Number Of Aircraft in the Instance", 
ylab = "Total Delay (minutes)") 

Answer 1

Основная проблема с подходом, где это говорит, что panel.abline() принимает формулу в качестве первого аргумента? (Это не так!) Что такое ?panel.abline, так это то, что он принимает аргументы h и v (в частности, для случая, который вы имеете в виду), и вы хотите аргумент h.

Заметим также, что решетка вызывает функции панели с аргументами x и y, содержащих x и y данные для каждой панели нарисованной. Вам не нужно выполнять подмножество и агрегирование извне для этой работы; это делается для вас. Следовательно, вы должны думать в терминах x и y при рисовании с функциями панели.

В этом случае вы хотите panel.abline(h = mean(y), ...) в своей функции панели. ... позволяет передавать дополнительные аргументы, переданные функции панели, для распространения до вызовов функций внутри функции панели. Следовательно, ваша анонимная функция панели также должна принимать аргумент .... Это показано ниже.

## code as factor 
temp <- transform(temp, NoOfAC = factor(NoOfAC)) 

## plot 
xyplot(maxDelay ~ NoOfAC |trExt + NoReleaseLevels, 
      data = temp, 
      panel = function(x, y, ...) { 
       panel.xyplot(x, y, ...) 
       panel.linejoin(x, y, horizontal = FALSE, ...) 
      }, 
     xlab = "Number Of Aircraft in the Instance", 
     ylab = "Total Delay (minutes)")