Displaying 1 result from an estimated 1 matches for "eigomega".
2013 Jan 18
0
problem that arises after using the new version of "BRugs"
...-c(1, (L - 1), L)] + rep(allKnots,
each = 3)[-c(1, 2, L)])/2
wts <- rep(diff(allKnots), each = 3) * rep(c(1, 4, 1)/6,
K + 7)
Bdd <- spline.des(allKnots, xtilde, derivs = rep(2, length(xtilde)),
outer.ok = TRUE)$design
Omega <- t(Bdd * wts) %*% Bdd
eigOmega <- eigen(Omega)
indsZ <- 1:(numIntKnots + 2)
UZ <- eigOmega$vectors[, indsZ]
LZ <- t(t(UZ)/sqrt(eigOmega$values[indsZ]))
indsX <- (numIntKnots + 3):(numIntKnots + 4)
UX <- eigOmega$vectors[, indsX]
L <- cbind(UX, LZ)
stabCheck <- t(crossprod(L, t(...