Advanced-R-Programming-Course/oop_output.txt at master · Chava93/Advanced-R-Programming-Course · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
> ## Read in the data
> library(readr)
> library(magrittr)
> source("oop_code.R")
>
> data <- read_csv("data/MIE.csv")
Parsed with column specification:
cols(
  id = col_integer(),
  visit = col_integer(),
  room = col_character(),
  value = col_double(),
  timepoint = col_integer()
)
> x <- make_LD(data)
> print(class(x))
[1] "LongitudinalData"
> print(x)
Longitudinal dataset with 10 subjects
>
> ## Subject 10 doesn't exist
> out <- subject(x, 10)
> print(out)
NULL
>
> out <- subject(x, 14)
> print(out)
Subject ID: 14
>
> out <- subject(x, 54) %>% summary
> print(out)
ID: 54
  visit  bedroom       den living room    office
1     0       NA        NA    2.792601 13.255475
2     1       NA 13.450946          NA  4.533921
3     2 4.193721  3.779225          NA        NA
>
> out <- subject(x, 14) %>% summary
> print(out)
ID: 14
  visit   bedroom family  room living room
1     0  4.786592           NA     2.75000
2     1  3.401442     8.426549          NA
3     2 18.583635           NA    22.55069
>
> out <- subject(x, 44) %>% visit(0) %>% room("bedroom")
> print(out)
ID: 44
Visit: 0
Room: bedroom
>
> ## Show a summary of the pollutant values
> out <- subject(x, 44) %>% visit(0) %>% room("bedroom") %>% summary
> print(out)
ID: 44
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
    8.0    30.0    51.0    88.8    80.0   911.0
>
> out <- subject(x, 44) %>% visit(1) %>% room("living room") %>% summary
> print(out)
ID: 44
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
   2.75   14.00   24.00   41.37   37.00 1607.00