Bmat: An internal function for jrfit
In kloke/jrfit: Joint ranking estimates

Description Usage Arguments Author(s) Examples

An internal function for jrfit.

1	Bmat(v1, v2, block)

`v1`	scalar
`v2`	scalar
`block`	vector

John Kloke <kloke@biostat.wisc.edu>

function (v1, v2, block) 
{
    ublock <- unique(block)
    m <- length(ublock)
    N <- length(block)
    B <- matrix(0, nrow = N, ncol = N)
    ind <- 1
    for (k in 1:m) {
        nk <- sum(block == ublock[k])
        Bk <- diag(rep(v1 - v2, nk)) + v2
        B[ind:(ind + nk - 1), ind:(ind + nk - 1)] <- Bk
        ind <- ind + nk
    }
    B
  }