estaparcial2

setwd("C:/Users/fiore/OneDrive/Escritorio/estaparcial")

library(rio)

## Warning: package 'rio' was built under R version 4.3.3

library(tidyverse)

## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.1.3     ✔ readr     2.1.4
## ✔ forcats   1.0.0     ✔ stringr   1.5.0
## ✔ ggplot2   3.4.3     ✔ tibble    3.2.1
## ✔ lubridate 1.9.2     ✔ tidyr     1.3.0
## ✔ purrr     1.0.2     
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors

library(haven)

## Warning: package 'haven' was built under R version 4.3.3

basedatos = import("DataPeru.xlsx")

str(basedatos)

## 'data.frame':    25 obs. of  8 variables:
##  $ DEPARTAMENTO       : chr  "AMAZONAS" "ÁNCASH" "APURÍMAC" "AREQUIPA" ...
##  $ UBIGEO             : chr  "010000" "020000" "030000" "040000" ...
##  $ buenEstado         : num  18.6 13.9 8.7 27.4 17 18 33.8 11.9 10.1 15.6 ...
##  $ contribuyentesSunat: num  75035 302906 103981 585628 151191 ...
##  $ peaOcupada         : num  130019 387976 140341 645001 235857 ...
##  $ pobUrbana          : num  205976 806065 243354 1383694 444473 ...
##  $ PobRural           : num  211389 333050 180905 76739 206467 ...
##  $ pobTotal           : num  417365 1139115 424259 1460433 650940 ...

#EJERCICIO 1
# Modelado de la relación
modelo <- lm(buenEstado ~ contribuyentesSunat + peaOcupada, data = basedatos)

# Interpretación del modelo
summary(modelo)

## 
## Call:
## lm(formula = buenEstado ~ contribuyentesSunat + peaOcupada, data = basedatos)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -10.589  -3.966  -1.347   1.907  21.518 
## 
## Coefficients:
##                       Estimate Std. Error t value Pr(>|t|)    
## (Intercept)          1.865e+01  2.694e+00   6.922 5.98e-07 ***
## contribuyentesSunat  1.786e-05  2.060e-05   0.867    0.395    
## peaOcupada          -1.596e-05  2.241e-05  -0.712    0.484    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 7.925 on 22 degrees of freedom
## Multiple R-squared:  0.1561, Adjusted R-squared:  0.07939 
## F-statistic: 2.035 on 2 and 22 DF,  p-value: 0.1546

# Prueba de significancia
anova(modelo)

## Analysis of Variance Table
## 
## Response: buenEstado
##                     Df  Sum Sq Mean Sq F value  Pr(>F)  
## contribuyentesSunat  1  223.76 223.758  3.5627 0.07235 .
## peaOcupada           1   31.84  31.838  0.5069 0.48395  
## Residuals           22 1381.71  62.805                  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

#EJERCICIO 2
# Modelando
modelo <- lm(peaOcupada ~ contribuyentesSunat + buenEstado, data = basedatos)

# Interpretacion
summary(modelo)

## 
## Call:
## lm(formula = peaOcupada ~ contribuyentesSunat + buenEstado, data = basedatos)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -91867 -58573 -11166  46174 155851 
## 
## Coefficients:
##                       Estimate Std. Error t value Pr(>|t|)    
## (Intercept)          1.155e+05  3.787e+04   3.049  0.00588 ** 
## contribuyentesSunat  9.206e-01  1.741e-02  52.872  < 2e-16 ***
## buenEstado          -1.412e+03  1.983e+03  -0.712  0.48395    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 74540 on 22 degrees of freedom
## Multiple R-squared:  0.9932, Adjusted R-squared:  0.9926 
## F-statistic:  1603 on 2 and 22 DF,  p-value: < 2.2e-16

# Significancia
anova(modelo)

## Analysis of Variance Table
## 
## Response: peaOcupada
##                     Df     Sum Sq    Mean Sq   F value Pr(>F)    
## contribuyentesSunat  1 1.7813e+13 1.7813e+13 3205.7524 <2e-16 ***
## buenEstado           1 2.8169e+09 2.8169e+09    0.5069  0.484    
## Residuals           22 1.2225e+11 5.5567e+09                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1