forked from mca91/EconometricsWithR
-
Notifications
You must be signed in to change notification settings - Fork 0
/
ex7_6.html
88 lines (73 loc) · 2.52 KB
/
ex7_6.html
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
<!DOCTYPE html>
<head>
<meta charset='utf-8'/>
<script async='async' src='https://cdn.datacamp.com/dcl-react-dev.js.gz'></script>
<style>
.DCexercise .datacamp-exercise {
border: 2px solid #3D678D;
border-radius: 10px 10px 10px 10px !important;
}
*[class*="lm_"], .ace_gutter, textarea[class*="ace_"], .ace_scroller {
background-color: #F0EFF0 !important;
}
div[class*="Editor-module__editor"], div[class*="dcl__Footer"] {
background-color: #3D678D !important;
}
*[class*="lm_"] {
border-radius: 10px 10px 0px 0px !important;
}
div[class*="dcl__Footer"] {
border-radius: 0px 0px 10px 10px !important;
}
.lm_content {
border-radius: 0px !important;
}
.lm_splitter {
background-color: #3D678D !important;
width: 3px !important;
}
.lm_drag_handle {
background: transparent !important;
}
div[class*="dcl__index-module"] {
outline-color: transparent !important;
border-radius: 10px !important;
}
button[class*="secondary-light"] {
background-color: #FF0000;
}
</style>
</head>
<body>
<div data-datacamp-exercise data-lang='r' data-height =320><code data-type="pre-exercise-code">
library(AER)
library(MASS)
mod <- lm(medv ~ lstat + crim + age, data = Boston)
</code>
<code data-type="sample-code">
# construct a 99% confidence set for the coefficients of `crim` and `lstat`
# conduct the corresponding F-test
</code>
<code data-type="solution">
# construct a 99% confidence set for the coefficients of `crim` and `lstat`
confidenceEllipse(mod, which.coef = c("crim", "lstat"), levels = 0.99)
# conduct the corresponding F-test
linearHypothesis(mod, c("crim = 0", "lstat = 0"))
</code>
<code data-type="sct">
test_or({
fun <- ex() %>% check_function("confidenceEllipse")
fun %>% check_arg("model") %>% check_equal()
fun %>% check_arg("which.coef") %>% check_equal()
fun %>% check_arg("levels") %>% check_equal()
}, {
fun <- ex() %>% override_solution("confidenceEllipse(mod, which.coef = c(\"lstat\", \"crim\"), levels = 0.99)") %>% check_function("confidenceEllipse")
fun %>% check_arg("model") %>% check_equal()
fun %>% check_arg("which.coef") %>% check_equal()
fun %>% check_arg("levels") %>% check_equal()
})
test_function_result("linearHypothesis")
success_msg("Correct! Since (0,0) is not an element of the 99% confidence set, the null hypothesis is rejected at the 1% level.")
</code></div>
</body>
</html>