Caracoles

Llamar la base de datos

library(readr)
huevos_mortalidad2 <- read_csv("huevos_mortalidad2.csv", 
    col_types = cols(huevos = col_number(), 
        mortalidad = col_number()))
head(huevos_mortalidad2)

## # A tibble: 6 × 6
##   fecha    dia   alimento  luz   huevos mortalidad
##   <chr>    <chr> <chr>     <chr>  <dbl>      <dbl>
## 1 14/09/26 L     Zanahoria con       31         NA
## 2 14/09/26 L     Harina    con      129         NA
## 3 14/09/26 L     Gelatina  con      120         NA
## 4 14/09/26 L     Polvo     con       25         NA
## 5 14/09/26 L     Espinaca  con       38         NA
## 6 14/09/26 L     Zanahoria sin       28         NA

Crear el grafico Crear un grafico de la diferencia entre el lunes y jueves

str(huevos_mortalidad2)

## spc_tbl_ [60 × 6] (S3: spec_tbl_df/tbl_df/tbl/data.frame)
##  $ fecha     : chr [1:60] "14/09/26" "14/09/26" "14/09/26" "14/09/26" ...
##  $ dia       : chr [1:60] "L" "L" "L" "L" ...
##  $ alimento  : chr [1:60] "Zanahoria" "Harina" "Gelatina" "Polvo" ...
##  $ luz       : chr [1:60] "con" "con" "con" "con" ...
##  $ huevos    : num [1:60] 31 129 120 25 38 28 75 85 79 142 ...
##  $ mortalidad: num [1:60] NA NA NA NA NA NA NA NA NA NA ...
##  - attr(*, "spec")=
##   .. cols(
##   ..   fecha = col_character(),
##   ..   dia = col_character(),
##   ..   alimento = col_character(),
##   ..   luz = col_character(),
##   ..   huevos = col_number(),
##   ..   mortalidad = col_number()
##   .. )
##  - attr(*, "problems")=<pointer: 0xba892c2e0>

huevos_mortalidad2$dia <- as.factor(huevos_mortalidad2$dia)
table(huevos_mortalidad2$dia)

## 
##  J  L 
## 30 30

aggregate(huevos ~ dia,
          data = huevos_mortalidad2,
          FUN = mean)

##   dia   huevos
## 1   J 119.9333
## 2   L 223.2333

par(mfrow=c(1,1))

boxplot(huevos ~ dia,
        data = huevos_mortalidad2,
        col = "lightblue2",
        xlab = "dias",
        ylab = "huevos",
        main = "ANOVA por dias")

Gracias por un curso tan sustancioso