search for: settheta

...mer(Reaction ~ Days + (Days|Subject), sleepstudy) pp <- fm1 at pp pp$theta ## [1] 0.96673279 0.01516906 0.23090960 address(pp$theta) ## something pp$Ptr ## <pointer: ...> xpp <- pp$copy() ## default is deep copy xpp$Ptr ## <pointer: (nil)> address(xpp$theta) ## same as above xpp$setTheta(c(0,0,0)) ## referenced through Ptr field xpp$Ptr ## now set to non-nil fm1 at pp$theta ## changes to (0,0,0). oops. So apparently when the xpp$theta object is copied into the external pointer, a reference/lazy copy is made. (xpp$theta itself is read-only, so I can't do the assignment tha...

...15.04 .. ..$ theta : num 0.499 .. ..$ u0 : num [1:4] 0 0 0 0 .. ..and 42 methods, of which 30 are possibly relevant: .. .. b, beta, CcNumer, copy#envRefClass, initialize, initializePtr, installPars, L, ldL2, ldRX2, linPred, P, ptr, RX, RXdiag, RXi, .. .. setBeta0, setDelb, setDelu, setTheta, solve, solveU, sqrL, u, unsc, updateDecomp, updateL, updateLamtUt, updateRes, updateXwts ..@ optinfo:List of 4 .. ..$ optimizer: chr "Nelder_Mead" .. ..$ control :List of 1 .. .. ..$ verbose: int 0 .. ..$ conv : num 0 .. ..$ warnings : list() #########################...

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.

...: num 0.499 > .. ..$ u0 : num [1:4] 0 0 0 0 > .. ..and 42 methods, of which 30 are possibly relevant: > .. .. b, beta, CcNumer, copy#envRefClass, initialize, initializePtr, installPars, L, ldL2, ldRX2, linPred, P, ptr, RX, RXdiag, RXi, > .. .. setBeta0, setDelb, setDelu, setTheta, solve, solveU, sqrL, u, unsc, updateDecomp, updateL, updateLamtUt, updateRes, updateXwts > ..@ optinfo:List of 4 > .. ..$ optimizer: chr "Nelder_Mead" > .. ..$ control :List of 1 > .. .. ..$ verbose: int 0 > .. ..$ conv : num 0 > .. ..$ warnings : list(...