#Estadistica descriptiva de datos de granja acuicola
#Se tiene datos de 12 estanques de 
#alimentacion y peso en semana 12
#12 gramos es el peso ideal
setwd("~/VeranoPyE")
library(readxl)
camarones3 <- read_excel("camarones3.xlsx")
View(camarones3)
#Variable de alimento
alimento <- camarones3$AlimentoDiario
#tabla de frecuencia 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(%)
##  [68.734,72.533) 10 0.83 83.33 10  83.33
##  [72.533,76.331)  0 0.00  0.00 10  83.33
##  [76.331,80.129)  0 0.00  0.00 10  83.33
##  [80.129,83.927)  0 0.00  0.00 10  83.33
##  [83.927,87.726)  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") #histograma de fecuencas acumuladas
plot(dist, type="cfp") #poligono de frecuencias acumuladas
plot(dist, type="fh") #histograma de frecuancia absoluta
plot(dist, type="fh", col="red") #histograma de frecuancia
sort(alimento)
##  [1] 69.42857 69.42857 69.42857 69.42857 69.42857 69.42857 69.42857 69.42857
##  [9] 69.42857 69.42857 86.85714 86.85714
#medidas de tendencia central
mean(alimento) #promedio, media, media aritmetica
## [1] 72.33333
median(alimento)
## [1] 69.42857
summary(alimento)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   69.43   69.43   69.43   72.33   69.43   86.86
#grafico de caja y bigote
boxplot(alimento)
#medidas de dispersion
#varianza, distancia media entre la media y todos los valores
var(alimento)
## [1] 46.0235
#la varianza esta elevada al cuadrado
#la raiz cuadrada de la varianza es la:
#desviacion estandar
sd(alimento)
## [1] 6.784062
#es el grado de desviacion de los datos respecto a su media
#grafica de dispersion
library(ggplot2)

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