#Estadística descriptiva de datos de granja acuicola
#Se tienen datos de 12 estanque 
#de alimentos y peso en semana 8
#12 gramos es el peso ideal

setwd("~/VERANOPYE")
library(readxl)
Camarones8 <- read_excel("Camarones8.xlsx")
View(Camarones8)

#Variable de alimento
alimento <- Camarones8$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(%)
##  [142.136,146.606) 1 0.08  8.33  1   8.33
##  [146.606,151.076) 3 0.25 25.00  4  33.33
##  [151.076,155.546) 3 0.25 25.00  7  58.33
##  [155.546,160.016) 3 0.25 25.00 10  83.33
##  [160.016,164.486) 2 0.17 16.67 12 100.00
#Partir la pantalla ne un arreglo de 3x2
par(mfrow=c(3,2))
hist(alimento, breaks="Sturges")
plot(dist, type="cfh") #Histograma de frecuencias acumuladas
plot(dist, type="cfp") #Polígono de frecuencias acumuladas
plot(dist, type="fh") #Histograma de frecuencias absolutas
plot(dist, type="fh", col="green") #Histograma de frecuencias absolutas, verde


sort(alimento)
##  [1] 143.5714 149.2857 150.0000 150.7143 151.4286 152.8571 155.0000 156.4286
##  [9] 157.8571 159.2857 161.4286 162.8571
#Medidas de tendencia central
mean(alimento) #promedio, media, media aritmética
## [1] 154.2262
median(alimento) #mediana
## [1] 153.9286
summary(alimento) 
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   143.6   150.5   153.9   154.2   158.2   162.9
#Grafico de caja y bigote
boxplot(alimento)

#Medidas de dispersión
#Varianza, distancia entre media entre la media y todos los valores
var(alimento)
## [1] 31.67517
#La varianza está elevada al cuadrado

#La raiz cuadrada de la varianza es la:
#Desviación estándar
sd(alimento)
## [1] 5.62807
#El grado de desviación de los datos con respecto a su media

#Gráfica de dispersión
#del alimento ndiario de la semana 12 en gramos
#y el peso de los camarones en gramos
library(ggplot2)

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

#CONCLUSIÓN DE HISTOGRAMAS
#Al 75% de los tanques se les proporcionó entre 146 y 160 gramos de alimento, estando cerca de la media que fue de 154.
#La moda y mediana se mantuvieron cerca entre sí, 154.22 (media) y 153.92 (mediana).

#CONCLUSIÓN GRÁFICA DE DISPERSIÓN:
#Los puntos se mantuvieron mayormente concentrados en el centro, como ya se dijo, la mayoría de los tanques (75%) se mantuvieron cerca de la media.
#Se alcanza a observar que 9 de los puntos se encuentran entre 146 y 160 gramos de alimentación.
#Se puede inferir que en la semana 8, 9 de los tanques estuvieron siendo alimentados de la misma forma pero obteniendo pesos diferentes y no todos los camarones de los distintos tanques no se desarrollaron de la misma forma.