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()