Andreu tarea 1 R y Estadistica descriptiva

Graficas, cajas y bigotes

qqnorm(AREVALO_ALCANTARA_ANDREU$nit)
qqline(AREVALO_ALCANTARA_ANDREU$nit)

boxplot(AREVALO_ALCANTARA_ANDREU$nit)

boxplot.stats(AREVALO_ALCANTARA_ANDREU$nit)

## $stats
## [1]  9.50 21.15 26.80 31.65 39.10
## 
## $n
## [1] 44
## 
## $conf
## [1] 24.29896 29.30104
## 
## $out
## numeric(0)

AREVALO_ALCANTARA_ANDREU$normit <- (AREVALO_ALCANTARA_ANDREU$nit - mean(AREVALO_ALCANTARA_ANDREU$nit) / sd(AREVALO_ALCANTARA_ANDREU$nit))
plot(AREVALO_ALCANTARA_ANDREU$normit, AREVALO_ALCANTARA_ANDREU$nit)

valores <- c(2, 13, 19, 1, 3, 17, 0, 20, 15)
matriz <- matrix(valores, nrow = 3, ncol = 3, byrow = TRUE)
print(matriz)

##      [,1] [,2] [,3]
## [1,]    2   13   19
## [2,]    1    3   17
## [3,]    0   20   15

diag(matriz)

## [1]  2  3 15

det(matriz)

## [1] -405

solve(matriz)

##             [,1]        [,2]        [,3]
## [1,]  0.72839506 -0.45679012 -0.40493827
## [2,]  0.03703704 -0.07407407  0.03703704
## [3,] -0.04938272  0.09876543  0.01728395

sum(matriz[3,]*3)

## [1] 105

indices <- c(19, 13, 1, 29, 25, 22, 8, 7, 39, 16, 34, 15, 41, 18, 38, 17, 40, 44, 14, 3)
conjunto <- AREVALO_ALCANTARA_ANDREU[indices,]
summary(conjunto)

##      Fuente           nit            normit      
##  Min.   : 1.00   Min.   :10.80   Min.   : 6.884  
##  1st Qu.:13.75   1st Qu.:23.40   1st Qu.:19.484  
##  Median :18.50   Median :26.85   Median :22.934  
##  Mean   :22.15   Mean   :26.86   Mean   :22.939  
##  3rd Qu.:35.00   3rd Qu.:32.73   3rd Qu.:28.809  
##  Max.   :44.00   Max.   :39.10   Max.   :35.184

indices2 <- c(19, 13, 1, 29, 25, 22, 8, 7, 39, 16, 34,
15, 41, 18, 38, 17, 40, 44, 14, 3)
AREVALO_ALCANTARA_ANDREU$nit2 <- AREVALO_ALCANTARA_ANDREU$nit
AREVALO_ALCANTARA_ANDREU$nit2[indices2] <- AREVALO_ALCANTARA_ANDREU$nit2[indices2] + 5
AREVALO_ALCANTARA_ANDREU$type <- "B"
AREVALO_ALCANTARA_ANDREU$type[indices2] <- "A"
AREVALO_ALCANTARA_ANDREU$type <- as.factor(AREVALO_ALCANTARA_ANDREU$type)
boxplot(nit2~type, data = AREVALO_ALCANTARA_ANDREU)

aggregate(nit~type, data = AREVALO_ALCANTARA_ANDREU, summary)

##   type nit.Min. nit.1st Qu. nit.Median nit.Mean nit.3rd Qu. nit.Max.
## 1    A 10.80000    23.40000   26.85000 26.85500    32.72500 39.10000
## 2    B  9.50000    20.27500   26.70000 25.84167    30.50000 35.50000

aggregate(nit2~type, data = AREVALO_ALCANTARA_ANDREU, summary)

##   type nit2.Min. nit2.1st Qu. nit2.Median nit2.Mean nit2.3rd Qu. nit2.Max.
## 1    A  15.80000     28.40000    31.85000  31.85500     37.72500  44.10000
## 2    B   9.50000     20.27500    26.70000  25.84167     30.50000  35.50000