Tutorial Bootstrap

Question 1

Was sind die 3 Schritte um ein Bootstrap in R zu berechnen?

Answer 1

1. Erstellen einer Funktion zur Berechnung des interessierenden Kennwerts
2. Erstellen eines Bootstrap-Objektes, das u. a. die Berechnungen des interessierenden Kennwertes in allen Bootstrap-Stichproben enthält
3. Erstellen eines Konfidenzintervalls auf Grundlage des Bootstrap-Objektes

Question 2

Mit welchen Befehl erstellt man ein Konfidenzintervall für den Mittelwert mittels Bootstrapping?

Answer 2

1. m.alter.fun <-function(piraten, i){
   mean(piraten$alter[i], na.rm = TRUE)}

2 m.alter.object <- boot(piraten, m.alter.fun, 2000)

3. boot.ci(m.alter.object, type = “bca”, conf = .90)

Question 3

Mit welchen Befehlen erstellt man ein Konfidenzintervall für die Standardabweichung mittels Bootstrapping?

Answer 3

1. sd.schiffe.fun <- function(piraten,i){sd(piraten$schiffe[i], na.rm = TRUE)}
2. sd.schiffe.object <- boot(piraten, sd.schiffe.fun, 2000)
3. boot.ci(sd.schiffe.object, type = “bca”, conf=.90)

Question 4

Mit welchen Befehlen erstellt man ein Konfidenzintervall für den Median mittels Bootstrapping?

Answer 4

1. me.alter.fun <- function (piraten,i){median(piraten$alter[i], na.rm=TRUE)}
2. me.alter.object <-boot(piraten, me.alter.fun,2000)

3.boot.ci.(me.alter.object. type = “bca”, conf=.90)

Question 5

Mit welchen Befehlen erstellt man ein Konfidenzintrvall für den Korrelationskoeffzienten mittels Bootstrapping?

Answer 5

1. r.schaetze.schiffe.fun <- function (piraten,i){cor(piraten$schaetze[i]piraten$schiffe[i], use = “pairwise.complete.obs”)}
2. r.schaetze.schiffe.object <-boot(piraten, r.schaetze.schiffe.fun,2000)

3.boot.ci(r.schaetze.schiffe.object, use “bca”, conf=.90)

Question 6

Mit welchen Befehlen erstellt man ein Konfidenzintrvall für den Determinationskoeffizienten mittels Bootstrapping?

Answer 6

1. r2.schiffe.alter.fun <-function(piraten,i){summary(lm(piraten$schiffe[i] ~piraten$alter[i]))$r.squared}
2. r2.schiffe.alter.object <-boot(piraten, r.schiffe.alter.fun, 2000)
3. boot.ci(r2.schiffe.alter.object, use = “bca”, conf=.90)

Question 7

Mit welchen Befehlen erstellt man ein Konfidenzintrvall für den Regressionskoeffizienten mittels Bootstrapping?

Answer 7

1. b.schiffe.alter.fun <- function(piraten,i){lm(piraten$schiffe[i] ~piraten$alter[i]$coefficients[2,1]}
2. b.schiffe.alter.object <-boot(piraten, b.schiffe.alter.fun,2000)
3. boot.ci(b.schiffe.alter.object, use = “bca”, conf=.90)