R/con2.parmsFUN.R

con2.parmsFUN<-
  function(jk,x,y,n){
    j = jk[1]
    k = jk[2]
    a <- con2.parms(x,y,n,j,k,1,1)
nr <- nrow(a$theta)
est <- a$theta[nr,  ]
b <- con2.est(x[j], x[k], est)
s2 <- 1/b$eta1
t2 <- 1/b$eta2
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.