#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("~/Stat")
library(readxl)
camarones11 <- read_excel("camarones11.xlsx")
View(camarones11)

#Variable de alimento
alimento <- camarones11$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(%)
##  [166.179,189.936) 2 0.17 16.67  2  16.67
##  [189.936,213.693) 2 0.17 16.67  4  33.33
##   [213.693,237.45) 3 0.25 25.00  7  58.33
##   [237.45,261.207) 1 0.08  8.33  8  66.67
##  [261.207,284.964) 4 0.33 33.33 12 100.00
#partir la pantalla en 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

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


sort(alimento)
##  [1] 167.8571 175.0000 192.8571 207.1429 228.5714 228.5714 235.7143 260.7143
##  [9] 275.0000 278.5714 282.1429 282.1429
#medidas de tendencia central
mean(alimento)  #promedio, media, media aritmética
## [1] 234.5238
median(alimento) #mediana 
## [1] 232.1429
summary(alimento)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   167.9   203.6   232.1   234.5   275.9   282.1
#grafico de caja y bigote
boxplot(alimento)

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

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

#grafica de dispersion
#del alimento diario en la semana 12 en gramos y el peso de los 
#camarones en gramos
library(ggplot2)

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

#Puedes encontrarlo en:
#https://rpubs.com/Manuel_Lopez/624451