diff --git a/elementary-statistics/.RData b/elementary-statistics/.RData
index 4c85a28..59b3fc3 100644
Binary files a/elementary-statistics/.RData and b/elementary-statistics/.RData differ
diff --git a/elementary-statistics/.Rhistory b/elementary-statistics/.Rhistory
index f727142..b19c9da 100644
--- a/elementary-statistics/.Rhistory
+++ b/elementary-statistics/.Rhistory
@@ -170,3 +170,67 @@ interval.freq = table(interval.cut)
 interval.cumfreq = cumsum(interval.freq)
 hist(interval.cumfreq)
 q()
+duration = faithful$eruptions
+breaks = seq(1.5, 5.5, by=0.5)
+duration.cut = cut(duration, breaks, right=FALSE)
+duration.freq = table(duration.cut)
+cumfreq0 = c(0, cumsum(duration.freq))
+plot(breaks, cumfreq0, main="Old Faithful Eruptions", xlab="Duration minutes", ylab="Cumulative Eruptions")
+lines(breaks, cumfreq0)
+intervals = faithful$waiting
+faithful$waiting
+cbind(faithful$waiting)
+faithful$waiting / nrow(faithful)
+faithful$waiting / nrow(faithful$waiting)
+intervals
+breaks
+intervals.cut = cut(breaks, intervals)
+help(cut)
+intervals.cut = cut(intervals, breaks)
+intervals.cut
+intervals.cut
+intervals.freq = table(intervals.cut)
+intervals.cut = cut(intervals, breaks, right=FALSE)
+intervals.freq = table(intervals.cut)
+cumfreq0 = c(0, cumsum(intervals.freq))
+cumfreq0
+intervals
+intervals.cut = cut(intervals, breaks, right=FALSE)
+intervals.cut
+intervals.freq = table(intervals.cut)
+intervals.freq
+cbind(intervals.freq)
+plot(breaks, cumfreq0)
+intervals
+cut(intervals, breaks)
+cut(intervals, breaks, right=FALSE)
+duration
+cut(duration, breaks)
+breaks = seq(0, 100, by=10)
+breaks
+cut(intervals, breaks)
+intervals.cut = cut(intervals, breaks)
+intervals.cut = cut(intervals, breaks, right=FALSE)
+intervals.cut
+intervals.freq = table(intervals.cut)
+cumfreq0 = c(0, cumsum(intervals.freq))
+cumfreq0
+plot(breaks, cumfreq0)
+lines(breaks, cumfreq0)
+intervals
+min(intervals)
+max(intervals)
+breaks = seq(40.0, 100.0, by=5.0)
+intervals.cut = cut(intervals, breaks, right=FALSE)
+intervals.cut
+floor?
+help(floor)
+intervals.freq = table(intervals.cut)
+intervals.freq
+cumfreq0 = c(0, cumsum(intervals.freq))
+cumfreq0
+plot(breaks, cumfreq0)
+lines(breaks, cumfreq0)
+plot(breaks, cumfreq0, main="Old Faithful Intervals", xlab="Interval minutes", ylab="Cumulative Intervals")
+lines(breaks, cumfreq0)
+q()