lawrence: Compute the Lawrence lifting of a configuration matrix

In dkahle/algstat: Algebraic Statistics

Compute the Lawrence lifting of a configuration matrix

Description

Compute the Lawrence lifting of a configuration matrix

Usage

lawrence(A, k = 1, extent = c("short", "tall"))

Arguments

A	the configuration matrix of the model
k	kth Lawrence lifting
extent	short or tall version

Value

an integer matrix

Author(s)

David Kahle [e-mail]

See Also

Aoki, Hara, and Takemura (2012), pp. 61-62

Examples

(A <- hmat(c(2,2), 1:2))
plot_matrix(A) # in package latte

lawrence(A)
plot_matrix(lawrence(A))

lawrence(A, 2)
plot_matrix(lawrence(A, 2))

lawrence(A, extent = "tall")
plot_matrix(lawrence(A, extent = "tall"))

plot_matrix(A)
plot_matrix(lawrence(A))
plot_matrix(lawrence(A, 2))
plot_matrix(lawrence(A, 3))
plot_matrix(lawrence(A, 20))

plot_matrix(lawrence(A))
plot_matrix(lawrence(A, extent = "tall"))


## Not run:  # requires 4ti2

J <- 10 # number of levels
A <- rbind(1L, 1:J)
markov(A)

lawrence(A)
plot_matrix(lawrence(A))
markov(lawrence(A))
# markov(lawrence(A, 2)) # takes unknown long time when J=10
zbasis(lawrence(A))


## End(Not run)

dkahle/algstat documentation built on May 23, 2023, 12:29 a.m.