source("../../TOOLS/sim_plotfuns.R")

load("weeds_simdat.RData")

## slight massage
L3 <- transform(L2,
                method=reorder(method,elapsed,mean),
                elapsed=elapsed+1e-3)

S <- simplots(L3[,1:7],parnames=c("logb1","logb2","logb3","logSigma"),
         truevals=data.frame(variable=names(simpars),value=simpars),
         has_conf=FALSE)

print(S$time+scale_y_log10()+coord_flip())
print(S$params+coord_flip())

## CG does badly
## logSigma is underestimated throughout: MLE is known to be biased