Displaying 2 results from an estimated 2 matches for "wcor".
Did you mean:
cor
2010 Jun 17
0
Modifyiing R working matrix within "gee" source code
...e(N)))
N <- as.double(N)
modvec <- as.integer(c(linkv, varfunv, corstrv))
if (v4.4compat)
compatflag <- 1
else compatflag <- 0
z <- .C("Cgee", x, y, id, N, offset, nobs, p, modvec, Mv,
estb = beta, nv = naivvar, rv = robvar, sc = phi, wcor = R,
tol, mc = maxcl, iter = maxiter, silent, err = errstate,
scale.fix, as.integer(compatflag), PACKAGE = "gee")
if (z$err != 0)
warning("Cgee had an error (code= ", z$err, "). Results suspect.")
if (min(eigen(z$wcor)$values) < 0) {...
2004 Dec 10
0
strange gee behavior
...I have a problem where the script locks up R when calling this function (it
starts execution and never finishes). I was able to track down the
offending data and determine that this particular dataset results in a
correlation matrix that is not positive definite and a line within gee with
eigen(z$wcor) is where things lock up . Specifically, the line within eigen
.Call("La_rs", x, only.values, if (capabilities("IEEE754")) "dsyevr" else
"dsyev", PACKAGE = "base")
is where it happens.
The issue is that I get different behavior when I run the fun...