Systematic sampling - variance estimation

• Tasso di umidità in cloni di Atriplex in vivaio sitemati su una griglia 3 x 11 ogni tripletta una stessa pianta madre

Accesso ai dati

library(googlesheets)
suppressMessages(library(dplyr))
gs_ls("DATI")

## Auto-refreshing stale OAuth token.

## # A tibble: 3 x 10
##                     sheet_title          author  perm version
##                           <chr>           <chr> <chr>   <chr>
## 1                          DATI  chiaramura1984    rw     old
## 2 DATI_ALTEZZE_CANALI_ESCA.xlsx          scotti    rw     old
## 3               DATI PEDOLOGICI antoniodesulo92    rw     old
## # ... with 6 more variables: updated <dttm>, sheet_key <chr>,
## #   ws_feed <chr>, alternate <chr>, self <chr>, alt_key <chr>

dati <- gs_key("1XeIA3mRQleBMCI9YwcocMFagIvPxS_mSfgZBj8p0Txc")

## Sheet successfully identified: "DATI"

gs_ws_ls(dati)

## [1] "DATI"                      "TabBase"                  
## [3] "Tabella pivot 1"           "Campionamento casuale"    
## [5] "Passo di campionamento"    "Campionamento sistematico"

tabBase <- gs_read(dati, ws = "TabBase")

## Accessing worksheet titled 'TabBase'.

## Parsed with column specification:
## cols(
##   prog = col_integer(),
##   colonna = col_character(),
##   riga = col_integer(),
##   loc = col_character(),
##   clone = col_character(),
##   `peso fresco` = col_number(),
##   `tasso di umidità` = col_character()
## )

Sistema i dati

tabBase <- tabBase %>%
  transmute(
    prog = as.factor(prog),
    colonna = as.factor(colonna),
    riga = as.factor(riga),
    loc = as.factor(loc) ,
    clone = as.factor(clone),
    umidita = as.numeric(gsub(",", ".", tabBase$`tasso di umidità`)))

Grafici

library(ggplot2)
ggplot(tabBase, aes(loc, umidita)) +
  geom_boxplot()

ggplot(tabBase, aes(riga, umidita)) +
  geom_boxplot() +
  geom_point()