R/AlphaCross.R

AlphaCross <-
function (L1, R1, n1, est2, k2)
{
  alpha = matrix(0, n1, k2, dimnames=list(1:n1, row.names(est2$alpha)))
  for (i in 1:n1)
  {
    tempint = L1[i] < est2$intmap[2,] & R1[i] >= est2$intmap[2,]
    alpha[i, tempint] = 1
  }
  alpha
}

Try the glrt package in your browser

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

glrt documentation built on May 2, 2019, 8:31 a.m.