For this chapter, you should do the simple linear regressions from the examples. For each, you should make sure you get the same results as in the example and understand the discussion of the example in the textbook. For this chapter, you should also use ggplot to create a scatter plot of the x and y variables used in the regression, and include a linear regression line. So I suggest doing this chapter after you complete the DataCamp course and your notes on [Intro to Data Visualization with ggplot2].
These are the packages you will need:
library(wooldridge)
library(ggplot2)
library(pander)For displaying the results of a single regression, you should use the pander package as shown in example 2.3 below. (Later in 380 you will be comparing several regression models side-by-side and will be using the stargazer package instead, but for a single regression model this works well).
I completed the first example ([Example 2.3: CEO Salary and Return on Equity]) for you to demonstrate what you need to do. For the subsequent examples you should fill in the code yourself. I've provided the template for you and included the code to load the data. I then added a comment that says "YOUR CODE GOES HERE" wherever you are supposed to add code. You should use the example I did for you as a guide. What you do later in 380 will not be as simple as copy/pasting code I give you and changing a few variable names, for for this chapter as you're first learning what to do, it's exactly that easy. Don't overthink what you're being asked to do.
That said, you are of course welcome to add more if it will help you. You can add more regressions. You can explore the changes of units I talk about in LN2.7. You can experiment with log transformations shown in the later parts of chapter 2 (we'll talk about these later with chapter 3). But the only things you're required to do for the BP are the 2 examples I left for you below.
Optionally, you can add notes here on chapter 2's (and LN2's) content.
Load the data from wooldridge package, estimate the regression, and display the results
# Load data
data(ceosal1)
# Estimate regression model
ex2.3 <- lm(salary ~ roe, data=ceosal1)
# Display model results
pander(summary(ex2.3)) Estimate Std. Error t value Pr(>|t|)
(Intercept) 963.2 213.2 4.517 1.053e-05
**roe** 18.5 11.12 1.663 0.09777
Observations Residual Std. Error
209 1367 0.01319 0.008421
Table: Fitting linear model: salary ~ roe
Display a scatter plot with regression line corresponding to this model
ggplot(data=ceosal1, aes(x=roe, y=salary)) + geom_point() + geom_smooth(method = "lm", se=FALSE)## `geom_smooth()` using formula 'y ~ x'
Load the data from wooldridge package, estimate the regression, and display the results
# Load data
data(wage1)
# Estimate regression model
ex2.4 <- lm(wage ~ educ, data=wage1)
# Display model results
pander(summary(ex2.4)) Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.9049 0.685 -1.321 0.1871
**educ** 0.5414 0.05325 10.17 2.783e-22
Observations Residual Std. Error
526 3.378 0.1648 0.1632
Table: Fitting linear model: wage ~ educ
Display a scatter plot with regression line corresponding to this model
ggplot(data=wage1, aes(x=educ, y=wage)) + geom_point() + geom_smooth(method = "lm", se=FALSE)## `geom_smooth()` using formula 'y ~ x'
Load the data from wooldridge package, estimate the regression, and display the results
# Load data
data(vote1)
# Estimate regression model
ex2.5 <- lm(voteA ~ shareA, data=vote1)
# Display model results
pander(summary(ex2.5)) Estimate Std. Error t value Pr(>|t|)
(Intercept) 26.81 0.8872 30.22 1.729e-70
Observations Residual Std. Error
173 6.385 0.8561 0.8553
Table: Fitting linear model: voteA ~ shareA
Display a scatter plot with regression line corresponding to this model
ggplot(data=vote1, aes(x=shareA, y=voteA)) + geom_point() + geom_smooth(method = "lm", se=FALSE)## `geom_smooth()` using formula 'y ~ x'
This example builds off of [Example 2.3: CEO Salary and Return on Equity] above. You don't need to do anything for these, but you should look at them to make sure you understand them.
First, you should understand how to calculate the OLS fitted values, the
After estimating a model, you can use the coef() function to get the values of
We estimated the model above and stored it into the variable ex2.3. We can access coef(ex2.3)["(Intercept)"] or coef(ex2.3)[1]. We can access the first (and only) roe, and coef(ex2.3)["roe"] or coef(ex2.3)[1].
ceosal1$salaryHat <- coef(ex2.3)["(Intercept)"] + coef(ex2.3)["roe"] * ceosal1$roeThe second way is by using the fitted() function.
ceosal1$salaryHat_fitted <- fitted(ex2.3)You can check that these two are the same by subtracting them and making sure they are all the same.
ceosal1$dif <- ceosal1$salaryHat - ceosal1$salaryHat_fitted
summary(ceosal1$dif)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## -2.046e-12 0.000e+00 0.000e+00 2.557e-14 2.274e-13 6.821e-13
Note that the min is -0.000000000002046363 and the max is 0.000000000000682121. Algebra with decimal numbers often results in small rounding errors, which is why these are 0 exactly, but they are effectively 0, demonstrating that the two methods of calculating the fitted
The second reason I'm including this example is to demonstrate plotting the fitted
ggplot(data=ceosal1, aes(x=roe, y=salary)) +
geom_point(aes(color="data",shape="data")) +
geom_point(aes(x=roe, y=salaryHat,color="predicted (yHat)",shape="predicted (yHat)")) +
geom_smooth(method = "lm",color="blue", se=FALSE,aes(linetype="OLS")) +
scale_color_manual(values = c("black", "red"),name="Values") +
scale_shape_manual(values = c(1,4),name="Values") +
scale_linetype_manual(values = c("solid","solid"), name="Regression Line")## `geom_smooth()` using formula 'y ~ x'
