Here is the correct function:
dar3 <- function(x)
{
x1 <- eigen(x)
if(any(x1$values < 0))
stop("Negative Eigenvalues")
x1e <- diag(x1$values)
x1e <- sqrt(x1e)
u1 <- cbind(x1$vectors)
u1i <- solve(u1)
w <- u1 %*% x1e %*% u1i
w
}
Sincerely,
Erin Hodgess
Assistant Professor
Dept. of Computer and Mathematical Sciences
University of Houston -- Downtown
e-mail: hodgess@uhdux2.dt.uh.edu
-----------------------------------------------------------------------
This message was distributed by s-news@wubios.wustl.edu. To unsubscribe
send e-mail to s-news-request@wubios.wustl.edu with the BODY of the
message: unsubscribe s-news