Displaying 2 results from an estimated 2 matches for "new_u".
Did you mean:
new_p
2011 Mar 25
0
Bounding ellipse for any set of points
...stions/1768197/bounding-ellipse
## for commented code
X <- Q %*% diag(u) %*% t(Q)
M <- diag(t(Q) %*% solve(X) %*% Q)
maximum <- max(M)
j <- which(M == maximum)
step_size = (maximum - d -1) / ((d+1)*(maximum-1))
new_u <- (1 - step_size) * u
new_u[j] <- new_u[j] + step_size
err <- sqrt(sum((new_u - u)^2))
count <- count + 1
if (count > max.iter) {
warning(paste("Iterated", max.iter, "times and still can't find
the boun...
2011 Jul 21
0
Bounding ellipse for any set of points
...; ## for commented code
> X <- Q %*% diag(u) %*% t(Q)
> M <- diag(t(Q) %*% solve(X) %*% Q)
> maximum <- max(M)
> j <- which(M == maximum)
> step_size = (maximum - d -1) / ((d+1)*(maximum-1))
> new_u <- (1 - step_size) * u
> new_u[j] <- new_u[j] + step_size
> err <- sqrt(sum((new_u - u)^2))
> count <- count + 1
> if (count > max.iter) {
> warning(paste("Iterated", max.iter, "times and still c...