BIBD - funzionalità renale del gatto
Per essere sicuri di aver capito davvero com’è fatto un BIBD, costruiamo da zero un insieme di dati che lo rispetti.
Supponiamo di voler testare la differenza tra quattro diversi cibi umidi nelle diete dei gatti sulla loro funzionalità renale. Tuttavia, il cibo per gatti è costoso, quindi testeremo solo 3 cibi per blocco per risparmiare. Il fattore di blocco è il color del gatto, perché non ci interessa come parte del nostro esperimento. L’esito sarà il livello di creatinina nel sangue, un indicatore della funzionalità (e disfunzione) renale sia nei gatti sia negli esseri umani.
Questo esercizio fa parte del corso
Progettazione Sperimentale in R
Istruzioni dell'esercizio
- La funzione personalizzata
lambda()è già stata caricata per te. Calcolalambdacont= 4,k= 3 er= 3 per verificare che un BIBD sia possibile. - Esegui il codice per assemblare l’insieme di dati. Puoi vedere l’ordine in cui i trattamenti (i cibi) vengono usati in ciascun blocco.
- Crea
cat_modelconaov()secondo la descrizione dell’esperimento sopra ed esamina i risultati consummary(). Il tipo di cibo umido fa davvero la differenza sui livelli di creatinina?
Esercizio pratico interattivo
Prova a risolvere questo esercizio completando il codice di esempio.
# Calculate lambda
___
# Build the data.frame
creatinine <- c(1.98, 1.97, 2.35, 2.09, 1.87, 1.95, 2.08, 2.01, 1.84, 2.06, 1.97, 2.22)
food <- as.factor(c("A", "C", "D", "A", "B", "C", "B", "C", "D", "A", "B", "D"))
color <- as.factor(rep(c("Black", "White", "Orange", "Spotted"), each = 3))
cat_experiment <- as.data.frame(cbind(creatinine, food, color))
# Create cat_model and examine with summary()
___ <- ___ (___ ~ ___ + color, data = ___)
___