---
title: "Regresión Lineal - Energía"
author: "Ricardo Palma"
date: "22/10/2021"
output: html_document
---

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
```

# Regresión Lineal

## Energı́a Solar en Las Vegas - Potrerillos Mendoza

```{r}
eficiency <- read.table("https://themys.sid.uncu.edu.ar/rpalma/R-cran/Solar.txt",header = TRUE)

names(eficiency)

```

```{r}
head(eficiency)
```


```{r}
pairs(eficiency)

```


Mejorando la apariencia de los gráficos

```{r}
library(psych)
pairs.panels(eficiency)
multi.hist(eficiency)
cor(eficiency)

```



## Construcción de la regresión lienal

```{r}
regresion <- lm(solar ~ gas, data = eficiency)
summary(regresion)

```


## Predicciones

```{r}
nuevas.gases <- data.frame(gas = seq(30, 50))
predict(regresion, nuevas.gases)

```

## Intervalo de confinaza de la predicción


```{r}
confint(regresion)

```


## Ploteo del modelo lineal con su intervalo de confianza


```{r}
options(tidy=TRUE, width=50)
 nuevas.gases <- data.frame(gas = seq(10, 90))
 plot(eficiency$gas, eficiency$solar, xlab = "gas", ylab = "solar")
 abline(regresion)
 ic <- predict(regresion, nuevas.gases, interval = "confidence")
 lines(nuevas.gases$gas, ic[, 2], lty = 2)
 lines(nuevas.gases$gas, ic[, 3], lty = 2)
 ic <- predict(regresion, nuevas.gases, interval ="prediction")
 lines(nuevas.gases$gas, ic[, 2], lty = 2, col = "red")
 lines(nuevas.gases$gas, ic[, 3], lty = 2, col = "red")


```

Intervalos de confianza de la respuesta media: ic es una matriz con tres
columnas: la primera es la prediccion, las otras dos son los extremos del intervalo

## Analisis de Varianza

```{r}
anova(regresion)

```


## Cálculo de los residuos (residuals)

```{r}
 residuos <- rstandard(regresion)
 valores.ajustados <- fitted(regresion)
 plot(valores.ajustados, residuos)
 qqnorm(residuos)
 qqline(residuos)

```

