R/rhoc.R

In meszre/tauPFC: Computes Robust Estimators for the PFC Model

rhoc=function (x,c)
{
  G0=1.38
  G1=0.55
  G2=-2.69
  G3=10.76
  G4=-11.66
  G5=4.04
  u =(x>c)
  v=(x<2/3*c)
  w=(1-u)*(1-v)
  v*G0*(x/c)^2+ w*(G1+G2*(x/c)^2+G3*(x/c)^4+G4*(x/c)^6+ G5*(x/c)^8)+ 1*u
}