Delete fisher_mix.R

R/fisher_mix.R

-##### Simulations of Fisher mixing
-###############################################################################
-#' Titre
-#'
-#' Description
-#'
-#' Long description (several lines)
-#'
-#' @param nsim number of simulations to perform
-#' @param p number of components in the mixture (number of variables)
-#' @param nobs total number of observations in the sample
-#'
-#' @return a vector with simulated random variables fitting the Fisher mixing
-#' distribution
-#'
-#' @seealso \code{\link{treediff}} to perform the test based on these
-#' simulations
-#' 
-#' @references To fulfill
-#'
-#' @examples
-#' x <- fisher_mix(1000, 500)
-#' \dontrun{hist(x, breaks = 100)}
-#'
-#' @export
-#' 
-#' @importFrom stats rf
-
-fisher_mix <- function(nsim, p, nobs) {
- all_fishers <- replicate(nsim, sum(rf(p, 1, 1 * (nobs - 2))) / p)
- 
- return(all_fishers)
-}