使用多个因子和2个分组变量来均衡箱形图的宽度，但是有些因子缺少其中一个分组

Question

我正在展示下面制作的箱形图的图片以及用于生成它的代码。

library(ggplot2)
ggplot(subset(Cyt_merged, !is.na(Hu.PDGF)), 
       aes(timepoint_coded, Hu.PDGF, fill=Group)) +  
  stat_boxplot(geom ='errorbar', position=position_dodge(width=.88)) +  
  geom_boxplot(width=1, outlier.size=.5, position=position_dodge(width=.88)) +
  ggtitle("Comparison of PDGF levels between groups over time") +  
  xlab("Timepoint") + ylab("PDGF level") +
  scale_fill_discrete(name="Group") + theme_bw() +
  theme(
    plot.title=element_text(family="Times", size=14.5, face="bold"), 
    axis.title.x=element_text(family="Times", size=12, face="bold"), 
    axis.title.y=element_text(family="Times", size=12, face="bold"), 
    legend.key=element_rect(colour=NA), legend.key.size=unit(3, "line") 
  )

如您所见，我实际上没有72小时或2周时间点的控件数据。我想使这些时间点上的患者的箱形图的大小与其他时间点上各组的箱形图的大小相同，并且我希望这些图以因子轴为中心。有关如何执行此操作的任何想法？