ejercicio-3.R

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

#Variable de alimento
alimento <- camarones1$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(%)
##   [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 en un arreglo de 3x2
par(mfrow=c(3,2))
hist(alimento,breaks = "Sturges")
plot(dist, tyoe="cfh")

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

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

## Warning in rect(brk[-length(brk)], 0, brk[-1], y, col = col, ...): "tyoe" is not
## a graphical parameter

## Warning in axis(1, at = round(brk, x.round), las = xlas, ...): "tyoe" is not a
## graphical parameter

plot(dist, tyoe="cfp")

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

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

## Warning in rect(brk[-length(brk)], 0, brk[-1], y, col = col, ...): "tyoe" is not
## a graphical parameter

## Warning in axis(1, at = round(brk, x.round), las = xlas, ...): "tyoe" is not a
## graphical parameter

plot(dist, tyoe="fh")

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

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

## Warning in rect(brk[-length(brk)], 0, brk[-1], y, col = col, ...): "tyoe" is not
## a graphical parameter

## Warning in axis(1, at = round(brk, x.round), las = xlas, ...): "tyoe" is not a
## graphical parameter

plot(dist, tyoe="fh",col="green")

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

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

## Warning in rect(brk[-length(brk)], 0, brk[-1], y, col = col, ...): "tyoe" is not
## a graphical parameter

## Warning in axis(1, at = round(brk, x.round), las = xlas, ...): "tyoe" is not a
## graphical parameter

#f=frecuencia absoluta
#rf= Frecuencia relativa
#rf(%)= frecuencia relativa porcentual
#cf= frecuencia acumulada
#cf(%)= frecuencia acumulada porcentual

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 aritmetica

## [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 dispersion 
#varianza, distancia media entre la media y todos los valores
var(alimento)

## [1] 90.05875

#la varianza esta elevada al cuadrado 

#la raiz cuadrada de la varianza es la:
#desviacion estandar
sd(alimento)

## [1] 9.489929

#el grado de desviacion de los datos con respecto a su media

#grafica de dispersion
library(ggplot2)

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