Displaying 2 results from an estimated 2 matches for "modelexpress".
2007 Mar 13
0
segfault with correlation structures in nlme
...fm10var.nlme,corr=corAR1(0.311))
>fm30var.nlme <- update(fm10var.nlme,corr=corARMA(p=0,q=2))
>fm60var.nlme <- update(fm10var.nlme,corr=corGaus(form=~Time))
*** caught segfault ***
address 0x1075e501, cause 'memory not mapped'
Traceback:
1: eval(expr, envir, enclos)
2: eval(modelExpression, env)
3: assign("modelValue", eval(modelExpression, env), envir = thisEnv)
4: function (newPars) { if (!missing(newPars)) { for (i in
names(ind)) { assign(i, clearNames(newPars[ind[[i]]]), envir
= env) } assign("modelValue", eval(modelExpres...
2001 Dec 03
0
problems with nmle
...date( fm1Indom.nlme,
random = pdDiag(A1 + lrc1 + A2 ~ 1) )
debugging nlme gives the following output:
Browse[1]> n
debug: modelResid <- ~eval(model, data.frame(data, getParsNlme(plist,
fmap, rmapRel, bmap, groups, beta, bvec, b, level, N)))[naPat]
Browse[1]> n
debug: ww <- eval(modelExpression[[2]], envir = nlEnv)
Browse[1]> n
Error: subscript out of bounds
Here the fm1Indom.nlme object is calculated with a formula somewhat
different from the one given
in the text, because that gives convergence problems:
fm1Indom.nlme <- nlme(fm1Indom.lis, random=pdDiag(A1+lrc1+A2+lrc2 ~...