#
# Lets look at the Florentine marriage data
#
data(florentine)
#
# test the summary_formula function
#
summary(flomarriage ~ edges + kstar(2))
m <- as.matrix(flomarriage)
summary(m ~ edges) # twice as large as it should be
summary(m ~ edges, directed=FALSE) # Now it's correct
Run the code above in your browser using DataLab