R/selsscale2.R

Defines functions selsscale2

Documented in selsscale2

# Sigma scaling 2
selsscale2 = function(fitvals, ns, selc, ...){
  n = length(fitvals)
  if(missing(ns)) ns=n
  if(missing(selc)) selc=2
  if(any(fitvals<0)) fitvals=fitvals-min(fitvals)
  favg = mean(fitvals)
  fsd = sd(fitvals)
  fstar = fitvals + (favg-selc*fsd)
  fstar[fstar<0]=0
  p = fstar/sum(fstar)
  matpool = sample(1:n, size=ns, prob=p, replace=TRUE)
  return(matpool)
}

Try the adana package in your browser

Any scripts or data that you put into this service are public.

adana documentation built on March 18, 2022, 6:03 p.m.