Displaying 2 results from an estimated 2 matches for "ddiag".
Did you mean:
diag
2002 Sep 10
2
Hat values for generalized additive models
...formatting etc. This is my first posting.)
diff -ur c:/usr/temp/work/mgcv/R/mgcv.r ./R/mgcv.r
--- c:/usr/temp/work/mgcv/R/mgcv.r 2002-05-17 03:04:26.000000000 +1000
+++ ./R/mgcv.r 2002-09-10 14:22:06.000000000 +1000
@@ -437,6 +437,7 @@
edf<-array(0,m) # estimated degrees of freedom
ddiag<-array(0,3*m) # array for diagonostics
idiag<-array(0,3) # array for diagnostics
+ hat<-array(0,n) # array for diagnostics
Vp[1,1]<-1.0
M$gcv.ubre<-1.0;
direct.mesh<-100 # number of points for overall s.p. initial direct
search
@@ -448,7 +449,8 @@
a...
2004 Jul 30
1
lme: problems with corARMA
Trying following example from Pinheiro and Bates in order to fit an
ARMA(1,1) model:
library(nlme)
fm1Ovary.lme<-lme(follicles~sin(2*pi*Time)+cos(*pi*Time),data=Ovary,random=p
dDiag(~sin(2*pi*Time)))
fm5Ovary.lme<-update(fm1Ovary.lme,corr=corARMA(p=1,q=1))
I get follwing error message:
Error in "coef<-.corARMA"(`*tmp*`, value = c(62.3428455941166,
62.3428517930051 :
Coefficient matrix not invertible
Does somebody know why it doesn't work ?
Step...