### 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"))