Displaying 2 results from an estimated 2 matches for "ldhs".
Did you mean:
dhs
2012 Sep 25
1
REML - quasipoisson
...{
S[i,j]=S1[i-1,j-1]
S[i+9,j+9]=S2[i-1,j-1]
}
}
beta<-m1$coef
#penalised deviance
Dp<-m1$dev+t(beta)%*%S%*%beta
F1<-Dp/(2*phi)
#2nd term
F2<-sum(ifelse(y==0,0,y*log(y)-y-log(factorial(y))))
#3rd term
X<-predict(m1,type="lpmatrix")
W<-diag(fitted(m1))
H<-t(X)%*%W%*%X
ldhs<-determinant(H+S,log=TRUE)$modulus[1]
eigS<-eigen(S,only.values=TRUE)$val
lds<-sum(log(eigS[1:16]))
F3<-(ldhs-lds)/2
#4th term
Mp=3
F4<-Mp/2*log(2*pi*phi)
F1-F2+F3-F4
m1$gcv
#reml score = formula
##quasipoisson with scale = 1
#fitting is identical to the poisson case
#F1, F3, and...
2012 Oct 01
0
[Fwd: REML - quasipoisson]
...1$coef
> #penalised deviance
> Dp<-m1$dev+t(beta)%*%S%*%beta
> F1<-Dp/(2*phi)
>
> #2nd term
> F2<-sum(ifelse(y==0,0,y*log(y)-y-log(factorial(y))))
>
> #3rd term
> X<-predict(m1,type="lpmatrix")
> W<-diag(fitted(m1))
> H<-t(X)%*%W%*%X
> ldhs<-determinant(H+S,log=TRUE)$modulus[1]
> eigS<-eigen(S,only.values=TRUE)$val
> lds<-sum(log(eigS[1:16]))
> F3<-(ldhs-lds)/2
>
> #4th term
> Mp=3
> F4<-Mp/2*log(2*pi*phi)
>
> F1-F2+F3-F4
> m1$gcv
> #reml score = formula
>
> ##quasipoisson with sca...