EIdist: Plot the empirical and exponential distributions of values

Description Usage Arguments Details Value Author(s) Examples

View source: R/EIdist.R

Description

Plot the empirical distribution of a set of values.

Usage

1
2
 EIdist(event_times,nbreaks=10,main="",xlab="",ylab="",
 xaxticks=NA,xaxlabs=NA)

Arguments

event_times

A numeric vector of times of occurrence of events.

nbreaks

The number of breaks desired (but not always delivered).

main

Title for the plot.

xlab

X axis label for the plot.

ylab

Y axis label for the plot.

xaxticks

Optional X tick mark positions.

xaxlabs

Optional X tick labels.

Details

EIdist calculates the intervals between the times passed in event_times and plots a histogram of their distribution. It adds a smoothed line based on the mids and counts. It then calls fitdistr from MASS to get an estimated rate for a model exponential distribution. A second line is added to the plot for comparison with the smoothed empirical distribution.

Value

nil

Author(s)

Jim Lemon

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
 ## Not run: 
  # won't pass check due to the call to fitdistr
  data(florida_hurr20)
  oldloc<-Sys.setlocale("LC_TIME", "C")
  fh_dates<-
   as.Date(paste(florida_hurr20$day,florida_hurr20$month,
   florida_hurr20$year,sep="-"),"%d-%B-%Y")
  fh_days<-as.numeric(fh_dates)
  fh_ints<-diff(fh_days)
  fh_counts<-tabulate(florida_hurr20$year[-1]-1899,nbins=100)
  require(MASS)
  EIdist(fh_days,main="Empirical distribution of hurricanes",
   xlab="Intervals in days",ylab="Frequencies")
  Sys.setlocale("LC_TIME", oldloc)
 
## End(Not run)

Example output

Loading required package: MASS
[1] "C"

eventInterval documentation built on May 29, 2017, 12:15 p.m.