R/hyper.sigcor.R

In hyper.fit: Generic N-Dimensional Hyperplane Fitting with Heteroscedastic Covariant Errors and Intrinsic Scatter

hyper.sigcor=function(N,parmDF){
  if(N<340){
    bias2unbias=(sqrt((N-parmDF)/2)*gamma((N-parmDF)/2))/gamma((N-parmDF+1)/2)
  }else{
    bias2unbias=1
  }
  samp2pop=sqrt(N/(N-parmDF))
  return=c(bias2unbias=bias2unbias,samp2pop=samp2pop,sampbias2popunbias=bias2unbias*samp2pop)
}