scale_colour_manual(values=legend)
g
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, colour = "0-9y"), size=1)+
geom_abline(aes(intercept=1.65, slope=-1.31, colour = "9-12y"), size=1)+
geom_abline(aes(intercept=1.65, slope=-0.78, colour = "12-15y"), size=1)+
geom_abline(aes(intercept=1.65, slope=-0.23, colour = "15-16y"), size=1)+
labs(x="MMT", y="Mean Underproduction Factor", colour="Project Age Group") +
scale_colour_manual(values=legend)
g
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38),colour = "#E69F00", size=1)+
geom_abline(aes(intercept=1.65, slope=-1.31), colour = "#56B4E9", size=1)+
geom_abline(aes(intercept=1.65, slope=-0.78), colour = "#D55E00", size=1)+
geom_abline(aes(intercept=1.65, slope=-0.23) , colour = "#CC79A7", size=1)+
labs(x="MMT", y="Mean Underproduction Factor", color="Project Age Group") +
scale_colour_manual(values=legend)
g
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38),colour = "#E69F00", size=1)+
geom_abline(aes(intercept=1.65, slope=-1.31), colour = "#56B4E9", size=1)+
geom_abline(aes(intercept=1.65, slope=-0.78), colour = "#D55E00", size=1)+
geom_abline(aes(intercept=1.65, slope=-0.23) , colour = "#CC79A7", size=1)+
labs(x="MMT", y="Mean Underproduction Factor", colour="Project Age Group") +
scale_colour_manual(values=legend)
g
colors <- c("0-9y":"#E69F00","9-12y":"#56B4E9", "12-15y":"#D55E00","15-16y":"#CC79A7")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38),colour = "#E69F00", size=1)+
geom_abline(aes(intercept=1.65, slope=-1.31), colour = "#56B4E9", size=1)+
geom_abline(aes(intercept=1.65, slope=-0.78), colour = "#D55E00", size=1)+
geom_abline(aes(intercept=1.65, slope=-0.23) , colour = "#CC79A7", size=1)+
labs(x="MMT", y="Mean Underproduction Factor") +
scale_colour_manual(name="Project Age Group", values=c("0-9y"="#E69F00","9-12y"="#56B4E9", "12-15y"="#D55E00","15-16y"="#CC79A7"))
g
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38),colour = "#E69F00", size=1)+
geom_abline(aes(intercept=1.65, slope=-1.31), colour = "#56B4E9", size=1)+
geom_abline(aes(intercept=1.65, slope=-0.78), colour = "#D55E00", size=1)+
geom_abline(aes(intercept=1.65, slope=-0.23) , colour = "#CC79A7", size=1)+
labs(x="MMT", y="Mean Underproduction Factor") +
scale_colour_manual(name="Project Age Group", values=c("0-9y"="#E69F00","9-12y"="#56B4E9", "12-15y"="#D55E00","15-16y"="#CC79A7"))
g
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, colour=new.age.factor), size=1)+
geom_abline(aes(intercept=1.65, slope=-1.38),colour = "#E69F00", size=1)+
geom_abline(aes(intercept=1.65, slope=-1.31), colour = "#56B4E9", size=1)+
geom_abline(aes(intercept=1.65, slope=-0.78), colour = "#D55E00", size=1)+
geom_abline(aes(intercept=1.65, slope=-0.23) , colour = "#CC79A7", size=1)+
labs(x="MMT", y="Mean Underproduction Factor") +
scale_colour_manual(name="Project Age Group", values=c("0-9y"="#E69F00","9-12y"="#56B4E9", "12-15y"="#D55E00","15-16y"="#CC79A7"))
g
data2 <- subset(data1, (data1$age / 365) < 14 )
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.5, colour=new.age.factor), size=1)+
geom_abline(aes(intercept=1.65, slope=-1.38),colour = "#E69F00", size=1)+
geom_abline(aes(intercept=1.65, slope=-1.31), colour = "#56B4E9", size=1)+
geom_abline(aes(intercept=1.65, slope=-0.78), colour = "#D55E00", size=1)+
geom_abline(aes(intercept=1.65, slope=-0.23) , colour = "#CC79A7", size=1)+
labs(x="MMT", y="Mean Underproduction Factor") +
scale_colour_manual(name="Project Age Group", values=c("0-9y"="#E69F00","9-12y"="#56B4E9", "12-15y"="#D55E00","15-16y"="#CC79A7"))
g
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.5, colour=new.age.factor), size=1)+
geom_abline(aes(intercept=1.65, slope=-1.7, colour=new.age.factor), size=1)+
geom_abline(aes(intercept=1.65, slope=-1.38),colour = "#E69F00", size=1)+
geom_abline(aes(intercept=1.65, slope=-1.31), colour = "#56B4E9", size=1)+
geom_abline(aes(intercept=1.65, slope=-0.78), colour = "#D55E00", size=1)+
geom_abline(aes(intercept=1.65, slope=-0.23) , colour = "#CC79A7", size=1)+
labs(x="MMT", y="Mean Underproduction Factor") +
scale_colour_manual(name="Project Age Group", values=c("0-9y"="#E69F00","9-12y"="#56B4E9", "12-15y"="#D55E00","15-16y"="#CC79A7"))
g
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38),colour = "#E69F00", size=1)+
geom_abline(aes(intercept=1.65, slope=-1.31), colour = "#56B4E9", size=1)+
geom_abline(aes(intercept=1.65, slope=-0.78), colour = "#D55E00", size=1)+
geom_abline(aes(intercept=1.65, slope=-0.23) , colour = "#CC79A7", size=1)+
labs(x="MMT", y="Mean Underproduction Factor") +
scale_colour_manual(name="Project Age Group", values=c("0-9y"="#E69F00","9-12y"="#56B4E9", "12-15y"="#D55E00","15-16y"="#CC79A7"))
g
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38),colour = "#E69F00", size=2)+
geom_abline(aes(intercept=1.65, slope=-1.31), colour = "#56B4E9", size=1)+
geom_abline(aes(intercept=1.65, slope=-0.78), colour = "#D55E00", size=1)+
geom_abline(aes(intercept=1.65, slope=-0.23) , colour = "#CC79A7", size=1)+
labs(x="MMT", y="Mean Underproduction Factor") +
scale_colour_manual(name="Project Age Group", values=c("0-9y"="#E69F00","9-12y"="#56B4E9", "12-15y"="#D55E00","15-16y"="#CC79A7"))
g
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38),colour = "#E69F00", size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31), colour = "#56B4E9", size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78), colour = "#D55E00", size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23) , colour = "#CC79A7", size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor") +
scale_colour_manual(name="Project Age Group", values=c("0-9y"="#E69F00","9-12y"="#56B4E9", "12-15y"="#D55E00","15-16y"="#CC79A7"))
g
colors <- c("0-9y"="#E69F00","9-12y":"#56B4E9", "12-15y":"#D55E00","15-16y":"#CC79A7")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="0-9y"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="9-12y"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="12-15y" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="15-16y") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor") +
scale_colour_manual(values=colors_legend)
g
#colors <- c("0-9y"="#E69F00","9-12y":"#56B4E9", "12-15y":"#D55E00","15-16y":"#CC79A7")
colors_legend <- c("0-9y"="red","9-12y":"green", "12-15y":"blue","15-16y":"orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="0-9y"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="9-12y"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="12-15y" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="15-16y") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor") +
scale_colour_manual(values=colors_legend)
#colors <- c("0-9y"="#E69F00","9-12y":"#56B4E9", "12-15y":"#D55E00","15-16y":"#CC79A7")
colors_legend <- c("0-9y"="red","9-12y":"green", "12-15y":"blue","15-16y":"orange")
#colors <- c("0-9y"="#E69F00","9-12y":"#56B4E9", "12-15y":"#D55E00","15-16y":"#CC79A7")
colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="0-9y"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="9-12y"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="12-15y" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="15-16y") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor") +
scale_colour_manual(values=colors_legend)
g
colors_legend <- c("0-9y"="#E69F00","9-12y":"#56B4E9", "12-15y":"#D55E00","15-16y":"#CC79A7")
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="0-9y"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="9-12y"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="12-15y" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="15-16y") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor") +
scale_colour_manual(values=colors_legend)
g
colors_legend <- c("0-9y"="#E69F00","9-12y":"#56B4E9", "12-15y":"#D55E00","15-16y":"#CC79A7")
colors_legend <- c("0-9y"="#E69F00","9-12y"="#56B4E9", "12-15y"="#D55E00","15-16y"="#CC79A7")
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="0-9y"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="9-12y"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="12-15y" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="15-16y") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor") +
scale_colour_manual(values=colors_legend)
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="0-9y"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="9-12y"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="12-15y" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="15-16y") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Project Age Group") +
scale_colour_manual(values=colors_legend)
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="0-9y"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="9-12y"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="12-15y" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="15-16y") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Age Group") +
scale_colour_manual(values=colors_legend)
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="0-9y"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="9-12y"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="12-15y" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="15-16y") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Age Group") +
scale_colour_manual(values=colors_legend, breaks = c("1", "2", "3", "4"))
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="0-9y"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="9-12y"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="12-15y" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="15-16y") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Age Group") +
scale_colour_manual(values=colors_legend, breaks=c("1", "2", "3", "4"))
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color=0), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color=1), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color=2 ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color=3) , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Age Group") +
scale_colour_manual(values=colors_legend)
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="b"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="c" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Age Group") +
scale_colour_manual(values=colors_legend)
g
colors_legend <- c("a"="#E69F00","b"="#56B4E9", "c"="#D55E00","d"="#CC79A7")
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="b"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="c" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Age Group") +
scale_colour_manual(values=colors_legend)
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="b"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="c" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Age Group") +
scale_colour_manual(values=colors_legend, labels=c("0-9y", "9-12y", "12-15y","15-16y"))
g
+ theme(legend.position='bottom')
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="b"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="c" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Age Group") +
scale_colour_manual(values=colors_legend, labels=c("0-9y", "9-12y", "12-15y","15-16y")) +
theme(legend.position='bottom')
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="b"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="c" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Age Group") +
scale_colour_manual(values=colors_legend, labels=c("0-9y", "9-12y", "12-15y","15-16y")) +
theme(legend.position=c(0.95, 0.95))
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="b"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="c" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Age Group") +
scale_colour_manual(values=colors_legend, labels=c("0-9y", "9-12y", "12-15y","15-16y")) +
theme(legend.position=c(1, -5))
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="b"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="c" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Age Group") +
scale_colour_manual(values=colors_legend, labels=c("0-9y", "9-12y", "12-15y","15-16y")) +
theme(legend.position=c(-5, 1))
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="b"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="c" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Age Group") +
scale_colour_manual(values=colors_legend, labels=c("0-9y", "9-12y", "12-15y","15-16y")) +
theme(legend.position="bottom")
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="b"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="c" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Age Group") +
scale_colour_manual(values=colors_legend, labels=c("0-9y", "9-12y", "12-15y","15-16y")) +
theme( legend.justification = c("right", "top"))
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="b"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="c" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Age Group") +
scale_colour_manual(values=colors_legend, labels=c("0-9y", "9-12y", "12-15y","15-16y")) +
theme( legend.justification = c("left", "bottom"))
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="b"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="c" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Age Group") +
scale_colour_manual(values=colors_legend, labels=c("0-9y", "9-12y", "12-15y","15-16y")) +
theme( legend.position = c(0.95, 0.95), legend.justification = c("left", "bottom"))
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="b"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="c" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Age Group") +
scale_colour_manual(values=colors_legend, labels=c("0-9y", "9-12y", "12-15y","15-16y")) +
theme( legend.position = c(0.75, 0.95), legend.justification = c("left", "bottom"))
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="b"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="c" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Age Group") +
scale_colour_manual(values=colors_legend, labels=c("0-9y", "9-12y", "12-15y","15-16y")) +
theme( legend.position = c(0.35, 0.95), legend.justification = c("left", "bottom"))
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="b"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="c" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Age Group") +
scale_colour_manual(values=colors_legend, labels=c("0-9y", "9-12y", "12-15y","15-16y")) +
theme( legend.position = c(0.35, 0.25), legend.justification = c("left", "bottom"))
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="b"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="c" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Age Group") +
scale_colour_manual(values=colors_legend, labels=c("0-9y", "9-12y", "12-15y","15-16y")) +
theme( legend.position = c(0.35, 0.15), legend.justification = c("left", "bottom"))
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="b"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="c" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Age Group") +
scale_colour_manual(values=colors_legend, labels=c("0-9y", "9-12y", "12-15y","15-16y")) +
theme( legend.position = c(0.25, 0.15), legend.justification = c("left", "bottom"))
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="b"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="c" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Age Group") +
scale_colour_manual(values=colors_legend, labels=c("0-9y", "9-12y", "12-15y","15-16y")) +
theme( legend.position = c(0.15, 0.15), legend.justification = c("left", "bottom"))
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="b"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="c" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Age Group") +
scale_colour_manual(values=colors_legend, labels=c("0-9y", "9-12y", "12-15y","15-16y")) +
theme( legend.position = c(0.15, 0.05), legend.justification = c("left", "bottom"))
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="b"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="c" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Age Group") +
scale_colour_manual(values=colors_legend, labels=c("0-9y", "9-12y", "12-15y","15-16y")) +
theme( legend.position = c(0.05, 0.05), legend.justification = c("left", "bottom"))
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.31, color="b"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.78, color="c" ), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-0.23, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Project Age Group") +
scale_colour_manual(values=colors_legend, labels=c("0-9y", "9-12y", "12-15y","15-16y")) +
theme( legend.position = c(0.05, 0.05), legend.justification = c("left", "bottom"))
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.72, slope=-1.38, color="b"), size=1.5)+
geom_abline(aes(intercept=1.65, slope=-1.38, color="c" ), size=1.5)+
geom_abline(aes(intercept=2.8, slope=-1.38, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Project Age Group") +
scale_colour_manual(values=colors_legend, labels=c("0-9y", "9-12y", "12-15y","15-16y")) +
theme( legend.position = c(0.05, 0.05), legend.justification = c("left", "bottom"))
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.72, slope=-1.38, color="b"), size=1.5)+
geom_abline(aes(intercept=2.255, slope=-1.38, color="c" ), size=1.5)+
geom_abline(aes(intercept=2.8, slope=-1.38, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Project Age Group") +
scale_colour_manual(values=colors_legend, labels=c("0-9y", "9-12y", "12-15y","15-16y")) +
theme( legend.position = c(0.05, 0.05), legend.justification = c("left", "bottom"))
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.72, slope=-1.38, color="b"), size=1.5)+
geom_abline(aes(intercept=2.25, slope=-1.38, color="c" ), size=1.5)+
geom_abline(aes(intercept=2.8, slope=-1.38, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Project Age Group") +
scale_colour_manual(values=colors_legend, labels=c("0-9y", "9-12y", "12-15y","15-16y")) +
theme( legend.position = c(0.05, 0.05), legend.justification = c("left", "bottom"))
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.72, slope=-1.38, color="b"), size=1.5)+
geom_abline(aes(intercept=2.25, slope=-1.38, color="c" ), size=1.5)+
geom_abline(aes(intercept=2.8, slope=-1.38, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Project Age Group") +
scale_colour_manual(values=colors_legend, labels=c("0-9y", "9-12y", "12-15y","15-16y")) +
theme( legend.position = c(0.05, 0.05), legend.justification = c("left", "bottom")) +
theme_bw()
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.72, slope=-1.38, color="b"), size=1.5)+
geom_abline(aes(intercept=2.25, slope=-1.38, color="c" ), size=1.5)+
geom_abline(aes(intercept=2.8, slope=-1.38, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Project Age Group") +
scale_colour_manual(values=colors_legend, labels=c("0-9y", "9-12y", "12-15y","15-16y")) +
theme_bw( legend.position = c(0.05, 0.05), legend.justification = c("left", "bottom"))
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.72, slope=-1.38, color="b"), size=1.5)+
geom_abline(aes(intercept=2.25, slope=-1.38, color="c" ), size=1.5)+
geom_abline(aes(intercept=2.8, slope=-1.38, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Project Age Group") +
scale_colour_manual(values=colors_legend, labels=c("0-9y", "9-12y", "12-15y","15-16y")) +
theme_bw(legend.position = c(0.05, 0.05), legend.justification = c("left", "bottom"))
g
#colors_legend <- c("0-9y"="red","9-12y"="green", "12-15y"="blue","15-16y"="orange")
g <- ggplot(data1, aes(x=mmt, y=up.fac.mean)) +
geom_point() +
geom_abline(aes(intercept=1.65, slope=-1.38, color="a"), size=1.5)+
geom_abline(aes(intercept=1.72, slope=-1.38, color="b"), size=1.5)+
geom_abline(aes(intercept=2.25, slope=-1.38, color="c" ), size=1.5)+
geom_abline(aes(intercept=2.8, slope=-1.38, color="d") , size=1.5)+
labs(x="MMT", y="Mean Underproduction Factor", color = "Project Age Group") +
scale_colour_manual(values=colors_legend, labels=c("0-9y", "9-12y", "12-15y","15-16y")) +
theme_bw()+
theme(legend.position = c(0.05, 0.05), legend.justification = c("left", "bottom"))
g
