A11lab

library(ggplot2)

## Warning: package 'ggplot2' was built under R version 4.0.5

library(dplyr)

## Warning: package 'dplyr' was built under R version 4.0.5

## 
## Attaching package: 'dplyr'

## The following objects are masked from 'package:stats':
## 
##     filter, lag

## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union

library(Metrics)

## Warning: package 'Metrics' was built under R version 4.0.5

library(readr)
burnwv <- read.csv("C:/Users/jmhp2/Downloads/single_regression/single_regression/burn_wv3.csv")

burnwv$R1 <- ((burnwv$Dc.B8-burnwv$Dc.B16)/(burnwv$Dc.B8+burnwv$Dc.B16))-((burnwv$Ap.B8-burnwv$Ap.B16)/(burnwv$Ap.B8+burnwv$Ap.B16))

R1lm <- lm(Total~R1, data = burnwv)
R1rmse <- rmse(burnwv$Total, R1lm$fitted.values)
R1r <- summary(R1lm)$r.squared

burnwv$R2 <- ((burnwv$Dc.B7-burnwv$Dc.B12)/(burnwv$Dc.B7+burnwv$Dc.B12))-((burnwv$Ap.B7-burnwv$Ap.B12)/(burnwv$Ap.B7+burnwv$Ap.B12))

R2lm <- lm(Total~R2, data = burnwv)
R2rmse <- rmse(burnwv$Total, R2lm$fitted.values)
R2r <- summary(R2lm)$r.squared

burnwv$R3 <- ((burnwv$Dc.B7-burnwv$Dc.B12)/(burnwv$Dc.B8+burnwv$Dc.B16))-((burnwv$Ap.B7-burnwv$Ap.B12)/(burnwv$Ap.B8+burnwv$Ap.B16))

R3lm <- lm(Total~R3, data = burnwv)
R3rmse <- rmse(burnwv$Total, R3lm$fitted.values)
R3r <- summary(R3lm)$r.squared

burnwv$R4 <- ((burnwv$Dc.B8-burnwv$Dc.B16)/(burnwv$Dc.B7+burnwv$Dc.B12))-((burnwv$Ap.B8-burnwv$Ap.B16)/(burnwv$Ap.B7+burnwv$Ap.B12))

R4lm <- lm(Total~R4, data = burnwv)
R4rmse <- rmse(burnwv$Total, R4lm$fitted.values)
R4r <- summary(R4lm)$r.squared

ggplot(burnwv)+
geom_point(aes(x=R1, y=R1rmse), size=1, color = "blue")

ggplot(burnwv)+
geom_point(aes(x=R2, y=R2rmse), size=1, color = "red")

ggplot(burnwv)+
geom_point(aes(x=R3, y=R3rmse), size=1, color = "purple")

ggplot(burnwv)+
geom_point(aes(x=R4, y=R4rmse), size=1, color = "orange")

#q1: R3 is the highest r^2 and the lowest rmse of all four.

ggplot(R3lm, aes(x=R3rmse, y=R3 ))+
geom_point(size = 1, color = "green")

A11lab_