Nothing
tests/testthat/test_calculations.R
In papaja: Prepare American Psychological Association Journal Articles with R Markdown
context("Calculations")
# data <- data.frame(
# "id" = rep(1:2, each = 8)
# , "F1" = rep(c(1, -1), each = 4)
# , "F2" = rep(c(1, -1), each = 2)
# , "F3" = c(1, -1)
# )
#
# data$dv <- data$F1 *rnorm(1) + data$F2 * rnorm(1) + data$F3 * rnorm(1) + rnorm(16) + rep(rnorm(2), each = 8)
#
# for (i in c("F1", "F2", "F3"))
# data[[i]] <- as.factor(data[[i]])
#
# library(Rmisc)
# Rmisc <- summarySEwithin(data = data, measurevar = "dv",
# betweenvars = NULL, withinvars = c("F1", "F2", "F3"), idvar = "id",
# na.rm = FALSE, conf.interval = 0.95, .drop = TRUE)
# Rmisc <- Rmisc[c("F1", "F2", "F3", "ci")]
# colnames(Rmisc)[4] <- "dv"
# dput(as.matrix(Rmisc))
test_that(
"Within-subjects confidence intervals: Three factors within-subjects design"
, {
data <- structure(
list(
id = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L)
, F1 = structure(c(2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L), .Label = c("-1", "1"), class = "factor")
, F2 = structure(c(2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L), .Label = c("-1", "1"), class = "factor")
, F3 = structure(c(2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L), .Label = c("-1", "1"), class = "factor")
, dv = c(
2.51108865147373, 2.79300316288011, 1.24498751760904, -3.2097171749705, -0.86829135419962, -3.17641293936803, -5.09655274729662, -5.67922281667759
, 4.99456366116516, 1.28920687940668, 2.80013312250661, -2.78384472239226, 2.36059467792542, -0.754370127863485, -0.863098014576623, -5.68140019794875
)
)
, .Names = c("id", "F1", "F2", "F3", "dv")
, row.names = c(NA, -16L)
, class = "data.frame"
)
Rmisc <- structure(c("-1", "-1", "-1", "-1", "1", "1", "1", "1", "-1",
"-1", "1", "1", "-1", "-1", "1", "1", "-1", "1", "-1", "1", "-1",
"1", "-1", "1", "10.9180149", "17.8493485", " 5.5466877", "11.0265669",
" 8.0108066", " 0.3410632", "21.1166383", " 5.9639200"), .Dim = c(8L, 4L), .Dimnames = list(NULL, c("F1", "F2", "F3", "dv")))
test <- wsci(data = data, id = "id", dv = "dv", factors = c("F1", "F2", "F3"), method = "Morey")
expect_equal(as.matrix(test), Rmisc)
summary_wsci <- summary(test) # test summary.papaja_wsci
expect_equal(
object = summary_wsci
, expected = structure(
list(
F1 = structure(c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L), .Label = c("-1", "1"), class = "factor")
, F2 = structure(c(1L, 1L, 2L, 2L, 1L, 1L, 2L, 2L), .Label = c("-1", "1"), class = "factor")
, F3 = structure(c(1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L), .Label = c("-1", "1"), class = "factor")
, mean = c(-5.68031150731317, -2.97982538093662, -1.96539153361576, 0.7461516618629, -2.99678094868138, 2.02256032005783, 2.04110502114339, 3.75282615631945)
, lower_limit = c(-16.5983264256805, -20.8291738335329, -7.51207923813347, -10.2804151932351, -11.0075875024338, 1.68149708762336, -19.0755332530318, -2.21109381019798)
, upper_limit = c(5.23770341105414, 14.8695230716597, 3.58129617090195, 11.7727185169609, 5.01402560507109, 2.36362355249229, 23.1577432953186, 9.71674612283687)
)
, row.names = c(NA, -8L)
, class = "data.frame"
)
)
}
)
# library(Rmisc)
# Rmisc <- summarySEwithin(data = data, measurevar = "dv",
# betweenvars = c("F1"), withinvars = c("F2", "F3"), idvar = "id",
# na.rm = FALSE, conf.interval = 0.95, .drop = TRUE)
# Rmisc <- Rmisc[c("F1", "F2", "F3", "ci")]
# dput(Rmisc)
test_that(
"Within-subjects confidence intervals: Three factors mixed design"
, {
data <- structure(
list(
id = c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L)
, F1 = structure(c(2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L), .Label = c("-1", "1"), class = "factor")
, F2 = structure(c(2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L), .Label = c("-1", "1"), class = "factor")
, F3 = structure(c(2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L), .Label = c("-1", "1"), class = "factor")
, dv = c(
2.51108865147373, 2.79300316288011, 1.24498751760904, -3.2097171749705, -0.86829135419962, -3.17641293936803, -5.09655274729662, -5.67922281667759
, 4.99456366116516, 1.28920687940668, 2.80013312250661, -2.78384472239226, 2.36059467792542, -0.754370127863485, -0.863098014576623, -5.68140019794875
)
)
, .Names = c("id", "F1", "F2", "F3", "dv")
, row.names = c(NA, -16L)
, class = "data.frame"
)
Rmisc <- structure(
list(
F1 = structure(c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L), .Label = c("-1", "1"), class = "factor")
, F2 = structure(c(1L, 1L, 2L, 2L, 1L, 1L, 2L, 2L), .Label = c("-1", "1"), class = "factor")
, F3 = structure(c(1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L), .Label = c("-1", "1"), class = "factor")
, ci = c(18.13976813201, 12.9325356320052, 0.355856736124444, 5.56308923612926, 2.30569581586717, 5.97857379893408, 16.4616119853719, 12.788734002305)
)
, .Names = c("F1", "F2", "F3", "ci")
, row.names = c(NA, -8L)
, class = "data.frame"
)
test <- wsci(data = data, id = "id", dv = "dv", factors = c("F1", "F2", "F3"), method = "Morey")
compare_data <- merge(test, Rmisc, by = c("F1", "F2", "F3"))
expect_equal(compare_data$dv, compare_data$ci)
}
)
test_that(
"Within-subjects confidence intervals: Only between-subejcts factors"
, {
data <- npk
data$id <- seq_len(nrow(data))
expect_error(
wsci(data = data, id = "id", dv = "yield", factors = "N", method = "Morey")
, regexp = "No within-subjects factor"
)
}
)
test_that(
"add_effect_sizes: Two-way between-subjects design"
, {
apa_variance_table <- structure(
list(
term = c("(Intercept)", "Gender", "Dosage", "Gender:Dosage")
, sumsq = c(3164.0625, 76.5625, 5.0625, 0.0625)
, df = c(1, 1, 1, 1)
, sumsq_err = c(311.25, 311.25, 311.25, 311.25)
, df.residual = c(12, 12, 12, 12)
, statistic = c(121.987951807229, 2.95180722891566, 0.195180722891566, 0.00240963855421686)
, p.value = c(1.21163629179605e-07, 0.111450651130857, 0.66649556577607, 0.961656681938102)
)
, .Names = c("term", "sumsq", "df", "sumsq_err", "df.residual", "statistic", "p.value")
, row.names = c(NA, -4L)
, class = c("apa_variance_table", "data.frame")
, correction = "none"
)
with_intercept <- papaja:::add_effect_sizes(apa_variance_table, es = c("pes", "ges", "es"), intercept = TRUE)
without_intercept <- papaja:::add_effect_sizes(apa_variance_table, es = c("pes", "ges", "es"), intercept = FALSE)[2:4, ]
# SS from car output:
SS_car <- c(3164.0625, 76.5625, 5.0625, 0.0625, 311.25)
es_car_with <- SS_car[1:4] / sum(SS_car)
es_car_without <- SS_car[2:4] / sum(SS_car[2:5])
tinylabels::variable_label(es_car_with) <-
tinylabels::variable_label(es_car_without) <- "$\\hat{\\eta}^2$"
pes_afex <- c(0.910439708659293, 0.197421434327156, 0.016004742145821, 0.000200762899016262)
tinylabels::variable_label(pes_afex) <- "$\\hat{\\eta}^2_p$"
ges_afex <- pes_afex
tinylabels::variable_label(ges_afex) <- "$\\hat{\\eta}^2_G$"
# Eta-squared
with_intercept <- papaja:::add_effect_sizes(apa_variance_table, es = "es", intercept = TRUE)
expect_identical(
object = with_intercept$estimate
# Expectation calculated via sums of squares from car
, expected = es_car_with
)
without_intercept <- papaja:::add_effect_sizes(apa_variance_table, es = "es", intercept = FALSE)[2:4, ]
expect_equal(
object = without_intercept$estimate
# Expectation calculated via sums of squares from car
, expected = es_car_without
)
# Partial eta-squared
with_intercept <- papaja:::add_effect_sizes(apa_variance_table, es = "pes", intercept = TRUE)
expect_equal(
object = with_intercept$estimate
, expected = pes_afex
)
without_intercept <- papaja:::add_effect_sizes(apa_variance_table, es = "pes", intercept = FALSE)[2:4, ]
expect_equal(
object = without_intercept$estimate
, expected = with_intercept$estimate[2:4]
)
# generalized eta-squared
with_intercept <- papaja:::add_effect_sizes(apa_variance_table, es = "ges", intercept = TRUE)
expect_equal(
object = with_intercept$estimate
, expected = ges_afex
)
without_intercept <- papaja:::add_effect_sizes(apa_variance_table, es = "ges", intercept = FALSE)[2:4, ]
expect_equal(
object = without_intercept$estimate
, expected = with_intercept$estimate[2:4]
)
}
)
test_that(
"add_effect_sizes: Two-way within-subjects design"
, {
apa_variance_table <- structure(
list(
sumsq = c(4177.2, 30.0000000000001, 9.80000000000002, 1.40000000000001)
, df = c(1, 1, 1.43869284635452, 1.62877229248125)
, sumsq_err = c(349.133333333333, 16.3333333333333, 26.8666666666667, 19.2666666666667)
, df.residual = c(4, 4, 5.7547713854181, 6.51508916992501)
, statistic = c(47.8579339316403, 7.34693877551021, 1.4590570719603, 0.290657439446368)
, p.value = c(0.00229109843354983, 0.053508296850918, 0.29341624721021, 0.714981769489886)
, term = c("(Intercept)", "Task", "Valence", "Task:Valence")
)
, .Names = c("sumsq", "df", "sumsq_err", "df.residual", "statistic", "p.value", "term")
, row.names = c(NA, -4L)
, class = c("apa_variance_table", "data.frame")
, correction = "GG"
)
# SS from car output:
SS_car <- c(4177.2, 30.0000000000001, 9.80000000000002, 1.40000000000001)
SS_err <- c(349.133333333333, 16.3333333333333, 26.8666666666667, 19.2666666666667)
es_with_intercept <- SS_car/sum(c(SS_car, SS_err))
es_without_intercept <- SS_car[2:4]/sum(c(SS_car[2:4], SS_err))
tinylabels::variable_label(es_with_intercept) <-
tinylabels::variable_label(es_without_intercept) <- "$\\hat{\\eta}^2$"
tinylabels::variable_label(SS_car) <- "$\\hat{\\eta}^2$"
pes_afex <- c(0.922866190441122, 0.64748201438849, 0.267272727272728, 0.0677419354838713)
ges_afex <- c(0.910303347280335, 0.0679347826086958, 0.0232558139534884, 0.00338983050847459)
tinylabels::variable_label(pes_afex) <- "$\\hat{\\eta}^2_p$"
tinylabels::variable_label(ges_afex) <- "$\\hat{\\eta}^2_G$"
# Eta-squared ----
with_intercept <- papaja:::add_effect_sizes(apa_variance_table, es = "es", intercept = TRUE)
without_intercept <- papaja:::add_effect_sizes(apa_variance_table, es = "es", intercept = FALSE)[2:4, ]
expect_equal(
object = with_intercept$estimate
# Expectation calculated via sums of squares from car
, expected = es_with_intercept
)
expect_equal(
object = without_intercept$estimate
# Expectation calculated via sums of squares from car
, expected = es_without_intercept
)
# Partial eta-squared ----
with_intercept <- add_effect_sizes(apa_variance_table, es = "pes", intercept = TRUE)
without_intercept <- add_effect_sizes(apa_variance_table, es = c("pes"), intercept = FALSE)[2:4, ]
expect_equal(
object = with_intercept$estimate
, expected = pes_afex
)
expect_equal(
object = without_intercept$estimate
, expected = with_intercept$estimate[2:4]
)
# Generalized eta-squared ----
with_intercept <- add_effect_sizes(apa_variance_table, es = "ges", intercept = TRUE)
without_intercept <- add_effect_sizes(apa_variance_table, es = "ges", intercept = TRUE)[2:4, ]
expect_equal(
object = with_intercept$estimate
, expected = ges_afex
)
expect_equal(
object = without_intercept$estimate
, expected = with_intercept$estimate[2:4]
)
}
)
test_that(
"Within-subjects confidence intervals: Cousineau vs. Morey method"
, {
aggregated_data <- stats::aggregate(yield ~ N + block, data = npk, FUN = mean)
object_1 <- wsci(data = aggregated_data, id = "block", dv = "yield", factors = c("N"), method = "Cousineau", level = .98)
object_2 <- wsci(data = aggregated_data, id = "block", dv = "yield", factors = c("N"), method = "Morey", level = .98)
expect_identical(
object = attributes(object_1)[
c("class", "names", "row.names", "Between-subjects factors", "Within-subjects factors"
, "Dependent variable", "Subject identifier", "Confidence level", "Method", "means")
] # default order changed from R 3.4 to 3.5
, expected = list(
class = c("papaja_wsci", "data.frame")
, names = c("N", "yield")
, row.names = c(1L, 2L)
, `Between-subjects factors` = "none"
, `Within-subjects factors` = "N"
, `Dependent variable` = "yield"
, `Subject identifier` = "block"
, `Confidence level` = 0.98
, Method = "Cousineau"
, means = aggregate(yield ~ N, data = aggregated_data, FUN = mean)
)
)
expect_identical(
object = object_1$yield
, expected = object_2$yield /sqrt(2)
)
}
)
test_that(
"Within-subjects confidence intervals: Handling of implicit and explicit NAs"
, {
# Implicit NAs
data <- npk
data$yield[2] <- NA
aggregated <- stats::aggregate(yield ~ N + P + block, data = data, FUN = mean)
expect_warning(
wsci(data = aggregated, id = "block", dv = "yield", factors = c("N", "P"))
, "Because of incomplete data, the following cases were removed from calculation of within-subjects confidence intervals:\nblock: 1"
)
# Explicit NAs
data <- npk
aggregated <- stats::aggregate(yield ~ N + P + block, data = data, FUN = mean)
aggregated$yield[5] <- NA
expect_warning(
wsci(data = aggregated, id = "block", dv = "yield", factors = c("N", "P"))
, "Because of NAs in the dependent variable, the following cases were removed from calculation of within-subjects confidence intervals:\nblock: 2"
)
}
)
test_that(
"Within-subjects confidence intervals: Throw error if data are not properly aggregated."
, {
expect_error(
wsci(data = npk, id = "block", factors = "N", dv = "yield")
, "More than one observation per cell. Ensure you aggregated multiple observations per participant/within-subjects condition combination."
)
}
)
test_that(
"conf_int(): Throw message if not enough observations"
, {
expect_message(conf_int(1), "Less than two non-missing values in at least one cell of your design: Thus, no confidence interval can be computed.")
expect_message(conf_int(NA), "Less than two non-missing values in at least one cell of your design: Thus, no confidence interval can be computed.")
}
)
test_that(
"add_mse(): Warn if information is not sufficient, error if x is not a data frame"
, {
object <- 1:4
expect_error(add_mse(object), "The first argument to add_mse() must inherit from class 'data.frame', but 'object' does not.", fixed = TRUE)
expect_warning(add_mse(data.frame()), "Mean-squared errors requested, but necessary information not available.", fixed = TRUE)
}
)
Try the papaja package in your browser
Any scripts or data that you put into this service are public.
papaja documentation built on Oct. 30, 2024, 9:09 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
context("Calculations")
# data <- data.frame(
# "id" = rep(1:2, each = 8)
# , "F1" = rep(c(1, -1), each = 4)
# , "F2" = rep(c(1, -1), each = 2)
# , "F3" = c(1, -1)
# )
#
# data$dv <- data$F1 *rnorm(1) + data$F2 * rnorm(1) + data$F3 * rnorm(1) + rnorm(16) + rep(rnorm(2), each = 8)
#
# for (i in c("F1", "F2", "F3"))
# data[[i]] <- as.factor(data[[i]])
#
# library(Rmisc)
# Rmisc <- summarySEwithin(data = data, measurevar = "dv",
# betweenvars = NULL, withinvars = c("F1", "F2", "F3"), idvar = "id",
# na.rm = FALSE, conf.interval = 0.95, .drop = TRUE)
# Rmisc <- Rmisc[c("F1", "F2", "F3", "ci")]
# colnames(Rmisc)[4] <- "dv"
# dput(as.matrix(Rmisc))
test_that(
"Within-subjects confidence intervals: Three factors within-subjects design"
, {
data <- structure(
list(
id = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L)
, F1 = structure(c(2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L), .Label = c("-1", "1"), class = "factor")
, F2 = structure(c(2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L), .Label = c("-1", "1"), class = "factor")
, F3 = structure(c(2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L), .Label = c("-1", "1"), class = "factor")
, dv = c(
2.51108865147373, 2.79300316288011, 1.24498751760904, -3.2097171749705, -0.86829135419962, -3.17641293936803, -5.09655274729662, -5.67922281667759
, 4.99456366116516, 1.28920687940668, 2.80013312250661, -2.78384472239226, 2.36059467792542, -0.754370127863485, -0.863098014576623, -5.68140019794875
)
)
, .Names = c("id", "F1", "F2", "F3", "dv")
, row.names = c(NA, -16L)
, class = "data.frame"
)
Rmisc <- structure(c("-1", "-1", "-1", "-1", "1", "1", "1", "1", "-1",
"-1", "1", "1", "-1", "-1", "1", "1", "-1", "1", "-1", "1", "-1",
"1", "-1", "1", "10.9180149", "17.8493485", " 5.5466877", "11.0265669",
" 8.0108066", " 0.3410632", "21.1166383", " 5.9639200"), .Dim = c(8L, 4L), .Dimnames = list(NULL, c("F1", "F2", "F3", "dv")))
test <- wsci(data = data, id = "id", dv = "dv", factors = c("F1", "F2", "F3"), method = "Morey")
expect_equal(as.matrix(test), Rmisc)
summary_wsci <- summary(test) # test summary.papaja_wsci
expect_equal(
object = summary_wsci
, expected = structure(
list(
F1 = structure(c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L), .Label = c("-1", "1"), class = "factor")
, F2 = structure(c(1L, 1L, 2L, 2L, 1L, 1L, 2L, 2L), .Label = c("-1", "1"), class = "factor")
, F3 = structure(c(1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L), .Label = c("-1", "1"), class = "factor")
, mean = c(-5.68031150731317, -2.97982538093662, -1.96539153361576, 0.7461516618629, -2.99678094868138, 2.02256032005783, 2.04110502114339, 3.75282615631945)
, lower_limit = c(-16.5983264256805, -20.8291738335329, -7.51207923813347, -10.2804151932351, -11.0075875024338, 1.68149708762336, -19.0755332530318, -2.21109381019798)
, upper_limit = c(5.23770341105414, 14.8695230716597, 3.58129617090195, 11.7727185169609, 5.01402560507109, 2.36362355249229, 23.1577432953186, 9.71674612283687)
)
, row.names = c(NA, -8L)
, class = "data.frame"
)
)
}
)
# library(Rmisc)
# Rmisc <- summarySEwithin(data = data, measurevar = "dv",
# betweenvars = c("F1"), withinvars = c("F2", "F3"), idvar = "id",
# na.rm = FALSE, conf.interval = 0.95, .drop = TRUE)
# Rmisc <- Rmisc[c("F1", "F2", "F3", "ci")]
# dput(Rmisc)
test_that(
"Within-subjects confidence intervals: Three factors mixed design"
, {
data <- structure(
list(
id = c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L)
, F1 = structure(c(2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L), .Label = c("-1", "1"), class = "factor")
, F2 = structure(c(2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L), .Label = c("-1", "1"), class = "factor")
, F3 = structure(c(2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L), .Label = c("-1", "1"), class = "factor")
, dv = c(
2.51108865147373, 2.79300316288011, 1.24498751760904, -3.2097171749705, -0.86829135419962, -3.17641293936803, -5.09655274729662, -5.67922281667759
, 4.99456366116516, 1.28920687940668, 2.80013312250661, -2.78384472239226, 2.36059467792542, -0.754370127863485, -0.863098014576623, -5.68140019794875
)
)
, .Names = c("id", "F1", "F2", "F3", "dv")
, row.names = c(NA, -16L)
, class = "data.frame"
)
Rmisc <- structure(
list(
F1 = structure(c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L), .Label = c("-1", "1"), class = "factor")
, F2 = structure(c(1L, 1L, 2L, 2L, 1L, 1L, 2L, 2L), .Label = c("-1", "1"), class = "factor")
, F3 = structure(c(1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L), .Label = c("-1", "1"), class = "factor")
, ci = c(18.13976813201, 12.9325356320052, 0.355856736124444, 5.56308923612926, 2.30569581586717, 5.97857379893408, 16.4616119853719, 12.788734002305)
)
, .Names = c("F1", "F2", "F3", "ci")
, row.names = c(NA, -8L)
, class = "data.frame"
)
test <- wsci(data = data, id = "id", dv = "dv", factors = c("F1", "F2", "F3"), method = "Morey")
compare_data <- merge(test, Rmisc, by = c("F1", "F2", "F3"))
expect_equal(compare_data$dv, compare_data$ci)
}
)
test_that(
"Within-subjects confidence intervals: Only between-subejcts factors"
, {
data <- npk
data$id <- seq_len(nrow(data))
expect_error(
wsci(data = data, id = "id", dv = "yield", factors = "N", method = "Morey")
, regexp = "No within-subjects factor"
)
}
)
test_that(
"add_effect_sizes: Two-way between-subjects design"
, {
apa_variance_table <- structure(
list(
term = c("(Intercept)", "Gender", "Dosage", "Gender:Dosage")
, sumsq = c(3164.0625, 76.5625, 5.0625, 0.0625)
, df = c(1, 1, 1, 1)
, sumsq_err = c(311.25, 311.25, 311.25, 311.25)
, df.residual = c(12, 12, 12, 12)
, statistic = c(121.987951807229, 2.95180722891566, 0.195180722891566, 0.00240963855421686)
, p.value = c(1.21163629179605e-07, 0.111450651130857, 0.66649556577607, 0.961656681938102)
)
, .Names = c("term", "sumsq", "df", "sumsq_err", "df.residual", "statistic", "p.value")
, row.names = c(NA, -4L)
, class = c("apa_variance_table", "data.frame")
, correction = "none"
)
with_intercept <- papaja:::add_effect_sizes(apa_variance_table, es = c("pes", "ges", "es"), intercept = TRUE)
without_intercept <- papaja:::add_effect_sizes(apa_variance_table, es = c("pes", "ges", "es"), intercept = FALSE)[2:4, ]
# SS from car output:
SS_car <- c(3164.0625, 76.5625, 5.0625, 0.0625, 311.25)
es_car_with <- SS_car[1:4] / sum(SS_car)
es_car_without <- SS_car[2:4] / sum(SS_car[2:5])
tinylabels::variable_label(es_car_with) <-
tinylabels::variable_label(es_car_without) <- "$\\hat{\\eta}^2$"
pes_afex <- c(0.910439708659293, 0.197421434327156, 0.016004742145821, 0.000200762899016262)
tinylabels::variable_label(pes_afex) <- "$\\hat{\\eta}^2_p$"
ges_afex <- pes_afex
tinylabels::variable_label(ges_afex) <- "$\\hat{\\eta}^2_G$"
# Eta-squared
with_intercept <- papaja:::add_effect_sizes(apa_variance_table, es = "es", intercept = TRUE)
expect_identical(
object = with_intercept$estimate
# Expectation calculated via sums of squares from car
, expected = es_car_with
)
without_intercept <- papaja:::add_effect_sizes(apa_variance_table, es = "es", intercept = FALSE)[2:4, ]
expect_equal(
object = without_intercept$estimate
# Expectation calculated via sums of squares from car
, expected = es_car_without
)
# Partial eta-squared
with_intercept <- papaja:::add_effect_sizes(apa_variance_table, es = "pes", intercept = TRUE)
expect_equal(
object = with_intercept$estimate
, expected = pes_afex
)
without_intercept <- papaja:::add_effect_sizes(apa_variance_table, es = "pes", intercept = FALSE)[2:4, ]
expect_equal(
object = without_intercept$estimate
, expected = with_intercept$estimate[2:4]
)
# generalized eta-squared
with_intercept <- papaja:::add_effect_sizes(apa_variance_table, es = "ges", intercept = TRUE)
expect_equal(
object = with_intercept$estimate
, expected = ges_afex
)
without_intercept <- papaja:::add_effect_sizes(apa_variance_table, es = "ges", intercept = FALSE)[2:4, ]
expect_equal(
object = without_intercept$estimate
, expected = with_intercept$estimate[2:4]
)
}
)
test_that(
"add_effect_sizes: Two-way within-subjects design"
, {
apa_variance_table <- structure(
list(
sumsq = c(4177.2, 30.0000000000001, 9.80000000000002, 1.40000000000001)
, df = c(1, 1, 1.43869284635452, 1.62877229248125)
, sumsq_err = c(349.133333333333, 16.3333333333333, 26.8666666666667, 19.2666666666667)
, df.residual = c(4, 4, 5.7547713854181, 6.51508916992501)
, statistic = c(47.8579339316403, 7.34693877551021, 1.4590570719603, 0.290657439446368)
, p.value = c(0.00229109843354983, 0.053508296850918, 0.29341624721021, 0.714981769489886)
, term = c("(Intercept)", "Task", "Valence", "Task:Valence")
)
, .Names = c("sumsq", "df", "sumsq_err", "df.residual", "statistic", "p.value", "term")
, row.names = c(NA, -4L)
, class = c("apa_variance_table", "data.frame")
, correction = "GG"
)
# SS from car output:
SS_car <- c(4177.2, 30.0000000000001, 9.80000000000002, 1.40000000000001)
SS_err <- c(349.133333333333, 16.3333333333333, 26.8666666666667, 19.2666666666667)
es_with_intercept <- SS_car/sum(c(SS_car, SS_err))
es_without_intercept <- SS_car[2:4]/sum(c(SS_car[2:4], SS_err))
tinylabels::variable_label(es_with_intercept) <-
tinylabels::variable_label(es_without_intercept) <- "$\\hat{\\eta}^2$"
tinylabels::variable_label(SS_car) <- "$\\hat{\\eta}^2$"
pes_afex <- c(0.922866190441122, 0.64748201438849, 0.267272727272728, 0.0677419354838713)
ges_afex <- c(0.910303347280335, 0.0679347826086958, 0.0232558139534884, 0.00338983050847459)
tinylabels::variable_label(pes_afex) <- "$\\hat{\\eta}^2_p$"
tinylabels::variable_label(ges_afex) <- "$\\hat{\\eta}^2_G$"
# Eta-squared ----
with_intercept <- papaja:::add_effect_sizes(apa_variance_table, es = "es", intercept = TRUE)
without_intercept <- papaja:::add_effect_sizes(apa_variance_table, es = "es", intercept = FALSE)[2:4, ]
expect_equal(
object = with_intercept$estimate
# Expectation calculated via sums of squares from car
, expected = es_with_intercept
)
expect_equal(
object = without_intercept$estimate
# Expectation calculated via sums of squares from car
, expected = es_without_intercept
)
# Partial eta-squared ----
with_intercept <- add_effect_sizes(apa_variance_table, es = "pes", intercept = TRUE)
without_intercept <- add_effect_sizes(apa_variance_table, es = c("pes"), intercept = FALSE)[2:4, ]
expect_equal(
object = with_intercept$estimate
, expected = pes_afex
)
expect_equal(
object = without_intercept$estimate
, expected = with_intercept$estimate[2:4]
)
# Generalized eta-squared ----
with_intercept <- add_effect_sizes(apa_variance_table, es = "ges", intercept = TRUE)
without_intercept <- add_effect_sizes(apa_variance_table, es = "ges", intercept = TRUE)[2:4, ]
expect_equal(
object = with_intercept$estimate
, expected = ges_afex
)
expect_equal(
object = without_intercept$estimate
, expected = with_intercept$estimate[2:4]
)
}
)
test_that(
"Within-subjects confidence intervals: Cousineau vs. Morey method"
, {
aggregated_data <- stats::aggregate(yield ~ N + block, data = npk, FUN = mean)
object_1 <- wsci(data = aggregated_data, id = "block", dv = "yield", factors = c("N"), method = "Cousineau", level = .98)
object_2 <- wsci(data = aggregated_data, id = "block", dv = "yield", factors = c("N"), method = "Morey", level = .98)
expect_identical(
object = attributes(object_1)[
c("class", "names", "row.names", "Between-subjects factors", "Within-subjects factors"
, "Dependent variable", "Subject identifier", "Confidence level", "Method", "means")
] # default order changed from R 3.4 to 3.5
, expected = list(
class = c("papaja_wsci", "data.frame")
, names = c("N", "yield")
, row.names = c(1L, 2L)
, `Between-subjects factors` = "none"
, `Within-subjects factors` = "N"
, `Dependent variable` = "yield"
, `Subject identifier` = "block"
, `Confidence level` = 0.98
, Method = "Cousineau"
, means = aggregate(yield ~ N, data = aggregated_data, FUN = mean)
)
)
expect_identical(
object = object_1$yield
, expected = object_2$yield /sqrt(2)
)
}
)
test_that(
"Within-subjects confidence intervals: Handling of implicit and explicit NAs"
, {
# Implicit NAs
data <- npk
data$yield[2] <- NA
aggregated <- stats::aggregate(yield ~ N + P + block, data = data, FUN = mean)
expect_warning(
wsci(data = aggregated, id = "block", dv = "yield", factors = c("N", "P"))
, "Because of incomplete data, the following cases were removed from calculation of within-subjects confidence intervals:\nblock: 1"
)
# Explicit NAs
data <- npk
aggregated <- stats::aggregate(yield ~ N + P + block, data = data, FUN = mean)
aggregated$yield[5] <- NA
expect_warning(
wsci(data = aggregated, id = "block", dv = "yield", factors = c("N", "P"))
, "Because of NAs in the dependent variable, the following cases were removed from calculation of within-subjects confidence intervals:\nblock: 2"
)
}
)
test_that(
"Within-subjects confidence intervals: Throw error if data are not properly aggregated."
, {
expect_error(
wsci(data = npk, id = "block", factors = "N", dv = "yield")
, "More than one observation per cell. Ensure you aggregated multiple observations per participant/within-subjects condition combination."
)
}
)
test_that(
"conf_int(): Throw message if not enough observations"
, {
expect_message(conf_int(1), "Less than two non-missing values in at least one cell of your design: Thus, no confidence interval can be computed.")
expect_message(conf_int(NA), "Less than two non-missing values in at least one cell of your design: Thus, no confidence interval can be computed.")
}
)
test_that(
"add_mse(): Warn if information is not sufficient, error if x is not a data frame"
, {
object <- 1:4
expect_error(add_mse(object), "The first argument to add_mse() must inherit from class 'data.frame', but 'object' does not.", fixed = TRUE)
expect_warning(add_mse(data.frame()), "Mean-squared errors requested, but necessary information not available.", fixed = TRUE)
}
)
Try the papaja package in your browser
Any scripts or data that you put into this service are public.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.