Extracting nb_species and K number of groups from model fit

This commit is contained in:
Louis Lacoste 2024-02-07 11:43:50 +01:00
parent aa13972a3d
commit 00d4f979f8

View file

@ -231,7 +231,6 @@ compute_trait_values <- function(
return(trait + trait_phylo + trait_error)
}
# TODO Récupérer K et nb_species des objets
#' Infere an ANOVA and a phyloanova
#'
#' @param y the vector of traits for which to fit the models
@ -263,6 +262,10 @@ pvalues_from_fits <- function(
#  For sanity test
match.arg(tested_method)
# Extracting values
nb_species <- nrow(model.frame(fit_anova))
K <- length(unique(model.frame(fit_anova)$groups))
is_invalid_value <- function(value) {
return(is.nan(value) ||
is.null(value) ||