Displaying 2 results from an estimated 2 matches for "thnew".
Did you mean:
t_new
2008 May 01
4
efficient code - yet another question
...- FALSE
while (!ende)
{
nr <- nr + 1
ph <- t(crossprod(th, Xh) * as.vector(1 /
crossprod(th)))
ph <- ph * as.vector(1/sqrt(crossprod(ph)))
thnew <- t(tcrossprod(t(ph), Xh) *
as.vector(1/(crossprod(ph))))
prec <- crossprod(th-thnew)
th <- thnew
if (prec <= (tol^2)) ende <- TRUE
if (it <= nr) ende <- TRUE # didn'...
2008 Oct 15
0
R-help Digest, Vol 67, Issue 31
...while (!ende)
> {
> nr <- nr + 1
> ph <- t(crossprod(th, Xh) * as.vector(1 /
> crossprod(th)))
> ph <- ph * as.vector(1/sqrt(crossprod(ph)))
> thnew <- t(tcrossprod(t(ph), Xh) *
> as.vector(1/(crossprod(ph))))
> prec <- crossprod(th-thnew)
> th <- thnew
> if (prec <= (tol^2)) ende <- TRUE
> if (it <= nr) ende...