Improve display of confidence intervals

1 files changed, 8 insertions, 4 deletions

diff --git a/R/mkinparplot.R b/R/mkinparplot.R
index 8e8de02..7a3a6d4 100644
--- a/R/mkinparplot.R
+++ b/R/mkinparplot.R
@@ -38,16 +38,20 @@ mkinparplot <- function(object) {
     xlim = switch(type,
                   state.optim = range(c(0, unlist(values)), 
                                       na.rm = TRUE, finite = TRUE),
-                  rates.optim = range(unlist(values), 
+                  rates.optim = range(c(0, unlist(values)), 
                                       na.rm = TRUE, finite = TRUE),
                   fractions.optim = range(c(0, 1, unlist(values)), 
                                           na.rm = TRUE, finite = TRUE))
-    stripchart(bpar[get(type)], xlim = xlim,
-               yaxt = "n", ylim = c(0.5, length(get(type)) + 0.5))
+    stripchart(values["Estimate", ], 
+               xlim = xlim,
+               ylim = c(0.5, length(get(type)) + 0.5),
+               yaxt = "n")
+    if (type %in% c("rates.optim", "fractions.optim")) abline(v = 0, lty = 2)
+    if (type %in% c("fractions.optim")) abline(v = 1, lty = 2)
     text(mean(xlim), 1:length(parnames), parnames, pos = 3, offset = 1)
     arrows(as.numeric(values["Lower", ]), 1:length(parnames), 
            as.numeric(values["Upper", ]), 1:length(parnames), 
-           code = 3, angle = 90, length = 0.1)
+           code = 3, angle = 90, length = 0.05)
   }
   par(oldpars)
 }