durations.hpd {BaM}R Documentation

durations.hpd

Description

Simple HPD calculator from Chapter 2 (page 51).

Usage

simple.hpd(support,fn.eval,start,stop,target=0.90,tol=0.01)

Arguments

support x-axis values
fn.eval function values at x-axis points
start starting point in the vectors
stop stoppng point in the vectors
target Desired X Level
tol Tolerance for round-off

Examples

## Not run: 
{
data(cabinet.duration)
  ruler <- seq(0.45,0.75,length=10000)
  g.vals <- round(dgamma(ruler,shape=sum(N), rate=sum(N*dur)),2)
  start.point  <- 1000; stop.point <- length(g.vals)
  simple.hpd(ruler,g.vals,start.point,stop.point)
}
## End(Not run)

[Package BaM version 0.96 Index]