#Estadistica descriptiva de una granja de camarones
#se tienen datos de 12 estanques durante 12 semanas
#aquí presentaré los datos de la sexta semana
setwd("~/PYE VERANO")
library(readxl)
camarones_semana6 <- read_excel("camarones_semana6.xlsx")
View(camarones_semana6)

#variable de alimento
alimento <- camarones_semana6$AlimentoDiario

#tabla de frecuencias de alimento
library(fdth)
## 
## Attaching package: 'fdth'
## The following objects are masked from 'package:stats':
## 
##     sd, var
dist <- fdt(alimento, breaks = "Sturges")
dist
##       Class limits  f   rf rf(%) cf  cf(%)
##  [125.871,127.967) 10 0.83 83.33 10  83.33
##  [127.967,130.063)  0 0.00  0.00 10  83.33
##  [130.063,132.159)  0 0.00  0.00 10  83.33
##  [132.159,134.254)  0 0.00  0.00 10  83.33
##   [134.254,136.35)  2 0.17 16.67 12 100.00
#partir la pantalla en un arreglo de 3x2
par(mfrow=c(3,2))
hist(alimento, breaks = "Sturges")
plot(dist, type="cfh")
plot(dist, type="cfp")
plot(dist, type="fh")
plot(dist, type="fh", col="blue")

#medidas de tendencia central
mean(alimento) #promedio, media
## [1] 128.4524
median(alimento) #mediana
## [1] 127.1429
summary(alimento) #da todos los valores automaticamente
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   127.1   127.1   127.1   128.5   127.1   135.0
#gráfico de caja y bigote
boxplot(alimento)

#medidas de dispersion
#varianza: distancia
var(alimento)
## [1] 9.353742
#desviacion estandar
sd(alimento)
## [1] 3.058389
#gráfica de dispersion
library(ggplot2)

ggplot(data=camarones_semana6)+
  geom_point(mapping=aes(x=AlimentoDiario, y=PesoActual))