#Estadistica descriptiva de datos de granja acuicola
#se tienen datos de 12 tanques
#alimentacion y peso en semana 10
#12 gramos es el peso ideal
setwd("~/PROBABILIDAD Y ESTADISTICA")
library(readxl)
camarones10 <- read_excel("camarones10.xlsx")
View(camarones10)
#variable de alimento
alimento <- camarones10$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")
plot(dist, type = "cfh")
plot(dist, type = "cfp")
plot(dist, type = "fh")
plot(dist, type = "fh", col = "red")
sort(alimento)
## [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
#medidas de tendencia central
mean(alimento) #promedio, media, media aritmetica
## [1] 195.9524
median(alimento) #mediana
## [1] 194.2857
summary(alimento)
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 177.1 190.9 194.3 196.0 203.0 214.3
#grafico de caja y bigote
boxplot(alimento)
#medidas de dispersion
#varianza
var(alimento)
## [1] 102.1645
#desviacion estandar
sd(alimento)
## [1] 10.10765
#grafica de dispersion
#de alimento diario en la semana 12 y el peso de los camarones en gramos
library(ggplot2)
ggplot(data = camarones10) +
geom_point(mapping = aes(x=AlimentoDiario, y=PesoActual))
