将 emmeans 输出范围指定为单个因子水平的范围

Question

我正在尝试使用 emmeans 的输出创建一个图形，为一个因子的 5 个水平绘制线条。我希望每个功能区的范围与出现该级别的 x 轴上的数据范围相对应，而不是在整个 x 轴上。即，某些因素仅在 x 轴的特定范围内具有数据，我不想推断超出这些范围。

在整个范围内推断的当前代码是：

newdata=emmeans(model, ~x|factor, at=list(factor=levels(data$factor), x=seq(min(data$x), max(data$x), len=100)), type='response') %>% as.data.frame

figure=ggplot(data, aes(y=y, x=x, color=factor, fill=factor))+
  geom_ribbon(data=newdata, aes(x=x, y=response,ymin=lower.CL, ymax=upper.CL), alpha=0.3, colour = NA)+
   geom_line(data=newdata, aes(x=x, y=response))
   
figure

Answer 1

我后来找到了一个庞大的解决方法：

#Build dataframes with max and min for each factor
factorvariable.1 <- c("factorvariable.1")
data.factorvariable.1=filter(data, factor %in% factorvariable.1)
factorvariable.1.range=range(data.factorvariable.1$x)%>% as.data.frame
factorvariable.1.range$factor=factorvariable.1
factorvariable.1.range$min.max=c('min','max')

factorvariable.2 <- c("factorvariable.2")
data.factorvariable.2=filter(data, factor %in% factorvariable.2)
factorvariable.2.range=range(data.factorvariable.2$x)%>% as.data.frame
factorvariable.2.range$factor=factorvariable.2
factorvariable.2.range$min.max=c('min','max')

Range=rbind(factorvariable.1.range,factorvariable.2.range)

Range <- spread(Range, min.max, .)

#filter emmeans data by max and min values
newdata=emmeans(model, ~x|factor, at=list(factor=levels(data$factor), x=seq(min(data$x), max(data$x), len=100)), type='response') %>% as.data.frame
newdata=merge(newdata, Range, by="factor")
newdata= newdata%>%filter(x>min)
newdata= newdata%>%filter(x<max)
newdata