Grafiken und Statistik in R

1.2 Blitzstart 1 Allgemeines
8
# Wie gibt man alle Spalten ohne die dritte an?
airquality[,-3] # alle Spalten ohne dritte: einfach mittels '-' (minus-Zeichen)
# Ozone Solar.R Temp Month Day
#1 41 190 67 5 1
#2 36 118 72 5 2
#3 12 149 74 5 3
#4 ... ... ... .. ..
# Wie gibt man Spalten 1, 3 und 6 an?
airquality[,c(1, 3, 6)] # Spalten 1, 3 und 6 durch Benutzen der Funktion c()
# 'c' steht für combine
# einfaches Zeichnen mit der Funktion plot() s.a. „Blasendiagramm“-Bsp. auf Seite 58
plot(airquality) # erzeugt eine pairs() - Grafik (s. auf Seite 57)
# Wie zeichnen wir nun einzelne Variablen? Erinnern wir uns an '$'
plot(airquality$Ozone~airquality$Temp) # allgemein: plot(y~x) WICHTIG: Tilde ~
# übersichtlicher wird es nach ausführen von attach()
attach(airquality) # Datensatz in den Suchpfad aufnehmen
par(no.readonly=TRUE) -> original # alte Grafikeinstellungen speichern
layout(matrix(c(1, 1, 2, 3), 2, 2) -> m); m # spezielles layout:
# [,1] [,2]
#[1,] 1 2
#[2,] 1 3
# Grafikfläche hat jetzt 3 Grafiken: li eine, re zwei übereinander s. auf Seite 31
plot(Ozone~Temp, # airquality$Ozone ist jetzt nicht mehr nötig, R weiß wo zu suchen ist
col=rainbow(5)[Month-4], # 5 Regenbogenfarben zeichnen bezüglich Monat
pch=16, # Punkttyp gefüllt s. auf Seite 30
# main="Luftdaten New York \n 1973", # Titelei - \n ist gleich Zeilenumbruch
xlab="Temperatur (°F)", # x-Achsenbeschriftung
ylab="Ozon ppb" # y-Achsenbeschriftung
) # Funktion plot() hier zu Ende
rug(Ozone, side=2) # gibt Anzahl der Beobachtungen als Randplot aus
� Zur Farbe: rainbow(5) liefert an sich nur eine Liste mit 5 Farben à la "#FF0000FF" "#CCFF00FF" usw, die intern fortlaufend
indiziert werden. rainbow(5)[1] ist also die erste Farbe: "#FF0000FF" , rainbow(5)[2] zweite Farbe: "#CCFF00FF" usw.
Damit die richtige Farbe zum richtigen Monat kommt (Monate 5, 6, 7, 8, 9), wird einfach in der eckigen Index-Klammer direkt nach
rainbow(5) 'Month-4' gerechnet. So haben wir jetzt nicht mehr eine Liste 5 5 5 5 .. 6 6 6 .. 7 7 7 ... 8 .. 9 9, sondern 1 1 1
1 .. 2 2 2 .. 3 3 3 ... 4 .. 5 5 – 1 bis 5 – und in der eckigen Index-Klammer [Month-4] von rainbow(5) wird also für Monat 5
die Zahl 1 zurückgegeben – und damit die erste Farbe, für Monat 6 die Zahl 2 – und damit die zweite Farbe – usw.
# Daten übersichtlicher mit Funktion boxplot(...) s. auf Seite 51
par(mar=c(0,4.1,4.1,2.1))
# Rand reduzieren bottom, left, top, right; s. S. 28
boxplot(Ozone~Month, # Month ist hier die Variable nach der gruppiert wird
col=rainbow(5), # 5 Regenbogenfarben zeichnen bezüglich Monat
xlab="", # x-Achsenbeschriftung
ylab="Ozon ppb", # y-Achsenbeschriftung
# main="Luftdaten New York \n 1973", # Titelei - \n ist gleich Zeilenumbruch
notch=TRUE, # Vertrauensintervall des Median
xaxt = "n" # x-Achse ohne 'ticks' und 'labels'
) # Funktion boxplot() hier zu Ende
...Fortsetzung umseitig