

In [4]:

    
require(ggplot2)
require(reshape2)









    



Loading required package: reshape2

SA Likelihood comparison

The following likelihood comparison shows the likelihood of the removal probability paramter $r$ calculated by the PF(SSA), PF(TL), and the analytical result computed by the SA package.

Parameter values used in this comparison:

$\beta=1.2$
$\mu=0.2$
$\psi=0.2$
$\rho=0.1$
$t_0=7$ (measured here as time before most recent sample.)



In [97]:

    
df <- read.table("BD_DensityMap_PsiNoRemove.log", header=T)
df <- melt(df, measure.vars = c("density","densitySA","densityTL"))

ggplot(df) + geom_line(aes(r, value, col=variable)) +
    geom_point(aes(r, value, col=variable))