[S] Sym Sq Root Errata

Erin Hodgess (hodgess@uhdux2.dt.uh.edu)
Mon, 18 May 1998 13:44:44 -0500

Dear S+ People
in line 3 of the previous function, there is a typo!!!!!

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

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