|
|
@ -282,3 +282,37 @@ interval.cumrelfreq = interval.cumfreq / nrow(faithful)
|
|
|
|
interval.cumrelfreq
|
|
|
|
interval.cumrelfreq
|
|
|
|
cbind(interval.cumfreq, interval.cumrelfreq)
|
|
|
|
cbind(interval.cumfreq, interval.cumrelfreq)
|
|
|
|
q()
|
|
|
|
q()
|
|
|
|
|
|
|
|
c(duration.freq)
|
|
|
|
|
|
|
|
duration.freq)
|
|
|
|
|
|
|
|
duration.freq
|
|
|
|
|
|
|
|
c(0, duration.freq)
|
|
|
|
|
|
|
|
cumfreq0
|
|
|
|
|
|
|
|
cumfreq0 = c(0, cumsum(duration.freq))
|
|
|
|
|
|
|
|
cumfreq0
|
|
|
|
|
|
|
|
cumrelfreq0 = cumfreq0 / nrow(faithful)
|
|
|
|
|
|
|
|
cumrelfreq0
|
|
|
|
|
|
|
|
breaks
|
|
|
|
|
|
|
|
plot(breaks, cumrelfreq0, main="Old Faithful Eruptions", xlab="Duration minutes", ylab="Cumulative Eruptions Proportion")
|
|
|
|
|
|
|
|
lines(breaks, cumrelfreq0)
|
|
|
|
|
|
|
|
library(actuar)
|
|
|
|
|
|
|
|
help("library")
|
|
|
|
|
|
|
|
library(actuar)
|
|
|
|
|
|
|
|
help(ogive)
|
|
|
|
|
|
|
|
Fn = ogive(breaks, duration.freq)
|
|
|
|
|
|
|
|
plot(Fn, main="Old Faithful Eruptions", xlab="Duration minutes", ylab="Cumulative Relative Frequency")
|
|
|
|
|
|
|
|
library(PBSModelling)
|
|
|
|
|
|
|
|
library(PBSModeling)
|
|
|
|
|
|
|
|
library(PBSmodelling)
|
|
|
|
|
|
|
|
plotCsum(duration)
|
|
|
|
|
|
|
|
intervals
|
|
|
|
|
|
|
|
intervals.cut
|
|
|
|
|
|
|
|
interval_breaks
|
|
|
|
|
|
|
|
interval
|
|
|
|
|
|
|
|
intervals
|
|
|
|
|
|
|
|
interval.cut
|
|
|
|
|
|
|
|
interval.freq
|
|
|
|
|
|
|
|
intervalFn = ogive(breaks, interval.freq)
|
|
|
|
|
|
|
|
intervalFn = ogive(interval_breaks, interval.freq)
|
|
|
|
|
|
|
|
plot(intervalFn, main="Old Faithful Eruptions", xlab="Intervals in minutes", ylab="Cumulative Relative Frequency")
|
|
|
|
|
|
|
|
plotCsum(interval)
|
|
|
|
|
|
|
|
q()
|
|
|
|