R/aux3.R

aux3 <- function(m,ind)
{
#
# Auxilliary function to re-structure the information matrix
# in the case where the variances of all components are constrained
# to be equal.
#
x <- apply(m[,ind],1,sum)
s <- sum(x[ind])
y <- x[-ind]
r <- m[-ind,-ind]
rbind(cbind(r,y),c(y,s))
}

Try the mixreg package in your browser

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

mixreg documentation built on May 2, 2019, 3:25 a.m.