tests/ras.R

tab <- matrix(round(runif(12, 0, 2)), 3, 4)

library(logmult)
data(gss8590)

tab <- gss8590[,,1]
row <- c(1, 1, 1, 1)
col <- c(1, 1, 1, 1, 1) * 4/5
tab2 <- ras(tab, row, col)
lambda1 <- logmult:::lambda(tab)
lambda2 <- logmult:::lambda(tab2)
stopifnot(all.equal(lambda1, lambda2))
nalimilan/logmult documentation built on March 28, 2022, 1:18 p.m.