U1A4.R

#estadistica descriptiva de granja acuicola
#se tienen datos de 12 estanques de
#alimentacion y peso en semana 6
#12 gramos es el peso ideal
setwd("~/verano")
library(readxl)
camarones6 <- read_excel("camarones6.xlsx")
View(camarones6)

#variable de alimento
alimento <- camarones6$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(%)
##  [125.871,127.967) 10 0.83 83.33 10  83.33
##  [127.967,130.063)  0 0.00  0.00 10  83.33
##  [130.063,132.159)  0 0.00  0.00 10  83.33
##  [132.159,134.254)  0 0.00  0.00 10  83.33
##   [134.254,136.35)  2 0.17 16.67 12 100.00

par(mfrow=c(3,2))
hist(alimento, breaks = "Sturges", xlab = "Temperatura en °C", ylab = "Frecuencia", main = paste("Histograma de Guasave"), border = (color = "blue") )
plot(dist, type = "cfh")
plot(dist, type = "cfp")
plot(dist, type = "fh")
plot(dist, type = "fh",col = "gold")


sort(alimento)

##  [1] 127.1429 127.1429 127.1429 127.1429 127.1429 127.1429 127.1429 127.1429
##  [9] 127.1429 127.1429 135.0000 135.0000

#medidas de tendencia central
mean(alimento) #promedio

## [1] 128.4524

median(alimento)

## [1] 127.1429

summary(alimento)

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   127.1   127.1   127.1   128.5   127.1   135.0

boxplot(alimento)

var(alimento)

## [1] 9.353742

#desviacion estandar
sd(alimento)

## [1] 3.058389

#grafica de dispercion
#del alimento diario en la semana 6 en gramos
#y el peso de los camarones en gramos

library(ggplot2)

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