polarolouis/anova-phylogenetique-projet-msv
1
0
Fork 0
mirror of https://github.com/Polarolouis/anova-phylogenetique-projet-msv.git synced 2026-06-17 10:15:25 +02:00
Code Issues Projects Releases Packages Wiki Activity

🐛Fixing bugs that lead to incorrect type I error evaluation

Browse source
This commit is contained in:
Louis Lacoste 2024-01-23 18:52:08 +01:00
parent 37e41bb2fc
commit 3908729ac2
1 changed files with 6 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

8
R/anovaComparison.R
View file

@ -142,6 +142,10 @@ simulate_matching_and_random <- function(
stoch_process, test_method,
risk_threshold = 0.05,
correct_hypothesis = "H1") {
# To be in the right configuration for typeI error
if (correct_hypothesis == "H0") {
base_values <- rep(0, length(base_values))
}
matching_phylo_traits <- compute_trait_values(
groups = phylo_matching_groups,
base_values = base_values, tree,
@ -165,7 +169,7 @@ simulate_matching_and_random <- function(
random_groups_pvalues <- phyloanova_anova_pvalues(
traits = random_groups_traits,
groups = random_groups, tree, stoch_process = stoch_process,
test_method = test_method, measurement_error = (sigma2_measure != 0)
test_method = test_method, measurement_error = TRUE
)
# Concatenate pvalues
@ -207,7 +211,7 @@ simulate_data <- function(
stoch_process = stoch_process,
test_method = test_method,
risk_threshold = risk_threshold,
correct_hypothesis = "H1"
correct_hypothesis = correct_hypothesis
)
}))