plot(df$Year, df$percentage_of_RNAi,
     ylim=range(c(df$percentage_of_RNAi-sdev, df$percentage_of_RNAi+sdev)),
     pch=19, xlab="Years", ylab="Mean of [%] of RNAi research",
     main="Mean per year of RNAi research n as [%] of total research in the EU", col = ifelse(RNAi_research_EU$Year > 2011, "blue", "red") )