# code for plotting phi vs length using # saved output for length, phi-hat, se-phi, etc from MARK to Excel # and stored as a csv file library(ggplot2) dat <- read.csv("TopModel_real_phi-length-grp.csv") head(dat) # rename 2nd column colnames(dat)[2] <- "phi" ggplot(dat, aes(x = length, y = phi, group = grp)) + geom_line(aes(color = grp)) + geom_ribbon(aes(ymin = LCI, ymax = UCI), alpha = 0.2) + scale_colour_brewer(palette = "Set2") + theme(legend.position = c(0.15, 0.8), legend.justification = c(1, 0)) + xlab("Body Length") + ylab(expression(hat(phi)))