Newer
Older
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
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
### R code from vignette source 'Exercise_1_trees_SOLUTION.Rnw'
###################################################
### code chunk number 1: Exercise_1_trees_SOLUTION.Rnw:49-50
###################################################
load("Covid_stress_data_subs.rdata")
###################################################
### code chunk number 2: Exercise_1_trees_SOLUTION.Rnw:55-56
###################################################
summary(dat_cvid)
###################################################
### code chunk number 3: Exercise_1_trees_SOLUTION.Rnw:63-65
###################################################
#install.packages("partykit")
library("partykit")
###################################################
### code chunk number 4: Exercise_1_trees_SOLUTION.Rnw:70-71
###################################################
my_ctree <- ctree(vaccWill ~ ., data = dat_cvid)
###################################################
### code chunk number 5: Exercise_1_trees_SOLUTION.Rnw:76-77
###################################################
plot(my_ctree)
###################################################
### code chunk number 6: Exercise_1_trees_SOLUTION.Rnw:90-93
###################################################
predicted <- predict(my_ctree, type = "response")
true <- dat_cvid$vaccWill
table(true, predicted)
###################################################
### code chunk number 7: Exercise_1_trees_SOLUTION.Rnw:104-106
###################################################
#install.packages("stablelearner")
library("stablelearner")
###################################################
### code chunk number 8: Exercise_1_trees_SOLUTION.Rnw:113-115
###################################################
set.seed(270312)
stab_my_ctree <- stabletree(my_ctree)
###################################################
### code chunk number 9: Exercise_1_trees_SOLUTION.Rnw:123-124 (eval = FALSE)
###################################################
## plot(stab_my_ctree, select = c("vaccAtt", "fearCorona", "compliance", "education"))
###################################################
### code chunk number 10: Exercise_1_trees_SOLUTION.Rnw:127-131 (eval = FALSE)
###################################################
## # show variable selection proportions
## barplot(stab_my_ctree)
## # illustrate variable selections of replications
## image(stab_my_ctree)
###################################################
### code chunk number 11: Exercise_1_trees_SOLUTION.Rnw:135-137
###################################################
# illustrate variable selections of replications
barplot(stab_my_ctree)
###################################################
### code chunk number 12: Exercise_1_trees_SOLUTION.Rnw:139-141
###################################################
# illustrate variable selections of replications
image(stab_my_ctree)
###################################################
### code chunk number 13: Exercise_1_trees_SOLUTION.Rnw:151-153
###################################################
# show graphical cutpoint analysis
plot(stab_my_ctree, select = c("vaccAtt", "fearCorona", "compliance", "education", "age", "trustInstitut"))