[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