ggplot(aes(x=as.factor(year),y=residual),data=op)+geom_violin()+
    geom_boxplot(width=0.1,fill="gray")+xlab("year")