черчения правильные значения для GAM модели в г

Question

У меня есть следующий пример (который находится в том же формате, моим данным):

set.seed(1) 
RandData <- runif(8760*2) 
Locations <- rep(c('UK','France'),each=8760) 

Date = seq(from=as.POSIXct("1991-01-01 00:00"), 
      to=as.POSIXct("1991-12-31 23:00"), length=8760) 

Final <- data.frame(Loc = Locations, 
        Doy = as.numeric(format(Date,format = "%j")), 
        Tod = as.numeric(format(Date,format = "%H")), 
        Temp = RandData) 
require(mgcv) 
mod1 <- gam(Temp ~ Loc + s(Doy) + s(Doy,by = Loc) + 
    s(Tod) + s(Tod,by = Loc),data = Final) 
plot(mod1,pages = 1, scale = 0) 

Из этой модели, исход которого отложена во-первых, показывает изменение средней температуры как функцию дня года, а затем величину, при которой температура в каждом месте изменяется от этого значения. То же самое показано на Время дня.

Проблема заключается в том, что yaxis не показывают правильный диапазон значений, т.е. первый график не показывает среднюю температуру для этих местоположений. Как можно изменить графики, чтобы показать (1) среднюю температуру для обоих мест и (2) температуру в каждом месте по сравнению со средним значением, а не какое-то произвольное число.

Если я не был ясен в выражении своих намерений, дайте мне знать, и я постараюсь предоставить лучший пример.

Answer 1

Графики модели гаммы отображают эффект самой гладкой. Если вам нужен комбинированный эффект нескольких частей модели, тогда создайте data.frame с соответствующими значениями для каждого значения, рассчитайте предсказание для этих значений с помощью прогноза (model, newdata) и создайте нужный сюжет.

dataset <- expand.grid(Loc = unique(Final$Loc), 
    Doy = pretty(Final$Doy), Tod = pretty(Final$Tod)) 
dataset$Temp <- predict(mod1, newdata = dataset) 
library(ggplot2) 
ggplot(dataset, aes(x = Doy, y = Temp, linetype = Loc)) + geom_line() + 
    facet_wrap(~Tod)