#Estadística descriptiva de datos de granja acuícola
#Se tienen datos de 12 estanques de
#alimentación y peso en semana 12
#12 gramos es el peso ideal

#Establecer folder de trabajo
setwd("~/Documents/R-PyE")
library(readxl)
camarones1 <- read_excel("camarones1.xlsx")
View(camarones1)

#variable de alimento
alimento <-camarones1$AlimentoDiario

#Tabla de frecuencias
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(%)
##   [268.714,274.85) 3 0.25 25.00  3  25.00
##   [274.85,280.986) 2 0.17 16.67  5  41.67
##  [280.986,287.121) 3 0.25 25.00  8  66.67
##  [287.121,293.257) 2 0.17 16.67 10  83.33
##  [293.257,299.393) 2 0.17 16.67 12 100.00
#Partir la pantalla de un arreglo 3x2
par(mfrow=c(3,2))
hist(alimento, breaks = "Sturges")
plot(dist, type = "cfh") #Histograma de frecuencias acumuladas
plot(dist, type = "cfp") #Poligono de frecuencias acumuladas
plot(dist, type = "fh")  #Histograma de frecuencias absolutas
plot(dist, type = "fh", col = "green") #Histogramas de frecuencuas absolutas con color

#f      = Frecuencia absoluta
#rf     = Frecuencia relativa
#rf(%)  = Frecuencia relativa porcentual
#cf     = Frecuencia acumulada
#cf(%)  = Frecuencia acumulada porcentual

sort(alimento)
##  [1] 271.4286 271.4286 271.4286 275.0000 278.5714 282.1429 282.1429 285.7143
##  [9] 289.2857 292.8571 296.4286 296.4286
#Medidas de tendencia central
mean(alimento) #Promedio, Media, Media Aritmética
## [1] 282.7381
median(alimento) #Mediana
## [1] 282.1429
summary(alimento)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   271.4   274.1   282.1   282.7   290.2   296.4
#Gráfico de caja y bigote
boxplot(alimento)

#Medidas de dispersión
#Varianza, Distancia media entre la media y todos los valores
var(alimento)
## [1] 90.05875
#La Varianza está elevada al cuadrado
#La Raíz Cuadrada de la Varianza es la Desviación Estándar
#Desviación estándar
sd(alimento)
## [1] 9.489929
#Es el grado de deviación de los datos con respecto a su media

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

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