Extrahiere die Bootstrap Replikla aus einem mzR Objekt

Falls ein mzR Objekt mit dem Parameter replicate = TRUE erzeugt wurde, lassen sich mit dieser Funktion alle Schätzwerte zurückgeben.

getReplicates(x)

Arguments

x	Ein Objekt der Klasse mzR. Typischerweise generiert durch GroupRate, Groupsize, Median, Mean oder Total.

Value

Tabelle mit Bootstrapreplika. Die Spalten entsprechen den Gruppierungsvariablen der Auswertung, falls vorhanden (siehe das Argument each in GroupRate, Median, etc.). Die Zeilen enstsprechen den Bootstrapgewichten (typischerweise 500 Stück).

Examples

library(dplyr)  ## Für %>%
#> 
#> Attaching package: ‘dplyr’
#> The following objects are masked from ‘package:data.table’:
#> 
#>     between, first, last
#> The following objects are masked from ‘package:stats’:
#> 
#>     filter, lag
#> The following objects are masked from ‘package:base’:
#> 
#>     intersect, setdiff, setequal, union
library(ggplot2)

######################## Beispiel 1: Durschnittliche Arbeitsstunden #############################

dat <- ImportData(year = 2014, quarter = 4)
#> Registered S3 method overwritten by 'openssl':
#>   method      from
#>   print.bytes Rcpp
#> Der Share ist bereits eingebunden!
#> re-encoding from CP1252
#> “/home/decill/mnt/mz_intern/2014/2014q4/dg7.mz2014q4.sav” wurde eingelesen.
#> Taking input= as a system command ('zcat /home/decill/mnt/mz_intern/2014/2014q4/mz2_2014q4_bootweights.csv.gz') and a variable has been used in the expression passed to `input=`. Please use fread(cmd=...). There is a security concern if you are creating an app, and the app could have a malicious user, and the app is not running in a secure envionment; e.g. the app is running as root. Please read item 5 in the NEWS file for v1.11.6 for more information and for the option to suppress this message.
#> '/home/decill/mnt/mz_intern/2014/2014q4/mz2_2014q4_bootweights.csv.gz'  wurde eingelesen.
mzObj <- Mean(dat, TFstring = "xerwstat==1&balt >= 15&balt <= 74",
              var = "estund*13+dtstd*13", replicates = TRUE)
replicates <- getReplicates(mzObj)$replicates

hist(replicates, main = "Durschnittlich geleistete Arbeitsstunden, 95% KI", freq = FALSE,
     col = "lightblue", xlab = "Durchschnittliche Wochenstunden in Replika")
abline(v = quantile(replicates, c(.025, 0.975)), col = "red", lwd = 2)
lines(density(replicates), col = "darkblue", lwd = 2)

########################### Beispiel 2: Wohnkosten nach Geschlecht ##############################

mzObj2 <- Mean(dat, TFstring = "xerwstat==1&balt >= 15&balt <= 74", var = "wkges",
               replicates = TRUE, each = "bsex")
getReplicates(mzObj2) %>% tidyr::gather(Geschlecht, wk) %>%
  mutate(Geschlecht = recode(Geschlecht, bsex_1 = "M", bsex_2 = "W")) %>%
  ggplot(aes(Geschlecht, wk, fill = Geschlecht)) +
  geom_boxplot(outlier.shape = NA, alpha = 0.5) +
  geom_jitter(aes(col = Geschlecht), alpha = 0.5) +
  ylab("Durchschnittliche Wohnkosten in Replika")

###################### Beispiel 3: Arbeitslosenquote nach Bundesland ############################

mzObj3 <- GroupSize(dat,TFstring="xerwstat==2&balt>=15&balt<=74", replicates = TRUE,
                    each = "xnuts2")
getReplicates(mzObj3) %>% tidyr::gather(Bundesland, unemployment) %>%
ggplot(aes(Bundesland, unemployment, fill = Bundesland)) +
  geom_boxplot(outlier.shape = NA, alpha = 0.5) +
  geom_jitter(aes(col = Bundesland), alpha = 0.05) +
  ylab("Arbeislosenzahlen in Bootstrap-Replika")