-
Notifications
You must be signed in to change notification settings - Fork 0
/
coeficiente_determinacao.R
76 lines (39 loc) · 2.18 KB
/
coeficiente_determinacao.R
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
## ----settings, echo=FALSE-----------------------------------------------------
library(knitr)
## ----importa-dadinhos---------------------------------------------------------
dadinhos <- read.csv("https://raw.githubusercontent.com/piLaboratory/BIE-0320/main/data/dadinho.csv")
## ----mostra-dadinhos, eval=FALSE----------------------------------------------
## dadinhos
## ---- echo=FALSE--------------------------------------------------------------
dadinhos$dif <- dadinhos$Y - mean(dadinhos$Y)
## ---- eval=FALSE--------------------------------------------------------------
## dadinhos$dif <- dadinhos$Y - mean(dadinhos$Y)
## dadinhos
## -----------------------------------------------------------------------------
V.total <- sum(dadinhos$dif^2)
## ---- eval=FALSE--------------------------------------------------------------
## V.total
## -----------------------------------------------------------------------------
dadinhos.lm <- lm(Y ~ X, data=dadinhos)
## ---- echo=FALSE--------------------------------------------------------------
dadinhos.cf <- coef(dadinhos.lm)
## ---- eval=FALSE--------------------------------------------------------------
## (dadinhos.cf <- coef(dadinhos.lm))
## -----------------------------------------------------------------------------
dadinhos$Y.pred <- predict(dadinhos.lm)
## -----------------------------------------------------------------------------
dadinhos$residuo <- dadinhos$Y - dadinhos$Y.pred
## ---- echo=FALSE--------------------------------------------------------------
kable(dadinhos)
## -----------------------------------------------------------------------------
V.resid <- sum(dadinhos$residuo^2)
## ---- eval=FALSE--------------------------------------------------------------
## V.resid
## ---- echo=FALSE--------------------------------------------------------------
V.expl <- V.total - V.resid
## ---- eval=FALSE--------------------------------------------------------------
## (V.expl <- V.total - V.resid)
## ---- eval=FALSE--------------------------------------------------------------
## V.expl/V.total
## ---- eval=FALSE--------------------------------------------------------------
## summary(dadinhos.lm)$r.squared