BIBD – funkcja nerek kotów
Żeby naprawdę dobrze zrozumieć, jak wygląda BIBD, zbudujemy zbiór danych zawierający taki układ od podstaw.
Przypuśćmy, że chcemy zbadać wpływ czterech różnych karm mokrych na funkcję nerek kotów. Karma mokra jest jednak droga, więc w każdym bloku przetestujemy tylko 3 rodzaje karmy, aby ograniczyć koszty. Czynnikiem blokującym jest color (kolor sierści kota), który nie jest przedmiotem naszego eksperymentu. Zmienną wynikową będzie zmierzony poziom kreatyniny we krwi – wskaźnik funkcji i dysfunkcji nerek zarówno u kotów, jak i u ludzi.
To ćwiczenie jest częścią kursu
Projektowanie eksperymentów w R
Instrukcje do ćwiczenia
- Niestandardowa funkcja
lambda()została już wczytana. Oblicz wartośćlambdadlat= 4,k= 3 ir= 3, aby sprawdzić, czy BIBD jest możliwy. - Uruchom kod, aby złożyć zbiór danych. Możesz zobaczyć kolejność, w jakiej zabiegi żywieniowe są stosowane w każdym bloku.
- Utwórz model
cat_modelza pomocą funkcjiaov()zgodnie z opisem eksperymentu powyżej i zbadaj wyniki za pomocą funkcjisummary(). Czy rodzaj karmy mokrej ma wpływ na poziom kreatyniny?
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
# 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 = ___)
___