In [ ]:

    

library(ggplot2)


    

In [ ]:

    

stats = read.csv("summary.csv")


    

In [ ]:

    

coverage = stats[stats$scale==0.01, ]
summary(coverage)


    

In [ ]:

    

ggplot(coverage, aes(x=coverage, y=spearman, group=measure)) +
    geom_line(aes(linetype=measure)) +
    xlab('Genome Coverage') +
    ylab(expression(paste("Spearman's ", rho))) +
    theme_bw()


    

In [ ]:

    

scale = stats[stats$coverage==30, ]
scale$scale = 1/scale$scale


    

In [ ]:

    

ggplot(scale, aes(x=scale, y=spearman, group=measure)) +
    geom_line(aes(linetype=measure)) +
    xlab('Divergence (1 variant per x bases)') +
    ylab(expression(paste("Spearman's ", rho))) +
    scale_x_reverse() +
    scale_x_log10() +
    theme_bw()


    

In [ ]:

    

?scale_x_reverse


    

In [ ]: