Displaying 3 results from an estimated 3 matches for "f4q".
Did you mean:
f4
2012 Sep 25
1
REML - quasipoisson
...t;-Dpq/(2*phiq)
#2nd term undefined
#3rd term
Xq<-predict(m3,type="lpmatrix")
Wq<-diag(fitted(m3))
Hq<-t(Xq)%*%Wq%*%Xq
ldhsq<-determinant(Hq+Sq,log=TRUE)$modulus[1]
eigSq<-eigen(Sq,only.values=TRUE)$val
ldsq<-sum(log(eigSq[1:16]))
F3q<-(ldhsq-ldsq)/2
#4th term
Mp=3
F4q<-Mp/2*log(2*pi*phiq)
F1q+F3q-F4q
m3$gcv
#quite different
#but if phiq is replaced by the Pearson estimate of the scale
P<-sum((y-fitted(m3))^2/fitted(m3))
phip<-P/(500-sum(m3$edf))
F1p<-Dpq/(2*phip)
F3p<-F3q
#third term independent of scale
F4p<-Mp/2*log(2*pi*phip)
F1p+F3p-F4p...
2012 Oct 01
0
[Fwd: REML - quasipoisson]
...q<-predict(m3,type="lpmatrix")
> Wq<-diag(fitted(m3))
> Hq<-t(Xq)%*%Wq%*%Xq
> ldhsq<-determinant(Hq+Sq,log=TRUE)$modulus[1]
> eigSq<-eigen(Sq,only.values=TRUE)$val
> ldsq<-sum(log(eigSq[1:16]))
> F3q<-(ldhsq-ldsq)/2
>
> #4th term
> Mp=3
> F4q<-Mp/2*log(2*pi*phiq)
>
> F1q+F3q-F4q
> m3$gcv
>
> #quite different
>
> #but if phiq is replaced by the Pearson estimate of the scale
> P<-sum((y-fitted(m3))^2/fitted(m3))
> phip<-P/(500-sum(m3$edf))
> F1p<-Dpq/(2*phip)
> F3p<-F3q
> #third term inde...
2009 Jul 23
1
[PATCH server] changes required for fedora rawhide inclusion.
...6Qd>!-#x{R
zW=%JkBu_w at fXW1!-qkc-P9nCbobn&LDbv=~Ud}$T<%FCsO^sajWJAogHwtpjSFqVG
znYK#$@rrbee-HCB{&l}Q#FO5&mH$NanQi$p!w^sF%qIl-xGVpOQ?2VdswLqr18XJ2
zhI0!i&yWa@{x3gr`u}6^|1$R0tXrKzUQ9gEkeX{Q5#&LBxp}hw*5oNyz?-c}bh*3%
zm}*SHt*oWf<2SjrYI)~J4QGeqZDDJQwQp#y94tAG_~5LhLSs&~!F4Q!J*dk7jJbv5
zr0f8~^hn;F5iilnc=+^US<#HVo&C|b*8`|xG?_lXXT(VB2EBVRon^S9>VV()^&Zv+
zfS413+*JFc`Z0N>`k!U^@Nci**01Y2!DAws1+QM$EuQn6 at 9ICZyThyEX@p31F-w3b
zSp;{Cmu1kte3$6V?kb8$(-YM27+?<CTY}JD0InHUQ at -GR)qtg#qw7n at So}7r5Wr;?
zp>5(awY>wF{bCO75F^>7Q7qt<1+d;kD-Usmb>}{1N1...