#Estadistica descriptiva de datos de granja acuicola
#Se tienen datos de 12 estanques de 
#alimentacion y peso en semana 12
#12 gramos es el peso ideal
setwd("~/PyE")
library(readxl)
camarones2 <- read_excel("camarones2.xlsx")#Esta va a ser la entrada de datos de los camarones
View(camarones2)

#variable de alimento
alimento <- camarones2$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(%)
##  [175.371,183.583) 1 0.08  8.33  1   8.33
##  [183.583,191.794) 3 0.25 25.00  4  33.33
##  [191.794,200.006) 4 0.33 33.33  8  66.67
##  [200.006,208.217) 3 0.25 25.00 11  91.67
##  [208.217,216.429) 1 0.08  8.33 12 100.00
#partir la pantalla en un arreglo de 3x2

par(mfrow = c(3, 2))
hist(alimento, breaks = "Sturges")
#Esta es la cantidad de camarones que ingieren cierta cantidad de alimento de forma ordenada 

plot(dist, type = "cfh")#histograma de frecuencia acumulada
#Esta tambien muestra la cantidad de camarones que ingieren cierta cantidad de alimento pero se le suma lo del parametro(o rango) anterior

plot(dist, type = "cfp")#poligono de frecuenciaas acumuladas
#Esta representa lo mismo que la anterior pero con un poligono

plot(dist, type = "fh")#Histogrma de frecuencias absolutas

plot(dist, type = "fh", col = "green")#histograma de frecuencias absolutas a color

#medidas de tendencia central

mean(alimento)#promedio, media, media aritmetica
## [1] 195.9524
median(alimento)#Mediana
## [1] 194.2857
summary(alimento)#resumen de una grafica de caja y bigotes(esta contiene los cuantiles la media y la mediana)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   177.1   190.9   194.3   196.0   203.0   214.3
sort(alimento)#ordena los datos
##  [1] 177.1429 186.4286 189.2857 191.4286 192.8571 192.8571 195.7143 196.4286
##  [9] 202.1429 205.7143 207.1429 214.2857
#grafico de caja y bigote
boxplot(alimento, horizontal = TRUE)
#(Para usos mas practicos se colocó en horizontal)
#la linea a negra que que esta en el gráfico representa la mediana
#tambien se puede observar que esta mas alargada del lado derecho esto se da dado que hay mas consentración y disperción
#de los datos en este lado de el grafico (no hay numeros extremos)

#medidas de dispersion
#varianza, distancia entre la media y todos los valores
var(alimento)
## [1] 102.1645
#la varianza esta elevada al cuadrado

#la raiz cuadrada de la varianza es la:
#desviasion estandar
sd(alimento)
## [1] 10.10765
#el grado de desviacion de los datos con respecto a su media
#en este caso se ve mas notoria la desviación a partir de la segunda mitad de los datos(ya estando previamente ordenados)

#grafica de dispersion
library(ggplot2)

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

#En esta grafica se puede observar que(en este caso) la cantidad de alimento ingerido por los camarones
#no influye mucho en el peso de los camarones dado que los que mas comen se mantienen con un peso por debajo
# de los 8.5gr