Nothing
context("gls models")
# ------------------------------------------------------------------------------
library(nlme)
tmp_ovary <- Ovary
tmp_ovary$x <- factor(sample(0:1, size = nrow(Ovary), replace = TRUE))
gls_fit_1 <- gls(logDens ~ Block + dilut,
data = Assay,
correlation = corCompSymm(form = ~ 1 | sample)
)
test_that("gls models", {
a_1 <- list(Block = "1", dilut = "1")
b_1 <- list(Block = "2", dilut = "1")
c_1 <- contrast(gls_fit_1, a_1, b_1)
expect_equivalent(
summary(gls_fit_1)$tTable["Block.L", "p-value"],
c_1$Pvalue
)
})
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.