############################# ## Einheit2 (Wiederholung) ## ############################# ################################################### ELECTION <- read.csv2(file = "election02.csv") ELECTION attach(ELECTION) barplot(PERCENT, names = as.character(PARTY), col = 2) VOTES sum(VOTES[1:4]) binom.test(105, 325) binom.test(105, 325, conf.level = 0.99)$conf.int ELECTION2 <- read.csv2(file = "election02_frequ.csv") ELECTION2 detach(ELECTION) remove(ELECTION, ELECTION2) ################################################### CMMRCIAL <- read.csv2("cmmrcial.csv") dim(CMMRCIAL) summary(CMMRCIAL) attach(CMMRCIAL) CEREAL2 <- factor(ifelse(CEREAL == "CC", "CC", "other")) tab <- table(CEREAL2, WERBUNG) tab prop.table(tab, 2) barplot(prop.table(tab, 2), legend = TRUE, xlim = c(0, 3)) chisq.test(tab) detach(CMMRCIAL) remove(tab, CMMRCIAL, CEREAL2) ################################################### PROMOT <- read.csv2("promot.csv") PROMOT tab <- xtabs(COUNT ~ TRY + PROMOT, data = PROMOT) tab prop.table(tab, 2) barplot(prop.table(tab, 2)) fisher.test(tab) tab tab[2:1, ] fisher.test(tab[2:1,]) or1 <- fisher.test(tab)$estimate or2 <- fisher.test(tab[2:1,])$estimate 1/or2 remove(PROMOT, tab, or1, or2) ################################################### AUTO <- read.csv2(file = "autoextras.csv") tab <- xtabs(COUNTS ~ KLIMA + ALTER, data = AUTO) tab prop.table(tab, 2) barplot(prop.table(tab, 2)) mosaicplot(t(tab)) chisq.test(tab) chisq.test(tab)$expected chisq.test(tab)$observed chisq.test(tab)$residuals mosaicplot(t(tab), shade = TRUE) remove(AUTO, tab) ############## ## Einheit3 ## ############## ## Metrische Merkmale ############################# NOEWASSER <- read.table("noewasser.tab", header = TRUE) dim(NOEWASSER) names(NOEWASSER) attach(NOEWASSER) mean(NITRAT) median(NITRAT) var(NITRAT) sd(NITRAT) mad(NITRAT) range(NITRAT) diff(range(NITRAT)) IQR(NITRAT) fivenum(NITRAT) summary(NITRAT) table(NITRAT) barplot(table(NITRAT)) hist(NITRAT) hist(NITRAT, prob = TRUE) source("tcltk.R") tkHist(NITRAT) boxplot(NITRAT) detach(NOEWASSER) remove(NOEWASSER) ## Lagemasse ###################################### LAGE <- read.table("lagemasse.tab", header = TRUE) dim(LAGE) names(LAGE) attach(LAGE) hist(TEST, prob = TRUE) hist(TEST, prob = TRUE, breaks = 1:9) abline(v = mean(TEST), col = 2) hist(KRANK, prob = TRUE) hist(KRANK, prob = TRUE, breaks = 0:76) abline(v = mean(KRANK), col = 2) hist(MAGAZIN, prob = TRUE) hist(MAGAZIN, prob = TRUE, breaks = 0:12) abline(v = mean(MAGAZIN), col = 2) mean(KRANK) mean(KRANK[-20]) hist(GEMUESE, prob = TRUE) boxplot(GEMUESE) boxplot(GEMUESE, range = 0) detach(LAGE) remove(LAGE) ################################################### AUTOS <- read.table("autos.tab", header = TRUE) dim(AUTOS) names(AUTOS) attach(AUTOS) hist(PS, prob = TRUE) qqnorm(PS) qqline(PS) detach(AUTOS) remove(AUTOS) ################################################### CARIDEX <- read.table("caridex.tab", header = TRUE) dim(CARIDEX) names(CARIDEX) attach(CARIDEX) t.test(CAVITIES) t.test(CAVITIES, conf.level = 0.99)$conf.int detach("CARIDEX") remove(CARIDEX) ################################################### DIET <- read.table("diet.tab", header = TRUE) dim(DIET) names(DIET) attach(DIET) t.test(GEWICHT, mu = 5, alternative = "less") detach("DIET") remove(DIET) ## 2 Metrische Merkmale ########################### ALCTOBAC <- read.table("alctobac.tab", header = TRUE) dim(ALCTOBAC) names(ALCTOBAC) attach(ALCTOBAC) plot(TOBACCO, ALCOHOL) cor(TOBACCO, ALCOHOL) outlierIndex <- which((ALCOHOL < 5) & (TOBACCO > 4)) outlierIndex ALCTOBAC[outlierIndex, ] cor(TOBACCO[-11], ALCOHOL[-11]) cor.test(TOBACCO, ALCOHOL) cor.test(TOBACCO[-11], ALCOHOL[-11]) detach(ALCTOBAC) remove(ALCTOBAC) ################################################### BICYCLE <- read.table("bicycle.tab", header = TRUE) dim(BICYCLE) names(BICYCLE) attach(BICYCLE) hist(KMNACH - KMVOR, prob = TRUE) abline(v = 0, col = 4) boxplot(KMNACH - KMVOR) abline(h = 0, col = 4) t.test(KMNACH, KMVOR, paired = TRUE, alternative = "less") detach("BICYCLE") remove(BICYCLE) ################################################### ALCATTR <- read.table("alcattrneu.tab", header = TRUE) dim(ALCATTR) names(ALCATTR) attach(ALCATTR) boxplot(NACHHER-VORHER) abline(h = 0, col = 4) wilcox.test(NACHHER, VORHER, paired = TRUE, alternative = "greater") detach("ALCATTR") remove(ALCATTR)