[Yuima-commits] r102 - in pkg/yuima: R man
noreply at r-forge.r-project.org
noreply at r-forge.r-project.org
Tue Jul 13 21:54:20 CEST 2010
Author: iacus
Date: 2010-07-13 21:54:20 +0200 (Tue, 13 Jul 2010)
New Revision: 102
Modified:
pkg/yuima/R/qmle.R
pkg/yuima/man/adaBayes.Rd
Log:
fixed examples in adaBayes
Modified: pkg/yuima/R/qmle.R
===================================================================
--- pkg/yuima/R/qmle.R 2010-07-13 19:27:41 UTC (rev 101)
+++ pkg/yuima/R/qmle.R 2010-07-13 19:54:20 UTC (rev 102)
@@ -421,11 +421,12 @@
}
}
}
+
-
- if(QL==-Inf){
+ if(!is.finite(QL)){
yuima.warn("quasi likelihood is too small to calculate.")
+ return(1e10)
}
if(print==TRUE){
yuima.warn(sprintf("NEG-QL: %f, %s", -QL, paste(names(param),param,sep="=",collapse=", ")))
Modified: pkg/yuima/man/adaBayes.Rd
===================================================================
--- pkg/yuima/man/adaBayes.Rd 2010-07-13 19:27:41 UTC (rev 101)
+++ pkg/yuima/man/adaBayes.Rd 2010-07-13 19:54:20 UTC (rev 102)
@@ -60,7 +60,7 @@
bayes1 <- adaBayes(yuima, start=param.init, lower=lower,upper=upper,prior=prior, method="nomcmc")
bayes1
-mle1 <- qmle(yuima, start=param.init, lower=lower, upper=upper, method="L-BFGS-B")
+mle1 <- qmle(yuima, start=param.init, lower=list(theta1=0,theta2=0), upper=list(theta1=1,theta2=1), method="L-BFGS-B")
mle1 at coef
}}
\keyword{ts}
More information about the Yuima-commits
mailing list