# NOT RUN {
# number of counties:
n <- nrow( UScounties.storder)
# }
# NOT RUN {
# make a precision matrix
Q <- diag.spam( n) + .2 * UScounties.storder + .1 * UScounties.ndorder
display( as.spam( chol( Q)))
# }
Run the code above in your browser using DataLab