1 require(Hmisc)
   2 
   3 x <- spss.get("data/ZA5240_v2-0-0.sav")
   4 ## y <- stata.get("data/ZA5240_v2-0-0.dta")
   5 
   6 ## the column V417 contains the net income,
   7 ## calculate the mean using the mean() function!
   8 ## What is the problem?
   9 
  10 
  11 
  12 
  13 
  14 
  15 
  16 ## summarize the net income using summary(), quantiles() and fivenum()!
  17 
  18 
  19 
  20 
  21 ## make a boxplot by using the following syntax!
  22 require(ggplot2)
  23 ggplot(x, aes(x=V86, y=V417)) +
  24     geom_boxplot()
  25 
  26 
  27 ggplot(x, aes(x=V86, y=V417, fill=V81)) +
  28     geom_boxplot()
  29 
  30 
  31 ######################################################################
  32 ####################### T-Test #######################################
  33 ######################################################################
  34 
  35 
  36 ## one sample
  37 set.seed(1)
  38 x <- rnorm(12)
  39 
  40 t.test(x,mu=0)
  41 
  42 t.test(x,mu=1)
  43 
  44 
  45 ## two sample Welch or Satterthwaite test 
  46 set.seed(1)
  47 x <- rnorm(12)
  48 y <- rnorm(12)
  49 g <- sample(c("A","B"),12,replace = T)
  50 t.test(x, y)
  51 t.test(x ~ g)
  52 t.test(x, y, var.equal = T)
  53 
  54 
  55 
  56 ######################################################################
  57 ####################### Exercises  ###################################
  58 ######################################################################
  59 
  60 ## do a t-test of income (V417): male against female (V81)!
  61 
  62 
  63 ## and compare the bmi (V279) in smokers and non-smokers (V272)
  64 ## and between the people with high and normal blood pressure (V242)
  65 
  66 
  67 
  68 
  69 
  70 ## bmi by smokers/ non smokers
  71 ggplot(x, aes(x=V272, y=V279)) +
  72     geom_boxplot()
  73 
  74 
  75 
  76 
  77 summary(x$V272)
  78 
  79 ## bmi by high/normal blood pressure
  80 
  81 
  82 
  83 
  84 
  85 ## Alter
  86 summary(x$V84)
  87