#Estadística descriptiva de datos de granja acuicola
#Se tienen datos de 12 estanque
#de alimentos y peso en semana 12
#12 gramos es el peso ideal
setwd("~/VERANOPYE")
library(readxl)
camarones1 <- read_excel("camarones1.xlsx")
View(camarones1)
#Variable de alimento
alimento <- camarones1$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(%)
## [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 ne 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
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
#Grafico de caja y bigote
boxplot(alimento)
#Medidas de dispersión
#Varianza, distancia entre media entre la media y todos los valores
var(alimento)
## [1] 90.05875
#La varianza está elevada al cuadrado
#La raiz cuadrada de la varianza es la:
#Desviación estándar
sd(alimento)
## [1] 9.489929
#El grado de desviación de los datos con respecto a su media
#Gráfica de dispersión
#del alimento ndiario de 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))
