Displaying 1 result from an estimated 1 matches for "sigsvd".
Did you mean:
signed
2006 Nov 30
1
data.frame within a function (PR#9294) (cont'd)
...are matrix")
}
if (length(mean) != nrow(sigma)) {
stop("mean and sigma have non-conforming size")
}
ev <- eigen(sigma, sym = TRUE)$values
if (!all(ev >= -sqrt(.Machine$double.eps) * abs(ev[1])))
warning("sigma is numerically not positive definite")
sigsvd <- svd(sigma)
retval <- t(sigsvd$v %*% (t(sigsvd$u) * sqrt(sigsvd$d)))
retval <- matrix(rnorm(n * ncol(sigma)), nrow = n) %*% retval
retval <- sweep(retval, 2, mean, "+")
retval
}
f<- function(obj){
update(obj,~ .+x,evaluate=T); cat("also ok\n")
}
#...