a grand renaming so that the most significant portion of the name comes first
This commit is contained in:
Dan Buch
367
r-practice/elementary-statistics/.Rhistory
Normal file
367
r-practice/elementary-statistics/.Rhistory
Normal file
@@ -0,0 +1,367 @@
|
||||
library(MASS)
|
||||
painters
|
||||
painters$School
|
||||
school = painters$School
|
||||
school.freq = table(school)
|
||||
school.freq
|
||||
cbind(school.freq)
|
||||
head(painters)
|
||||
composition = painters$Composition
|
||||
composition.freq = table(composition)
|
||||
cbind(composition.freq)
|
||||
school
|
||||
school.freq
|
||||
school.freq[1]
|
||||
help(max)
|
||||
school[,][1]
|
||||
school[1,][1]
|
||||
school[1][1]
|
||||
school[1]
|
||||
school.freq[1,]
|
||||
school.freq[1]
|
||||
school.freq[1][1]
|
||||
school.freq[1]
|
||||
help(typeof)
|
||||
typeof(school.freq[1])
|
||||
max(school.freq)
|
||||
max(school)
|
||||
max(school.freq)
|
||||
school.relfreq = school.freq / nrow(painters)
|
||||
school.relfreq
|
||||
settings = options(digits=1)
|
||||
school.relfreq
|
||||
options(settings)
|
||||
cbind(school.relfreq)
|
||||
composition.freq
|
||||
composition.relfreq = composition.freq / nrow(composition)
|
||||
composition.relfreq = composition.freq / nrow(painters)
|
||||
composition.relfreq
|
||||
cbind(composition.relfreq)
|
||||
# bar graph exercise
|
||||
library(MASS)
|
||||
school = painters$School
|
||||
school.freq = table(school)
|
||||
barplot(school.freq)
|
||||
colors = c("red", "yellow", "green", "violet", "orange", "blue", "pink", "cyan")
|
||||
barplot(schools.freq, col=colors)
|
||||
barplot(school.freq, col=colors)
|
||||
composition = painters$Composition
|
||||
composition.freq = table(composition)
|
||||
count.fields(composition.freq)
|
||||
ncol(composition.freq)
|
||||
composition.freq
|
||||
nrow(composition.freq)
|
||||
nrow(school.freq)
|
||||
colors
|
||||
barplot(composition.freq, col=colors)
|
||||
q()
|
||||
library(MASS)
|
||||
school = painters$school
|
||||
school
|
||||
school = painters$School
|
||||
school
|
||||
school.freq = table(school)
|
||||
pip(school.freq)
|
||||
pie(school.freq)
|
||||
colors = c("red", "yellow", "green", "violet", "orange", "blue", "pink", "cyan")
|
||||
pie(school.freq, col=colors)
|
||||
composition = painters$Composition
|
||||
composition.freq = table(composition)
|
||||
pip(composition.freq, col=colors)
|
||||
pie(composition.freq, col=colors)
|
||||
library(MASS)
|
||||
school = painters$School
|
||||
c_school = school == "C"
|
||||
c_painters = painters[c_school,]
|
||||
mean(c_painters$Composition)
|
||||
pie(tapply(painters$Composition, painters$School, mean)
|
||||
pie(tapply(painters$Composition, painters$School, mean))
|
||||
barchart(tapply(painters$Composition, painters$School, mean))
|
||||
barplot(tapply(painters$Composition, painters$School, mean))
|
||||
q()
|
||||
head(faithful)
|
||||
duration = faithful$eruptions
|
||||
range(duration)
|
||||
breaks = seq(1.5, 5.5, by=0.5)
|
||||
breaks
|
||||
duration.cut = cut(duration, breaks, right=FALSE)
|
||||
duration.freq = table(duration.cut)
|
||||
duration.freq
|
||||
cbind(duration.freq)
|
||||
barplot(duration.freq)
|
||||
nrow(faithful)
|
||||
faithful$eruptions / nrow(faithful)
|
||||
help(hist)
|
||||
hist(faithful$eruptions)
|
||||
faithful)
|
||||
head(faithful)
|
||||
hist(faithful$waiting)
|
||||
hist(faithful)
|
||||
help(barplot)
|
||||
VADeaths
|
||||
help(demo)
|
||||
demo()
|
||||
demo("graphics")
|
||||
hist(faithful$waiting)
|
||||
interval = faithful$waiting
|
||||
interval
|
||||
range(interval)
|
||||
breaks = seq(43, 96, by=1)
|
||||
breaks
|
||||
interval.cut = cut(interval, breaks, right=FALSE)
|
||||
interval.freq = table(interval.cut)
|
||||
barplot(interval.freq)
|
||||
breaks = seq(43, 96, by=5)
|
||||
breaks
|
||||
interval.cut = cut(interval, breaks, right=FALSE)
|
||||
interval.freq = table(interval.cut)
|
||||
barplot(interval.freq)
|
||||
help(points)
|
||||
points(interval.freq)
|
||||
barplot(interval.freq)
|
||||
hist(faithful$eruptions, right=FALSE)
|
||||
duration = faithful$eruptions
|
||||
hist(duration, right=FALSE)
|
||||
colors = c("red", "yellow", "green", "violet", "orange", "blue", "pink", "cyan")
|
||||
hist(duration, right=FALSE, col=colors, main="Old Faithful Eruptions", xlab="Duration minutes")
|
||||
waiting_period = faithful$waiting
|
||||
waiting_period = faithful$waiting
|
||||
hist(waiting_period, right=FALSE, col=colors, main="Old Faithful Eruptions", xlab="Interval")
|
||||
hist(waiting_period, right=FALSE, col=colors, main="Old Faithful Eruptions", xlab="Interval in minutes")
|
||||
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)
|
||||
duration.relfreq = duration.freq / nrow(faithful)
|
||||
duration.relfreq
|
||||
pie(duration.relfreq)
|
||||
hist(duration.relfreq)
|
||||
hist(duration.relfreq)
|
||||
duration.freq
|
||||
duration = faithful$eruptions
|
||||
breaks = seq(1.5, 5.5, by=0.5)
|
||||
breaks
|
||||
typeof(breaks)
|
||||
duration.cut = cut(duration, breaks, right=FALSE)
|
||||
duration.freq = table(duration.cut)
|
||||
duration.cumfreq = cumsum(duration.freq)
|
||||
duration.cumfreq
|
||||
hist(cbind(duration.cumfreq))
|
||||
interval = faithful$waiting
|
||||
help(function)
|
||||
help("function")
|
||||
help("function")
|
||||
calc_cumfreq <- function(records)
|
||||
breaks = seq(1.5, 5.5, by=0.5)
|
||||
records.cut = cut(duration, breaks, right=FALSE)
|
||||
records.freq = table(duration.cut)
|
||||
records.cumfreq = cumsum(records.freq)
|
||||
return(records)
|
||||
calc_cumfreq <- function(records)
|
||||
breaks = seq(1.5, 5.5, by=0.5)
|
||||
records.cut = cut(records, breaks, right=FALSE)
|
||||
help("function")
|
||||
duration.cumfreq
|
||||
interval = faithful$eruptions
|
||||
breaks
|
||||
breaks.cut = cut(interval, breaks, right=FALSE)
|
||||
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()
|
||||
duration = faithful$eruptions
|
||||
breaks = seq(1.5, 5.5, by=0.5)
|
||||
duration.cut = cut(duration, breaks, right=FALSE)
|
||||
duration.freq = table(duration.cut)
|
||||
duration.cumfreq = cumsum(duration.freq)
|
||||
duration.cumrelfreq = duration.cumfreq / nrow(faithful)
|
||||
duration.cumrelfreq
|
||||
old = options(digits=2)
|
||||
duration.cumrelfreq
|
||||
options(old)
|
||||
duration.cumrelfreq
|
||||
cbind(duration.cumfreq, duration.cumrelfreq)
|
||||
interval = faithful$waiting
|
||||
min(interval)
|
||||
max(interval)
|
||||
ceiling(max(interval))
|
||||
help(ceiling)
|
||||
interval * 0.10
|
||||
max(interval) * 0.1)
|
||||
ceiling(max(interval) * 0.1)
|
||||
interval_breaks = seq(floor(min(interval) * 0.1) * 10.0, ceiling(max(interval) * 0.1) * 10.0, right=FALSE)
|
||||
interval_breaks
|
||||
interval_breaks = seq(floor(min(interval) * 0.1) * 10.0, ceiling(max(interval) * 0.1) * 10.0, by=nrow(interval) / 10.0, right=FALSE)
|
||||
interval_breaks = seq(floor(min(interval) * 0.1) * 10.0, ceiling(max(interval) * 0.1) * 10.0, by=(nrow(faithful) / 10.0), right=FALSE)
|
||||
interval_breaks
|
||||
interval_breaks = seq(floor(min(interval) * 0.1) * 10.0, ceiling(max(interval) * 0.1) * 10.0, by=(nrow(faithful) / 5.0), right=FALSE)
|
||||
interval_breaks
|
||||
interval_breaks = seq(floor(min(interval) * 0.1) * 10.0, ceiling(max(interval) * 0.1) * 10.0, by=(nrow(faithful) / 20.0), right=FALSE)
|
||||
interval_breaks
|
||||
interval_breaks = seq(floor(min(interval) * 0.1) * 10.0, ceiling(max(interval) * 0.1) * 10.0, by=(nrow(faithful) / 30.0), right=FALSE)
|
||||
interval_breaks
|
||||
interval_breaks = seq(floor(min(interval) * 0.1) * 10.0, ceiling(max(interval) * 0.1) * 10.0, by=5.0 right=FALSE)
|
||||
interval_breaks = seq(floor(min(interval) * 0.1) * 10.0, ceiling(max(interval) * 0.1) * 10.0, by=5.0, right=FALSE)
|
||||
interval_breaks
|
||||
length(interval_breaks)
|
||||
type(interval_breaks)
|
||||
typeof(interval_breaks)
|
||||
interval
|
||||
interval.freq
|
||||
interval.cut = cut(interval, interval_breaks, right=FALSE)
|
||||
interval.cut
|
||||
interval.freq = table(interval.cut)
|
||||
interval.cumfreq = cumsum(interval.freq)
|
||||
interval.cumfreq
|
||||
interval.cumrelfreq = interval.cumfreq / nrow(faithful)
|
||||
interval.cumrelfreq
|
||||
cbind(interval.cumfreq, interval.cumrelfreq)
|
||||
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()
|
||||
duration
|
||||
stem(duration)
|
||||
interval
|
||||
stem(interval)
|
||||
interval
|
||||
help(stem)
|
||||
q()
|
||||
waiting = faithful$waiting
|
||||
head(cbind(duration, waiting))
|
||||
plot(duration, waiting, xlab="Eruption duration", ylab="Time waited")
|
||||
par()
|
||||
plot(1:10, 1:10, type="n")
|
||||
help("font")
|
||||
??font
|
||||
# the above was a brief excursion into font settings and such
|
||||
abline(lm(waiting - duration))
|
||||
plot(duration, waiting, xlab="Eruption duration", ylab="Time waited")
|
||||
abline(lm(waiting - duration))
|
||||
abline(lm(waiting ~ duration))
|
||||
q()
|
||||
mean(duration)
|
||||
mean(interval)
|
||||
??points
|
||||
median(duration)
|
||||
median(interval)
|
||||
quantile(duration)
|
||||
quantile(interval)
|
||||
quantile(duration, c(.32, .57, .98))
|
||||
quantile(interval, c(.17, .43, .67, .85))
|
||||
q()
|
||||
q()
|
||||
max(duration)
|
||||
min(duration)
|
||||
max(interval) - min(interval)
|
||||
IQR(duration)
|
||||
IQR(interval)
|
||||
boxplot(duration, horizontal=TRUE)
|
||||
boxplot(interval, horizontal=TRUE)
|
||||
q
|
||||
q()
|
||||
var(duration)
|
||||
var(interval)
|
||||
max(interval) - min(interval)
|
||||
sd(duration)
|
||||
sd(interval)
|
||||
waiting
|
||||
interval == waiting
|
||||
cov(duration, waiting)
|
||||
q()
|
||||
Reference in New Issue
Block a user