#Estadistica desccriptiva de datos de granja acuicola
#Se tienen datos de 12 estanques 
#alimentacion Y peso en semana 8
#12 gramos es el peso ideal
setwd("~/ProbYEstadistica")
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 en un arreglo de 3x2
par(mfrow=c(3,2))
hist(alimento,breaks="Sturges")
plot(dist, type=("cfh"))
plot(dist, type=("cfp"))
plot(dist, type=("fh"))
plot(dist, type ="fh", col="orange")

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
mean(alimento) #promedio, media, media aritmetica
## [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 dispersion 
#varianza, distanica media ente la media y todos los valores 
var(alimento)
## [1] 31.67517
#suma de la varianza y la media (90)
#la raiz cuadrada de la varianza
#deviacion estandar 
sd(alimento )
## [1] 5.62807
#el grado  de desviacion de los datos con respecto 
#grafica de dispersion
#del alimento diario en la semana 8 en gramos y el peso 
#de los camarones en gramos 
library(ggplot2)

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

#Este documento en linea en el URL