R/genInits.R

In statdivlab/happi: a hierarchical approach to pangenomics inference

#' Function for generating initial starts for estimated beta
#'
#' @param num_covariate number of covariates p (includes intercept)
#' @param nstarts number of starts
#' @param seed numeric seed for random initializations
#' @param norm_sd standard deviation for the Normal distribution used to generate random initializations
#'
#' @return matrix of initializations
#'
#' @examples
#' genInits(num_covariate = 2, nstarts = 1, seed = 88)
#' @export
genInits <- function(num_covariate, 
                     nstarts = 1, 
                     seed,
                     norm_sd) {
  
  init_start <- rep(0,num_covariate) # fixed initial start parameters
  inits <- rbind(init_start)

if (nstarts > 1) { 
    set.seed(seed)
    for (i in 2:nstarts) {
    inits <- rbind(inits, stats::rnorm(length(init_start), init_start, norm_sd))
    # make sure to set.seed() before running multiple starts
    } 
  }
  
  return(inits)
}