Vignette Type - in class exercise 1

Data

Read in the data from the website.

dta <- read.csv("https://stats.idre.ucla.edu/stat/data/binary.csv")

Keep only gpa and gre variables

Create new data object.

dta <- dta[, c("gpa", "gre")] 

Plot

plot(dta, xlab = "GPA", ylab = "GRE")
grid()

Analysis

options(digits = 4, show.signif.stars = FALSE)
summary(m0 <- lm(gpa ~ gre, data = dta))

## 
## Call:
## lm(formula = gpa ~ gre, data = dta)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -1.0867 -0.2244 -0.0002  0.2481  0.7618 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)
## (Intercept) 2.645898   0.091310    29.0  < 2e-16
## gre         0.001266   0.000152     8.3  1.6e-15
## 
## Residual standard error: 0.352 on 398 degrees of freedom
## Multiple R-squared:  0.148,  Adjusted R-squared:  0.146 
## F-statistic: 68.9 on 1 and 398 DF,  p-value: 1.6e-15

End