Nothing
#devtools::test("dae")
context("analysis")
cat("#### Test for Tukey.1df using Kirk example\n")
test_that("Tukey1dfKirk", {
skip_on_cran()
library(dae)
#'## Dataset from Kirk (Table 7.2-1)
kirk721 <- data.frame(su = factor(rep(1:8,each=4)),
A = factor(rep(1:4,8)),
y = c( 3,4,4,3, 2,4,4,5, 2,3,3,6, 3,3,3,5,
1,2,4,7, 3,3,6,6, 4,4,5,10, 6,5,5,8)
)
#'### Add a unit factor
kirk721$unit <- kirk721$A
#'## Test tukey.1df - su + A
#'### No error term
fit <- aov(y ~ su + A, kirk721)
kirk1df <- tukey.1df(fit,kirk721,"Residuals")
testthat::expect_lt(abs(kirk1df$Tukey.SS - 1.773891), 1e-5)
testthat::expect_lt(abs(kirk1df$Tukey.F - 1.279581), 1e-5)
testthat::expect_lt(abs(kirk1df$Tukey.p - 0.271392), 1e-4)
testthat::expect_lt(abs(kirk1df$Devn.SS - 27.72611), 1e-5)
fit <- aov(y ~ su + A + Error(su), kirk721)
kirk1df <- tukey.1df(fit,kirk721)
testthat::expect_lt(abs(kirk1df$Tukey.F - 1.279581), 1e-5)
fit <- aov(y ~ su + A + Error(su/unit), kirk721)
kirk1df <- tukey.1df(fit,kirk721, "su:unit")
testthat::expect_lt(abs(kirk1df$Tukey.F - 1.279581), 1e-5)
#'### su not outside Error
fit <- aov(y ~ A + Error(su), kirk721)
kirk1df <- tukey.1df(fit,kirk721)
testthat::expect_gt(abs(kirk1df$Tukey.F - 1.279581), 1e-01)
fit <- aov(y ~ A + Error(su/unit), kirk721)
kirk1df <- tukey.1df(fit,kirk721, "su:unit")
testthat::expect_gt(abs(kirk1df$Tukey.F - 1.279581), 1e-01)
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.