Desempleo_Homicidios

library(readxl)
caso <- read_excel("C:/Users/User/Desktop/caso.xlsx")

attach(caso)

#a)

plot(desempleo,homicidios)

#b)

cor(desempleo, homicidios)

## [1] 0.9608183

#c)

modelo= lm(homicidios~desempleo)
summary(modelo)

## 
## Call:
## lm(formula = homicidios ~ desempleo)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -19.335 -11.928  -4.618   6.006  62.193 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -628.936     35.846  -17.55   <2e-16 ***
## desempleo     63.751      2.983   21.37   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 18.06 on 38 degrees of freedom
## Multiple R-squared:  0.9232, Adjusted R-squared:  0.9212 
## F-statistic: 456.6 on 1 and 38 DF,  p-value: < 2.2e-16

#d)

par(mfrow=c(2,2))
plot(modelo)

#e)

exp(predict(modelo,newdata = list(desempleo=11)))

##            1 
## 2.580847e+31