search for: lamtut

...eturns the external pointer, regenerating if necessary' if (length(theta)) { if (.Call(isNullExtPtr, Ptr)) initializePtr() } Ptr }, ## ditto initializePtr = function() { Ptr <<- .Call(merPredDCreate, as(X, "matrix"), Lambdat, LamtUt, Lind, RZX, Ut, Utr, V, VtV, Vtr, Xwts, Zt, beta0, delb, delu, theta, u0) ... } merPredDCreate in turn just copies the relevant bits into a new C++ class object: /* see src/external.cpp */ SEXP merPredDCreate(SEXP Xs, SEXP Lambdat, SEXP LamtUt, SEXP Lind,...

.....@ i : int [1:4] 0 1 2 3 .. .. .. ..@ p : int [1:5] 0 1 2 3 4 .. .. .. ..@ Dim : int [1:2] 4 4 .. .. .. ..@ Dimnames:List of 2 .. .. .. .. ..$ : NULL .. .. .. .. ..$ : NULL .. .. .. ..@ x : num [1:4] 0.499 0.499 0.499 0.499 .. .. .. ..@ factors : list() .. ..$ LamtUt :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots .. .. .. ..@ i : int [1:20] 0 0 0 0 0 1 1 1 1 1 ... .. .. .. ..@ p : int [1:21] 0 1 2 3 4 5 6 7 8 9 ... .. .. .. ..@ Dim : int [1:2] 4 20 .. .. .. ..@ Dimnames:List of 2 .. .. .. .. ..$ : NULL .....

Javier, Creo que aquí aplica la ley de Linus que dice: "Dado un número suficientemente elevado de ojos, todos los errores se convierten en obvios". La persona que revisa y encuentra un error no necesariamente tiene que ser la misma que la que lo escribe. Una motivación muy importante al compartir un código es la de recibir los beneficios del control de calidad por parte de tus pares.

.... .. .. ..@ p : int [1:5] 0 1 2 3 4 > .. .. .. ..@ Dim : int [1:2] 4 4 > .. .. .. ..@ Dimnames:List of 2 > .. .. .. .. ..$ : NULL > .. .. .. .. ..$ : NULL > .. .. .. ..@ x : num [1:4] 0.499 0.499 0.499 0.499 > .. .. .. ..@ factors : list() > .. ..$ LamtUt :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots > .. .. .. ..@ i : int [1:20] 0 0 0 0 0 1 1 1 1 1 ... > .. .. .. ..@ p : int [1:21] 0 1 2 3 4 5 6 7 8 9 ... > .. .. .. ..@ Dim : int [1:2] 4 20 > .. .. .. ..@ Dimnames:List of 2 > .....