如何在R中的x轴上使用两个分类变量制作散点图

我正在尝试在x轴上使用两个分类变量在R中创建散点图。对于箱线图，我知道如何执行此操作（请参见下面的代码的第一部分），但是以某种方式，我无法使其用于散点图。我已经尝试了几件事，但是当我绘制点时，它们总是重叠并且不再显示第二个分类变量。抖动也不起作用，因为我希望我的类别聚类并且不要随机分散它们。有谁知道如何做到这一点？在下面，您可以找到一些示例数据和我尝试过的一些图形，包括注释。第一张图给我的东西与我想要的东西相似，但随后是一个箱线图而不是散点图。第二张图给出了一个散点图（人为第二个分类变量创建数字），

为了使其更加复杂，我还想显示一条带有所有散点图的平均值线。与在R中使用ggplot2在具有均值段的分类散点图中执行的操作类似。如何添加呢？

感谢你的帮助！

time = c(rep('t1',12),rep('t2',12))
Origin =  c(rep('I1B',4),rep('I1C',4),rep('J4A',4),rep('I1B',4),rep('I1C',4),rep('J4A',4))
LB_FR = runif(24)

df = data.frame(time,Origin,LB_FR)

#does not work with geom_point
ggplot(df, aes(x = time, y = LB_FR, fill = Origin)) + geom_boxplot() + ggtitle('LB_FR')

#create df_2 with numbers instead of categories for Origin
df_2 = df
for (r in 1:nrow(df)){
  if (df$Origin[r] == 'I1B') df_2[r,'OriginNr'] = 1
  if (df$Origin[r] == 'I1C') df_2[r,'OriginNr'] = 2
  if (df$Origin[r] == 'J4A') df_2[r,'OriginNr'] = 3
}

# indices for time
t1 = df_2$time=="t1"
t2 = df_2$time=="t2"

plot(df_2$OriginNr,df$LB_FR, 
     xlim = c(0,4), ylim = c(0,1), bty = 'n',
     main = 'LB_FR', ylab = 'Fraction remaining', xlab = 'Origin', type = 'n')
points(df_2$OriginNr[t1],df_2$LB_FR[t1],col='red')
points(df_2$OriginNr[t2],df_2$LB_FR[t2],col='blue')
legend(0.1,0.9,legend=c('month 0-6','month 6-12'),pch=1,col=c('red','blue'),bty='n',cex=1.2)