Description of the data, data_quiz5

Q1 Import data

Hint: The data is posted in Moodle. Look for data_quiz5.csv under the Data Files section.

myRegressionData <- read.csv("data_quiz5.csv")

Q2 Review data

Hint: Use head() to display the first six rows.

head(myRegressionData, 6)
##     country continent lifeExp      pop gdpPercap
## 1   Albania    Europe  76.423  3600523  5937.030
## 2   Algeria    Africa  72.301 33333216  6223.367
## 3 Argentina  Americas  75.320 40301927 12779.380
## 4 Australia   Oceania  81.235 20434176 34435.367
## 5   Austria    Europe  79.829  8199783 36126.493
## 6   Bahrain      Asia  75.635   708573 29796.048

Q3 Visualize data

Hint: Create a scatter plot to examine the relationship between GDP per capita (mapped to y-axis) and life expectancy (mapped to x-axis).

library(tidyverse)

ggplot(myRegressionData,
       aes(x = lifeExp,
           y = gdpPercap))+
  geom_point() +
  geom_smooth(method="lm")

Q4 Build a regression model to predict GDP per capita using life expectancy.

gdpPercap_lm <- lm(gdpPercap ~ lifeExp,
                   data = myRegressionData)
summary(gdpPercap_lm)
## 
## Call:
## lm(formula = gdpPercap ~ lifeExp, data = myRegressionData)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -17319.8  -4512.4    -63.2   3443.1  24014.4 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -215340.5    18057.2  -11.93   <2e-16 ***
## lifeExp        3075.6      237.6   12.94   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 7578 on 81 degrees of freedom
## Multiple R-squared:  0.6741, Adjusted R-squared:  0.6701 
## F-statistic: 167.5 on 1 and 81 DF,  p-value: < 2.2e-16
options(scipen = 999)

Q5 Is the coefficient of life expectancy statistically significant at 5%?

The coefficient of lifeExp is statistically significant at 5% becasue when you look at the pvalue it is less that 5%.

Q6 Interpret the coefficient of life expectancy.

Hint: Discuss both its sign and magnitude.

The coefficient of life Exp is every time the lifexp goes up by 1 year the GDP increases by 3076 US Dollars. You can see a positive relationship between LifeExp and GDP.

Q7 Your friend suggests that the more populous a country, the higher its standard living (GDP per capita) is. Create a new model below by adding an additional predictor to the regression model above to test this hypothesis. Is the new variable statistically significant? What would you say to your friend regarding his/her claim?

Hint: Make your argument using the relevant test results, such as p-value.

gdpPercap_lm <- lm(gdpPercap ~ lifeExp + pop,
                   data = myRegressionData)
summary(gdpPercap_lm)
## 
## Call:
## lm(formula = gdpPercap ~ lifeExp + pop, data = myRegressionData)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -17337  -4536    -82   3463  23993 
## 
## Coefficients:
##                       Estimate         Std. Error t value
## (Intercept) -215081.8386328746   18328.1060846602 -11.735
## lifeExp        3072.6060571578     240.7532777810  12.762
## pop              -0.0000006064       0.0000056600  -0.107
##                        Pr(>|t|)    
## (Intercept) <0.0000000000000002 ***
## lifeExp     <0.0000000000000002 ***
## pop                       0.915    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 7625 on 80 degrees of freedom
## Multiple R-squared:  0.6741, Adjusted R-squared:  0.666 
## F-statistic: 82.75 on 2 and 80 DF,  p-value: < 0.00000000000000022
options(scipen = 999)

I would have to say to my friend that the claim is unreasonable because the pvalue for popuoation is not statistically significant to effecting lifeExp and GDP Percap.

Q8 Hide the messages, but display the code and its results on the webpage.

Q9 Display the title and your name correctly at the top of the webpage.

Q10 Use the correct slug.