R/phi.R

In leonardobfn/ThesiR:

phi <- function(zm,alfastart,BB,TT=TT,Lmm=Lmm){
  #BB=B;zm=z
  CB   <-  ((sin(pi*alfastart*BB)/sin(pi*BB))^(1/(1-alfastart)))*((sin(pi*(1-alfastart)*BB))/sin(pi*alfastart*BB))
  #print(CB)
  aux_CB <- CB*exp(-CB*zm^(-alfastart/(1-alfastart)))

  # aux_CB[aux_CB==0] <- 1

  phi0 <- mean(aux_CB,na.rm=F)
  phi1 <- mean(log(CB)*aux_CB,na.rm=F)
  phi2 <- mean(CB*aux_CB,na.rm=F)
  #phiB <- mean(BB*aux_CB)
  return(data.frame(phi0,phi1,phi2))
}