completeFun <- function(data, desiredCols) { completeVec <- complete.cases(data[, desiredCols]) return(data[completeVec, ]) } conditions <- make_conditions(model.3way, "period") me<-marginal_effects(model.3way, "time:treatment", conditions = conditions) me$`time:treatment`$`effect1__`[me$`time:treatment`$period=="1pre" & me$`time:treatment`$time>10]<-NA me$`time:treatment`$`effect2__`[me$`time:treatment`$period=="1pre" & me$`time:treatment`$time>10]<-NA me$`time:treatment`$`estimate__`[me$`time:treatment`$period=="1pre" & me$`time:treatment`$time>10]<-NA me$`time:treatment`$`se__`[me$`time:treatment`$period=="1pre" & me$`time:treatment`$time>10]<-NA me$`time:treatment`$`lower__`[me$`time:treatment`$period=="1pre" & me$`time:treatment`$time>10]<-NA me$`time:treatment`$`upper__`[me$`time:treatment`$period=="1pre" & me$`time:treatment`$time>10]<-NA me$`time:treatment`$`effect1__`[me$`time:treatment`$period=="2post" & me$`time:treatment`$time<11]<-NA me$`time:treatment`$`effect2__`[me$`time:treatment`$period=="2post" & me$`time:treatment`$time<11]<-NA me$`time:treatment`$`estimate__`[me$`time:treatment`$period=="2post" & me$`time:treatment`$time<11]<-NA me$`time:treatment`$`se__`[me$`time:treatment`$period=="2post" & me$`time:treatment`$time<11]<-NA me$`time:treatment`$`lower__`[me$`time:treatment`$period=="2post" & me$`time:treatment`$time<11]<-NA me$`time:treatment`$`upper__`[me$`time:treatment`$period=="2post" & me$`time:treatment`$time<11]<-NA me$`time:treatment`<-completeFun(me$`time:treatment`, "effect1__") plot1<-plot(me) plot1 + facet_wrap("cond__", scales = "free")