R/Str.R

In xinghuq/HierD: R package for calculating hierarchical diversity across metrixs and scales, from genes to ecosystem

#### Deifine hierarchical structure for genetic data

Str=function(nreg,r,npops){
  Str=data.frame(matrix(data=0,ncol=npops, nrow=3)) # 
  Str[1,]=c(rep("ecosystem",times=npops))
  # Str[2,]=c(rep(paste0("region",1:nreg),each=(npops/nreg)))
  #r=vector(r) # eg. c(r1,r2,r3,r4) number of elements should be equal to nreg
  if (length(r)!= nreg) stop("number of regions should be equal to the number of  r  defined")
  Str[2,]=c(rep(paste0("region",1:nreg),r))
  Str[3,]=c(paste("pop",1:npops))
  Str=as.matrix(Str)
  return(Str)
}