You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Copy file name to clipboardExpand all lines: inst/tutorials/statistical_models_ws/statistical_models_ws.Rmd
+10-10Lines changed: 10 additions & 10 deletions
Original file line number
Diff line number
Diff line change
@@ -622,12 +622,12 @@ Unlike the fruit fly data set, no pre-manipulation is needed so let's view the d
622
622
623
623
```{r gapminder-package, exercise = TRUE}
624
624
# load and check the data
625
-
library(<package>) # not installed on this machine
625
+
library(<package>)
626
626
<function>(gapminder)
627
627
```
628
628
629
629
```{r gapminder-package-solution}
630
-
library(gapminder) # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder
630
+
library(gapminder) # Data from Gapminder
631
631
head(gapminder)
632
632
```
633
633
@@ -650,7 +650,7 @@ gapminder %>%
650
650
```
651
651
652
652
```{r setup-gapminder-plot, include = FALSE}
653
-
library(gapminder) # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder
653
+
library(gapminder) # Data from Gapminder
654
654
```
655
655
656
656
Although there is very high variance, we do see a certain trend with mean life expectancy increasing over time. Similarly, we can naively hypothesize that life expectancy is higher where the per-capita GDP is higher. In R, this is represented with the formula: `lifeExp ~ gdpPercap`.
@@ -673,7 +673,7 @@ gapminder %>%
673
673
```
674
674
675
675
```{r setup-gapminder-plot2, include = FALSE}
676
-
library(gapminder) # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder
676
+
library(gapminder) # Data from Gapminder
677
677
```
678
678
679
679
In the cases above, we plotted a relationship between a continuous *dependent variable* and a continuous *explanatory* variable.
library(gapminder) # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder
716
+
library(gapminder) # Data from Gapminder
717
717
```
718
718
719
719
Now, use the function `summary()` to see all the relevant results.
@@ -727,7 +727,7 @@ summary(lifeExp_model1)
727
727
```
728
728
729
729
```{r setup-gapminder-model1-sum, include = FALSE}
730
-
library(gapminder) # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder
730
+
library(gapminder) # Data from Gapminder
731
731
lifeExp_model1 <- lm(lifeExp ~ year,
732
732
data = gapminder)
733
733
```
@@ -763,7 +763,7 @@ gapminder %>%
763
763
```
764
764
765
765
```{r setup-gapminder-line1, include = FALSE}
766
-
library(gapminder) # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder
766
+
library(gapminder) # Data from Gapminder
767
767
```
768
768
769
769
### Residuals
@@ -782,7 +782,7 @@ plot(lifeExp_model1)
782
782
```
783
783
784
784
```{r setup-gapminder-res, include = FALSE}
785
-
library(gapminder) # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder
785
+
library(gapminder) # Data from Gapminder
786
786
lifeExp_model1 <- lm(lifeExp ~ year,
787
787
data = gapminder)
788
788
```
@@ -875,7 +875,7 @@ par(mfrow=c(2,2))
875
875
```
876
876
877
877
```{r setup-ex2, include = FALSE}
878
-
library(gapminder) # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder
878
+
library(gapminder) # Data from Gapminder
879
879
```
880
880
881
881
```{r ex2-hint-1}
@@ -963,7 +963,7 @@ plot(lifeExp_model2)
963
963
```
964
964
965
965
```{r setup-mlr-one-res, include = FALSE}
966
-
library(gapminder) # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder # Data from Gapminder
0 commit comments