#Estadistica descriptiva de datos de granja acuicola
#Se tienen datos de 12 estanques de
#alimentacion y peso en semana 6
#12 gramos es el peso normal
setwd("~/VERANOPYE")
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
#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")
#f= frecuencia absoluta
#rf= frecuencia relativa
#rf(%)= frecuencia relativa porcentual
#cf= frecuencia acumulada
#cf(%)= frecuencia acumulada porcentual
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, media, media aritmetica
## [1] 128.4524
median(alimento) #mediana
## [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
#grafico de caja y bigote
boxplot(alimento)
#medidas de dispersion
#varianza
var(alimento)
## [1] 9.353742
#la varianza esta elevada al cuadrado
#la raiz cuadrada de la varianza es la:
#desviacion estandar
sd(alimento)
## [1] 3.058389
#el grado de desviacion de los datos con respecto a su media
#grafica de dispersion
#del alimento diario en la semana 12 en gramos y el peso de los
#camarones en gramos
library(ggplot2)
ggplot(data = camarones6) +
geom_point(mapping = aes(x=AlimentoDiario, y=PesoActual))
