R/Hbeta2.fun.R

Hbeta2.fun <-
function(ord.z,ord.bz,Rt,r,p,ord.wt){
    a = matrix(0,nrow=p,ncol=p)
    for(j in 1:p){
      for(k in 1:p){
         a[j,k] = -sum(ord.wt*as.numeric(ord.z[,j])*as.numeric(ord.z[,k])*Rt*exp(-ord.bz)/(1+r*Rt*exp(-ord.bz)))
      }
    }
    return(a)
}

Try the TransModel package in your browser

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

TransModel documentation built on May 2, 2019, 6:13 a.m.