diff options
| -rw-r--r-- | R/seqB.R | 6 | ||||
| -rw-r--r-- | man/seqB.Rd | 6 |
2 files changed, 8 insertions, 4 deletions
@@ -59,12 +59,14 @@ #' abline(h=1/(20/0.01+1), col="red") #' abline(v=res1$Rhat, col="blue") #' abline(v=Rpost$supp[loc], col="purple") -#' legend("topright", legend=c("prior", "posterior", "posterior mean (Rhat)", "posterior mode"), col=c("red", "black", "blue", "purple"), lty=1) +#' legend("topright", legend=c("prior", "posterior", "posterior mean (Rhat)", "posterior mode"), +#' col=c("red", "black", "blue", "purple"), lty=1) #' plot(Rpost$supp, Rpost$pmf, col="black", type="l", xlim=c(0.5, 1.5), xlab="", ylab="") #' abline(h=1/(20/0.01+1), col="red") #' abline(v=res1$Rhat, col="blue") #' abline(v=Rpost$supp[loc], col="purple") -#' legend("topright", legend=c("prior", "posterior", "posterior mean (Rhat)", "posterior mode"), col=c("red", "black", "blue", "purple"), lty=1) +#' legend("topright", legend=c("prior", "posterior", "posterior mean (Rhat)", "posterior mode"), +#' col=c("red", "black", "blue", "purple"), lty=1) #' #' ## ========================================================= ## #' ## Compute Rhat using only the first five weeks of data ## diff --git a/man/seqB.Rd b/man/seqB.Rd index 0628a89..184f476 100644 --- a/man/seqB.Rd +++ b/man/seqB.Rd @@ -71,12 +71,14 @@ plot(Rpost$supp, Rpost$pmf, col="black", type="l", xlab="", ylab="") abline(h=1/(20/0.01+1), col="red") abline(v=res1$Rhat, col="blue") abline(v=Rpost$supp[loc], col="purple") -legend("topright", legend=c("prior", "posterior", "posterior mean (Rhat)", "posterior mode"), col=c("red", "black", "blue", "purple"), lty=1) +legend("topright", legend=c("prior", "posterior", "posterior mean (Rhat)", "posterior mode"), + col=c("red", "black", "blue", "purple"), lty=1) plot(Rpost$supp, Rpost$pmf, col="black", type="l", xlim=c(0.5, 1.5), xlab="", ylab="") abline(h=1/(20/0.01+1), col="red") abline(v=res1$Rhat, col="blue") abline(v=Rpost$supp[loc], col="purple") -legend("topright", legend=c("prior", "posterior", "posterior mean (Rhat)", "posterior mode"), col=c("red", "black", "blue", "purple"), lty=1) +legend("topright", legend=c("prior", "posterior", "posterior mean (Rhat)", "posterior mode"), + col=c("red", "black", "blue", "purple"), lty=1) ## ========================================================= ## ## Compute Rhat using only the first five weeks of data ## |