HelpersMG (version 5.1)

modeled.hist: Return the theoretical value for the histogram bar

Description

Return the theoretical value for the histogram bar based on a model of distribution.

Usage

modeled.hist(breaks, FUN, ..., sum = 1)

Arguments

breaks

Vector with the breaks; it can be obtained directly from hist()

FUN

Function to be used to integrate the density, ex. pnorm

...

Parameters to be used by FUN

sum

Total numbers in the histogram; 1 for emperical frequencies

Value

A list with x (the center of the bar) and y components

Details

modeled.hist returns the theoretical value for the histogram bar based on a model of distribution.

Examples

# NOT RUN {
n <- rnorm(100, mean=10, sd=2)
breaks <- 0:20
hist(n, breaks=breaks)

s <- modeled.hist(breaks=breaks, FUN=pnorm, mean=10, sd=2, sum=100)

points(s$x, s$y, pch=19)
lines(s$x, s$y)

n <- rlnorm(100, meanlog=2, sdlog=0.4)
b <- hist(n, ylim=c(0, 70))

s <- modeled.hist(breaks=b$breaks, FUN=plnorm, meanlog=2, sdlog=0.4, sum=100)

points(s$x, s$y, pch=19)
lines(s$x, s$y)
# }