BIBD - funcția renală la pisici

Pentru a înțelege cu adevărat cum arată un BIBD, să construim de la zero un set de date care să conțină un astfel de design.

Să presupunem că vrem să testăm efectul a patru tipuri diferite de hrană umedă pentru pisici asupra funcției renale. Hrana pentru pisici este însă costisitoare, așa că vom testa doar 3 tipuri per bloc pentru a economisi. Factorul de blocare este color (culoarea pisicii), deoarece acesta nu face parte din experimentul nostru. Rezultatul măsurat va fi nivelul de creatinină din sânge, un indicator al funcției renale și al disfuncțiilor renale atât la pisici, cât și la oameni.

Acest exercițiu face parte din cursul

Design Experimental în R

Vezi cursul

Instrucțiuni pentru exercițiu

Funcția personalizată lambda() a fost deja încărcată. Calculează lambda cu t = 4, k = 3 și r = 3 pentru a verifica dacă un BIBD este posibil.
Rulează codul pentru a asamble setul de date. Poți vedea ordinea în care sunt utilizate tratamentele cu hrană în fiecare bloc.
Creează cat_model cu aov() conform descrierii experimentului de mai sus și examinează rezultatele cu summary(). Tipul de hrană umedă influențează nivelurile de creatinină?

Exercițiu interactiv practic

Încearcă acest exercițiu completând acest cod de exemplu.

# 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 = ___)
___

Editează și rulează codul