R/p.estFUN.R

p.estFUN <- function(jk, x, y, n){
  j = jk[1]
  k = jk[2]
  a <- p.est(x,y,n,j,k)
  s2 <- a$sigma2
  t2 <- a$tau2
  return(p.ll(n, j, k, s2, t2))
}

Try the hcp package in your browser

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

hcp documentation built on May 2, 2019, 2:37 a.m.