Simple Linear Regression is a statistical method used to predict one dependent variable based on a single independent variable.
2025-03-17
Simple Linear Regression is a statistical method used to predict one dependent variable based on a single independent variable.
The simple linear regression model is represented as:
\[ y = \beta_0 + \beta_1 x + \varepsilon \]
The coefficient of determination \(R^2\) is calculated using the formula:
\[ R^2 = 1 - \frac{SSR}{SST} \]
We will use the mtcars dataset:
mpg (Miles Per Gallon): Fuel efficiencywt (Weight): Vehicle weighthead(mtcars)
## mpg cyl disp hp drat wt qsec vs am gear carb ## Mazda RX4 21.0 6 160 110 3.90 2.620 16.46 0 1 4 4 ## Mazda RX4 Wag 21.0 6 160 110 3.90 2.875 17.02 0 1 4 4 ## Datsun 710 22.8 4 108 93 3.85 2.320 18.61 1 1 4 1 ## Hornet 4 Drive 21.4 6 258 110 3.08 3.215 19.44 1 0 3 1 ## Hornet Sportabout 18.7 8 360 175 3.15 3.440 17.02 0 0 3 2 ## Valiant 18.1 6 225 105 2.76 3.460 20.22 1 0 3 1
library(ggplot2)
ggplot(mtcars, aes(x = wt, y = mpg)) +
geom_point(color = "darkblue") +
labs(title = "Relationship between Vehicle Weight and MPG",
x = "Weight (wt)",
y = "Miles Per Gallon (mpg)") +
theme_minimal()
model <- lm(mpg ~ wt, data = mtcars) summary(model)
## ## Call: ## lm(formula = mpg ~ wt, data = mtcars) ## ## Residuals: ## Min 1Q Median 3Q Max ## -4.5432 -2.3647 -0.1252 1.4096 6.8727 ## ## Coefficients: ## Estimate Std. Error t value Pr(>|t|) ## (Intercept) 37.2851 1.8776 19.858 < 2e-16 *** ## wt -5.3445 0.5591 -9.559 1.29e-10 *** ## --- ## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 ## ## Residual standard error: 3.046 on 30 degrees of freedom ## Multiple R-squared: 0.7528, Adjusted R-squared: 0.7446 ## F-statistic: 91.38 on 1 and 30 DF, p-value: 1.294e-10
ggplot(mtcars, aes(x = wt, y = mpg)) +
geom_point(color = "darkblue") +
geom_smooth(method = "lm", se = TRUE, color = "red") +
labs(title = "Regression Line",
x = "Weight (wt)",
y = "Miles Per Gallon (mpg)")
## `geom_smooth()` using formula = 'y ~ x'
library(plotly)
x <- mtcars$wt
y <- mtcars$hp
z <- mtcars$mpg
plot_ly(x = ~x, y = ~y, z = ~z, type = "scatter3d", mode = "markers") %>%
layout(scene = list(xaxis = list(title = 'Weight (wt)'),
yaxis = list(title = 'Horsepower (hp)'),
zaxis = list(title = 'Miles Per Gallon (mpg)')))