######################################### #### Fortgeschrittene Versuchsplanung ### ######## R-Code zu Übungsblatt 2 ######## ######## Unkommentierte Version ######### ######################################### ## Aufgabe 4 ## library(MASS) d1 <- c(0,0,1,1,2,2) d2 <- c(0,0,0,2,2,2) X_d1 <- cbind(rep(1,6),d1,d1^2, deparse.level = 0) X_d2 <- cbind(rep(1,6),d2,d2^2, deparse.level = 0) L1 <- t(c(1,0,0)) L2 <- t(c(0,0,1)) # Lösungen: L1%*%ginv(crossprod(X_d1))%*%t(L1) L2%*%ginv(crossprod(X_d1))%*%t(L2) L1%*%ginv(crossprod(X_d2))%*%t(L1) ## Aufgabe 5 ## library(agricolae) # Optional zur Vergleichbarkeit: set.seed(23) design.crd(c("t1","t2"),c(15,10)) design.crd(c("t1","t2","t3"),rep(10,3))$book