[1] Skript zum R-Crashkurs

###Datenimport###

# Eigener Datensatz kann über "Import Dataset" [siehe rechts "Environment"] geladen werden

# Wichtige Anmerkung: Kommas werden in Punkte umgewandelt, damit statistische Analysen funktionieren. Hierzu wurde auch der Generator für den Import gezeigt.

###Beispieldatensatz laden###

data("ToothGrowth")

###Datensatz kennenlernen/analysieren###

ls(ToothGrowth)

str(ToothGrowth)

View(ToothGrowth)

by(ToothGrowth$len, INDICES = list(ToothGrowth$supp, ToothGrowth$dose), summary)

plot (ToothGrowth$dose,ToothGrowth$len)

plot(ToothGrowth)

###Bestimmung von Datentypen###

# Anmerkung: Es gibt folgende Datentypen: numeric, complex, logical, character, raw

mode(ToothGrowth$len)

###Neue Variable generieren/löschen###

attach(ToothGrowth)

ToothGrowth$neuevariable <- ToothGrowth$len/ToothGrowth$dose

View(ToothGrowth)

ToothGrowth$neuevariable <- NULL

###Basisanalysen###

###Welch Two Sample Test Nr. 1: Ist Orangensaft effektiver als Ascorbinsäure? [Zahnwachstum]###

t.test(len~supp, data=ToothGrowth)

# Interpretation: Erster t-Test zeigt, dass generell OS effektiver als AS ist.Null-Hypothese (Differenz = 0) kann verworfen werden -> Indikation für H1. Signifikanz/p-Wert jedoch schwach (-> t-Test für dose = 2.0 zeigt Ursache hierfür)

t.test(len~supp, data=ToothGrowth[ToothGrowth$dose==0.5,])

t.test(len~supp, data=ToothGrowth[ToothGrowth$dose==1.0,])

t.test(len~supp, data=ToothGrowth[ToothGrowth$dose==2.0,])

###Welch Two Sample Test Nr. 2: Besteht eine Korrelation zwischen Zahnlänge und Dosierung?###

# Schritt 1: Dosierungen definieren, die verglichen werden sollen, z.B. dose1 = 0,5 und 1 mg

dose1 <- subset(ToothGrowth, dose %in% c(0.5, 1.0))

dose2 <- subset(ToothGrowth, dose %in% c(0.5, 2.0))

dose3 <- subset(ToothGrowth, dose %in% c(1.0, 2.0))

# Schritt 2: Tests für dose 1 bis 3

t.test(len ~ dose, paired = F, var.equal = F, data = dose1)

t.test(len ~ dose, paired = F, var.equal = F, data = dose2)

t.test(len ~ dose, paired = F, var.equal = F, data = dose3)

# Interpretation: Bei allen t-Tests kann die H0 verworfen werden -> Indikation für H1

###Korrelationsanalyse###

cor (ToothGrowth$len, ToothGrowth$dose, method = c("pearson"))

cor (ToothGrowth$len, ToothGrowth$dose, method = c("kendall"))

cor (ToothGrowth$len, ToothGrowth$dose, method = c("spearman"))

###Regressionsanalyse###

# Anmerkung: Datensatz ToothGrowth ist nicht geeignet für eine Regressionsanalyse -> Neuer Datensatz "swiss"

###Datensatz kennenlernen/analysieren###

# Folgende Variablen  sind enthalten:

# Fertility  Ig, "common standardized fertility measure"

# Agriculture % of males involved in agriculture as occupation

# Examination % draftees receiving highest mark on army examination

# Education % education beyond primary school for draftees.

# Catholic % "catholic" [as opposed to "protestant"].

# Infant.Mortality % live births who live less than 1 year.

data("swiss")

View(swiss)

plot(swiss)

cor(swiss)

###Einfache lineare Regression###

###Modellwahl###

# Anmerkung: Die Auswahl der Variablen, die ins Modell eingehen sollen, kann vorab per Akaikes Informationskriterium [AIC] erfolgen. Es gibt eine ganze Reihe weiterer Kriterien, die angelegt werden können.

step (lm(swiss$Fertility ~ swiss$Agriculture + swiss$Examination + swiss$Education + swiss$Catholic + swiss$Infant.Mortality))

regfert <- lm(formula = swiss$Fertility ~ swiss$Agriculture + swiss$Education + swiss$Catholic + swiss$Infant.Mortality)

###Ergebnisse###

regfert$coefficients

summary(regfert)

###Modellannahmen testen###

plot(regfert)

[2] Datensatzbeschreibungen "Toothgrowth" und "Swiss" [in Englisch]

Toothgrowth

• In 1947 data measuring the effect of vitamin C on tooth growth in guinea pigs  was collected.

• The response is the length of odontoblasts [cells responsible for tooth growth] in 60 guinea pigs. 

• Each animal received one of three dose levels of vitamin C [0.5, 1, and 2 mg/day] by one of two delivery methods, [orange juice or ascorbic acid (a form of vitamin C and coded as VC].

-> Data

Swiss

• Switzerland, in 1888, was entering a period known as the demographic transition; i.e., its fertility was beginning to fall from the high level typical of underdeveloped countries. 

• The data collected are for 47 French-speaking “provinces” at about 1888.

• Here, all variables are scaled to [0, 100], where in the original, all but "Catholic" were scaled to [0, 1].

-> Data

[3] Hinweise zu Statistik-Software

Open Source:

• R [freie Programmiersprache und zugleich Software-Umgebung]

• RStudio [integrierte Entwicklungsumgebung und grafische Benutzeroberfläche für R; Bedienung mittels syntaxbasierten "Scripts"]

• Weitere freie Softwares [z. T. auch R-basiert]

Proprietäre Software:

• Stata [Umfassende Software für empirische Analysen; Bedienung über Menü oder mittels syntaxbasierten "Do-Files"]

• SPSS [Umfassende Software für empirische Analysen; i. d. R. Bedienung über Menü]

• EViews [Einsatzschwerpunkt: Zeitreihenanalyse]

• MATLAB [Einsatzschwerpunkt: Numerische Berechnungen, insb. Matrizen]