U1A4.R

#Estadistica descriptiva de datos de granja acuicola
#Se tienen datos de 12 estanques de 
#Alimentacion y peso en la semana 3
#12 gramos es el peso ideal
setwd("~/PyE")
library(readxl)
camarones3 <- read_excel("camarones3.xlsx")
View(camarones3)

#Variable
alimento<- camarones3$AlimentoDiario

#tabla de frecuecnias de alimento
library(fdth)

## 
## Attaching package: 'fdth'

## The following objects are masked from 'package:stats':
## 
##     sd, var

dist <- fdt(alimento, braks="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, braks="Sturges")

## Warning in plot.window(xlim, ylim, "", ...): "braks" is not a graphical
## parameter

## Warning in title(main = main, sub = sub, xlab = xlab, ylab = ylab, ...): "braks"
## is not a graphical parameter

## Warning in axis(1, ...): "braks" is not a graphical parameter

## Warning in axis(2, ...): "braks" is not a graphical parameter

plot(dist, type="cfh") 
plot(dist, type="cfp") 
plot(dist, type="fh") 
plot(dist, type="fh", col="green") 

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

## [1] 72.33333

median(alimento)#Mediana

## [1] 69.42857

summary(alimento)#Resumen

##    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
var(alimento)

## [1] 46.0235

#desviacion estandar
sd(alimento)

## [1] 6.784062

#Grafico de dispersion 
#del alimento diario en la semana 3 
library(ggplot2) 

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