## ----setup, include = FALSE--------------------------------------------------- knitr::opts_chunk$set( collapse = TRUE, comment = "#>" ) ## ----------------------------------------------------------------------------- library(gexp) crd <- gexp(mu = 15, err = matrix(rep(0, 6), nrow = 6), r = 3, fe = list(alpha = c(1, -2))) summary(crd) ## ----------------------------------------------------------------------------- crd_mv <- gexp(mu = c(15, 6), err = matrix(c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), ncol = 2), r = 3L, fe = list(alpha = matrix(c(1, -2, 2, 3), ncol = 2))) summary(crd_mv) ## ----------------------------------------------------------------------------- rcbd <- gexp(mu = 15, r = 1, err = matrix(rep(0, 6), nrow = 6), fe = list(alpha = c(1, -2)), blke = c(2, 4, 6), design = 'RCBD') summary(rcbd) ## ----------------------------------------------------------------------------- lsd <- gexp(mu = 15, r = 1, err = matrix(rep(0, 9), nrow = 9), fe = list(alpha = c(1, -2, 3)), rowe = c(2, 3, 4), cole = c(6, 7, 8), design = 'LSD') summary(lsd) ## ----------------------------------------------------------------------------- FE_crd <- gexp(mu = 15, r = 2, err = matrix(rep(0, 12), nrow = 12), fe = list(alpha = c(1, -2), tau = c(1, -1, 1)), inte = c(3, 1, 1, -5, 1, 1), type = 'FE') summary(FE_crd) ## ----------------------------------------------------------------------------- SPE_crd <- gexp(mu = 15, err = matrix(rep(0, 12), nrow = 12), errp = matrix(rep(0, 6), nrow = 6), r = 3, fe = list(alpha = c(1, -2), tau = c(1, -1)), inte = c(3, 1, 1, -5), type = 'SPE') summary(SPE_crd) ## ----------------------------------------------------------------------------- level <- seq(0, 15, 5) cont_crd <- matrix(c(level, level^2, level^3), ncol = 3) crd_l <- gexp(mu = 2, r = 3, err = matrix(rep(0,12), nrow = 12), fe = list(f1 = c(3, 0, 0)), fl = list(dose = level), contrasts = list(dose = cont_crd)) summary(crd_l) ## ----echo=FALSE, results='hide'----------------------------------------------- with(crd_l$dfm, plot(Y1 ~ dose, axes = FALSE, ylim = c(0,50), xlab = 'dose')) axis(1, pos = 0) axis(2, pos = 0) abline(h = 2, col = 'red') abline(2, 3, col = 'red', lty = 2) text(-0.5,2,'2',col = 'red', xpd = TRUE) ## ----------------------------------------------------------------------------- crd_q <- gexp(mu = 2, r = 3, err = matrix(rep(0,12), nrow = 12), fe = list(f1 = c(3, 4, 0)), fl = list(dose = level), contrasts = list(dose = cont_crd)) summary(crd_q) ## ----echo=FALSE, results='hide'----------------------------------------------- with(crd_q$dfm, plot(Y1 ~ dose, axes = FALSE, xlab = 'dose')) axis(1, pos = 0) axis(2, pos = 0) curve(2+3*x+4*x^2, col = 'red', lty = 2, add=T) text(-0.5,2,'2',col = 'red', xpd = TRUE) ## ----------------------------------------------------------------------------- level2 <- seq(2,8,2) cont_crd2 <- matrix(c(level2,level2^2,level2^3), ncol = 3) crd_hb <- gexp(mu = 1, r = 2, err = matrix(rep(0,32), nrow = 32), fe = list(f1 = c(3, 0, 0), f2 = c(2, 0, 0)), fl = list(N = level, P = level2), inte = c(6,0,0,0,0,0,0,0,0), contrasts = list(N = cont_crd, P = cont_crd2), type = 'FE') summary(crd_hb) ## ----echo=FALSE, results='hide'----------------------------------------------- z <- matrix(crd_hb$dfm$Y1, ncol = 4, byrow = FALSE)[seq(1,8,2),] persp(level, level2, z, theta=30, xlab='N', ylab='P', col='red') ## ----------------------------------------------------------------------------- crd_hb2 <- gexp(mu = 1, r = 2, err = matrix(rep(0,32), nrow = 32), fe = list(f1 = c(3, 2, 0), f2 = c(2, 3, 0)), fl = list(N = level, P = level2), inte = c(2,1,0,1,8,0,0,0,0), contrasts = list(N = cont_crd, P = cont_crd2), type = 'FE') summary(crd_hb2) ## ----echo=FALSE, results='hide'----------------------------------------------- z1 <- matrix(crd_hb2$dfm$Y1, ncol=4, byrow=FALSE)[seq(1,8,2),] persp(level, level2, z1, theta=30, xlab='N', ylab='P', col='red') ## ----------------------------------------------------------------------------- crd_design <- gexp() plot(crd_design) ## ----------------------------------------------------------------------------- rcbd_designs <- gexp(r = 1, design = 'RCBD') ## ----fig.width = 7, fig.height = 7-------------------------------------------- plot(rcbd_designs) ## ----------------------------------------------------------------------------- rcbd_designc <- gexp(r = 2, design = 'RCBD') ## ----fig.width = 7, fig.height = 7-------------------------------------------- plot(rcbd_designc) ## ----------------------------------------------------------------------------- lsd_design <- gexp(r = 1, design = 'LSD') ## ----fig.width = 7, fig.height = 7-------------------------------------------- plot(lsd_design) ## ----------------------------------------------------------------------------- fat_design <- gexp(r = 3, type = 'FE') ## ----fig.width = 7, fig.height = 7-------------------------------------------- plot(fat_design) ## ----------------------------------------------------------------------------- split_design <- gexp(r = 3, type = 'SPE') summary(split_design) ## ----fig.width = 7, fig.height = 7-------------------------------------------- plot(split_design) ## ----eval=FALSE--------------------------------------------------------------- # crd_plot <- gexp(r = 3, # fe = list(f1 = rep(1,2))) # plot(crd_plot, # dynamic = TRUE)