Nothing
tests/testthat/test.R
In multiplestressR: Additive and Multiplicative Null Models for Multiple Stressor Data
test_that("effect_size_additive works as expected", {
expect_equal(effect_size_additive(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.05)$Interaction_Effect_Size[1], -0.4949691273)
expect_equal(effect_size_additive(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.10)$Interaction_CI_Lower[1], -2.0407155118)
expect_equal(effect_size_additive(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.75)$Null_Model[1], "Additive")
expect_equal(effect_size_additive(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.075,
Small_Sample_Correction = FALSE)$Interaction_Effect_Size[1], -0.5287170223)
expect_equal(effect_size_additive(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.01,
Small_Sample_Correction = FALSE)$Interaction_CI_Upper[1], 2.0583386246)
expect_equal(effect_size_additive(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.50,
Small_Sample_Correction = FALSE)$Null_Model[1], "Additive")
})
test_that("effect_size_multiplicative works as expected", {
expect_equal(effect_size_multiplicative(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.05)$Interaction_Effect_Size[1], -0.0705007812)
expect_equal(effect_size_multiplicative(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.05)$Interaction_CI_Lower[1], -0.2353032640)
expect_equal(effect_size_multiplicative(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.75)$Null_Model[1], "Multiplicative")
})
test_that("critical_effect_size_additive works as expected", {
expect_equal(critical_effect_size_additive(Control_N = 4,
StressorA_N = 4,
StressorB_N = 4,
StressorsAB_N = 4), 1.9397344512)
expect_equal(critical_effect_size_additive(Control_N = 3,
StressorA_N = 3,
StressorB_N = 3,
StressorsAB_N = 3,
Small_Sample_Correction = FALSE,
Significance_Level = 0.10), 2.0163641416)
expect_equal(critical_effect_size_additive(Control_N = 46,
StressorA_N = 23,
StressorB_N = 73,
StressorsAB_N = 4,
Small_Sample_Correction = TRUE,
Significance_Level = 0.50), 0.3850781858)
})
test_that("classify_interactions works as expected", {
expect_equal(classify_interactions(effect_size_additive(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.05),
assign_reversals = TRUE,
remove_directionality = FALSE)$Interaction_Classification[1], "Null")
expect_equal(classify_interactions(effect_size_additive(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.5),
assign_reversals = FALSE,
remove_directionality = TRUE)$Interaction_Classification[1], "Null")
expect_equal(classify_interactions(effect_size_multiplicative(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.01),
assign_reversals = TRUE,
remove_directionality = TRUE)$Interaction_Classification[1], "Null")
expect_equal(classify_interactions(effect_size_multiplicative(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.005),
assign_reversals = FALSE,
remove_directionality = FALSE)$Interaction_Classification[1], "Null")
})
test_that("summary_plots works as expected", {
expect_equal(summary_plots(classify_interactions(effect_size_additive(Control_N = multiplestressR::survival$Sample_Size_Control,
Control_SD = multiplestressR::survival$Standard_Deviation_Control,
Control_Mean = multiplestressR::survival$Mean_Control,
StressorA_N = multiplestressR::survival$Sample_Size_Temperature,
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature,
StressorA_Mean = multiplestressR::survival$Mean_Temperature,
StressorB_N = multiplestressR::survival$Sample_Size_pH,
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH,
StressorB_Mean = multiplestressR::survival$Mean_pH,
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH,
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH,
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH,
Significance_Level = 0.05),
assign_reversals = TRUE,
remove_directionality = FALSE),
Significance_Level = 0.05)[[1]]$data[1,3], 0.764)
expect_equal(summary_plots(classify_interactions(effect_size_additive(Control_N = multiplestressR::survival$Sample_Size_Control,
Control_SD = multiplestressR::survival$Standard_Deviation_Control,
Control_Mean = multiplestressR::survival$Mean_Control,
StressorA_N = multiplestressR::survival$Sample_Size_Temperature,
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature,
StressorA_Mean = multiplestressR::survival$Mean_Temperature,
StressorB_N = multiplestressR::survival$Sample_Size_pH,
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH,
StressorB_Mean = multiplestressR::survival$Mean_pH,
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH,
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH,
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH,
Significance_Level = 0.05),
assign_reversals = TRUE,
remove_directionality = FALSE),
Significance_Level = 0.05)[[2]]$data$Sample_Size_Median[123], 6)
expect_equal(summary_plots(classify_interactions(effect_size_additive(Control_N = multiplestressR::survival$Sample_Size_Control,
Control_SD = multiplestressR::survival$Standard_Deviation_Control,
Control_Mean = multiplestressR::survival$Mean_Control,
StressorA_N = multiplestressR::survival$Sample_Size_Temperature,
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature,
StressorA_Mean = multiplestressR::survival$Mean_Temperature,
StressorB_N = multiplestressR::survival$Sample_Size_pH,
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH,
StressorB_Mean = multiplestressR::survival$Mean_pH,
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH,
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH,
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH,
Significance_Level = 0.05),
assign_reversals = TRUE,
remove_directionality = FALSE),
Significance_Level = 0.05)[[3]]$labels$y, "Density")
expect_equal(summary_plots(classify_interactions(effect_size_additive(Control_N = multiplestressR::survival$Sample_Size_Control,
Control_SD = multiplestressR::survival$Standard_Deviation_Control,
Control_Mean = multiplestressR::survival$Mean_Control,
StressorA_N = multiplestressR::survival$Sample_Size_Temperature,
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature,
StressorA_Mean = multiplestressR::survival$Mean_Temperature,
StressorB_N = multiplestressR::survival$Sample_Size_pH,
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH,
StressorB_Mean = multiplestressR::survival$Mean_pH,
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH,
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH,
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH,
Significance_Level = 0.05),
assign_reversals = TRUE,
remove_directionality = FALSE),
Significance_Level = 0.05)[[4]]$labels$y, "Inverse Variance")
expect_equal(summary_plots(classify_interactions(effect_size_multiplicative(Control_N = multiplestressR::survival$Sample_Size_Control,
Control_SD = multiplestressR::survival$Standard_Deviation_Control,
Control_Mean = multiplestressR::survival$Mean_Control,
StressorA_N = multiplestressR::survival$Sample_Size_Temperature,
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature,
StressorA_Mean = multiplestressR::survival$Mean_Temperature,
StressorB_N = multiplestressR::survival$Sample_Size_pH,
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH,
StressorB_Mean = multiplestressR::survival$Mean_pH,
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH,
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH,
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH,
Significance_Level = 0.05),
assign_reversals = TRUE,
remove_directionality = FALSE),
Significance_Level = 0.05)[[5]]$labels$y, "Standard Error")
})
Try the multiplestressR package in your browser
Any scripts or data that you put into this service are public.
multiplestressR documentation built on Oct. 26, 2021, 9:06 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.
test_that("effect_size_additive works as expected", {
expect_equal(effect_size_additive(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.05)$Interaction_Effect_Size[1], -0.4949691273)
expect_equal(effect_size_additive(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.10)$Interaction_CI_Lower[1], -2.0407155118)
expect_equal(effect_size_additive(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.75)$Null_Model[1], "Additive")
expect_equal(effect_size_additive(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.075,
Small_Sample_Correction = FALSE)$Interaction_Effect_Size[1], -0.5287170223)
expect_equal(effect_size_additive(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.01,
Small_Sample_Correction = FALSE)$Interaction_CI_Upper[1], 2.0583386246)
expect_equal(effect_size_additive(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.50,
Small_Sample_Correction = FALSE)$Null_Model[1], "Additive")
})
test_that("effect_size_multiplicative works as expected", {
expect_equal(effect_size_multiplicative(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.05)$Interaction_Effect_Size[1], -0.0705007812)
expect_equal(effect_size_multiplicative(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.05)$Interaction_CI_Lower[1], -0.2353032640)
expect_equal(effect_size_multiplicative(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.75)$Null_Model[1], "Multiplicative")
})
test_that("critical_effect_size_additive works as expected", {
expect_equal(critical_effect_size_additive(Control_N = 4,
StressorA_N = 4,
StressorB_N = 4,
StressorsAB_N = 4), 1.9397344512)
expect_equal(critical_effect_size_additive(Control_N = 3,
StressorA_N = 3,
StressorB_N = 3,
StressorsAB_N = 3,
Small_Sample_Correction = FALSE,
Significance_Level = 0.10), 2.0163641416)
expect_equal(critical_effect_size_additive(Control_N = 46,
StressorA_N = 23,
StressorB_N = 73,
StressorsAB_N = 4,
Small_Sample_Correction = TRUE,
Significance_Level = 0.50), 0.3850781858)
})
test_that("classify_interactions works as expected", {
expect_equal(classify_interactions(effect_size_additive(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.05),
assign_reversals = TRUE,
remove_directionality = FALSE)$Interaction_Classification[1], "Null")
expect_equal(classify_interactions(effect_size_additive(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.5),
assign_reversals = FALSE,
remove_directionality = TRUE)$Interaction_Classification[1], "Null")
expect_equal(classify_interactions(effect_size_multiplicative(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.01),
assign_reversals = TRUE,
remove_directionality = TRUE)$Interaction_Classification[1], "Null")
expect_equal(classify_interactions(effect_size_multiplicative(Control_N = multiplestressR::survival$Sample_Size_Control[1],
Control_SD = multiplestressR::survival$Standard_Deviation_Control[1],
Control_Mean = multiplestressR::survival$Mean_Control[1],
StressorA_N = multiplestressR::survival$Sample_Size_Temperature[1],
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature[1],
StressorA_Mean = multiplestressR::survival$Mean_Temperature[1],
StressorB_N = multiplestressR::survival$Sample_Size_pH[1],
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH[1],
StressorB_Mean = multiplestressR::survival$Mean_pH[1],
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH[1],
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH[1],
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH[1],
Significance_Level = 0.005),
assign_reversals = FALSE,
remove_directionality = FALSE)$Interaction_Classification[1], "Null")
})
test_that("summary_plots works as expected", {
expect_equal(summary_plots(classify_interactions(effect_size_additive(Control_N = multiplestressR::survival$Sample_Size_Control,
Control_SD = multiplestressR::survival$Standard_Deviation_Control,
Control_Mean = multiplestressR::survival$Mean_Control,
StressorA_N = multiplestressR::survival$Sample_Size_Temperature,
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature,
StressorA_Mean = multiplestressR::survival$Mean_Temperature,
StressorB_N = multiplestressR::survival$Sample_Size_pH,
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH,
StressorB_Mean = multiplestressR::survival$Mean_pH,
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH,
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH,
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH,
Significance_Level = 0.05),
assign_reversals = TRUE,
remove_directionality = FALSE),
Significance_Level = 0.05)[[1]]$data[1,3], 0.764)
expect_equal(summary_plots(classify_interactions(effect_size_additive(Control_N = multiplestressR::survival$Sample_Size_Control,
Control_SD = multiplestressR::survival$Standard_Deviation_Control,
Control_Mean = multiplestressR::survival$Mean_Control,
StressorA_N = multiplestressR::survival$Sample_Size_Temperature,
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature,
StressorA_Mean = multiplestressR::survival$Mean_Temperature,
StressorB_N = multiplestressR::survival$Sample_Size_pH,
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH,
StressorB_Mean = multiplestressR::survival$Mean_pH,
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH,
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH,
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH,
Significance_Level = 0.05),
assign_reversals = TRUE,
remove_directionality = FALSE),
Significance_Level = 0.05)[[2]]$data$Sample_Size_Median[123], 6)
expect_equal(summary_plots(classify_interactions(effect_size_additive(Control_N = multiplestressR::survival$Sample_Size_Control,
Control_SD = multiplestressR::survival$Standard_Deviation_Control,
Control_Mean = multiplestressR::survival$Mean_Control,
StressorA_N = multiplestressR::survival$Sample_Size_Temperature,
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature,
StressorA_Mean = multiplestressR::survival$Mean_Temperature,
StressorB_N = multiplestressR::survival$Sample_Size_pH,
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH,
StressorB_Mean = multiplestressR::survival$Mean_pH,
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH,
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH,
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH,
Significance_Level = 0.05),
assign_reversals = TRUE,
remove_directionality = FALSE),
Significance_Level = 0.05)[[3]]$labels$y, "Density")
expect_equal(summary_plots(classify_interactions(effect_size_additive(Control_N = multiplestressR::survival$Sample_Size_Control,
Control_SD = multiplestressR::survival$Standard_Deviation_Control,
Control_Mean = multiplestressR::survival$Mean_Control,
StressorA_N = multiplestressR::survival$Sample_Size_Temperature,
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature,
StressorA_Mean = multiplestressR::survival$Mean_Temperature,
StressorB_N = multiplestressR::survival$Sample_Size_pH,
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH,
StressorB_Mean = multiplestressR::survival$Mean_pH,
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH,
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH,
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH,
Significance_Level = 0.05),
assign_reversals = TRUE,
remove_directionality = FALSE),
Significance_Level = 0.05)[[4]]$labels$y, "Inverse Variance")
expect_equal(summary_plots(classify_interactions(effect_size_multiplicative(Control_N = multiplestressR::survival$Sample_Size_Control,
Control_SD = multiplestressR::survival$Standard_Deviation_Control,
Control_Mean = multiplestressR::survival$Mean_Control,
StressorA_N = multiplestressR::survival$Sample_Size_Temperature,
StressorA_SD = multiplestressR::survival$Standard_Deviation_Temperature,
StressorA_Mean = multiplestressR::survival$Mean_Temperature,
StressorB_N = multiplestressR::survival$Sample_Size_pH,
StressorB_SD = multiplestressR::survival$Standard_Deviation_pH,
StressorB_Mean = multiplestressR::survival$Mean_pH,
StressorsAB_N = multiplestressR::survival$Sample_Size_Temperature_pH,
StressorsAB_SD = multiplestressR::survival$Standard_Deviation_Temperature_pH,
StressorsAB_Mean = multiplestressR::survival$Mean_Temperature_pH,
Significance_Level = 0.05),
assign_reversals = TRUE,
remove_directionality = FALSE),
Significance_Level = 0.05)[[5]]$labels$y, "Standard Error")
})
Try the multiplestressR 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.