On 11.11.2011 18:48, Gyanendra Pokharel wrote:
Following is my code, can some one help on the error at the bottom?
mh<-function(iterations,alpha,beta){
+ data<-read.table("epidemic.txt",header = TRUE)
+ attach(data, warn.conflicts = F)
+ k<-97
+ d<- (sqrt((x-x[k])^2 + (y-y[k])^2))
+ p<- 1-exp(-alpha*d^(-beta))
+ p.alpha<-1 - exp(-3*d^(-beta))
+ p.beta<- 1 - exp(alpha*d^(-2))
+ iterations<-1000
+ mu.lambda<- c(0,0);s.lambda<- c(100,100)
+ prop.s<- c(0.1,0.1)
+ lambda<- matrix(nrow=iterations, ncol=2)
+ acc.prob<-0
+ current.lambda<- c(0,0)
+ for(t in 1:iterations){
+ prop.lambda<- rnorm(2,current.lambda,prop.s)
+ a<- p.beta/p.alpha
*(dnorm(prop.lambda,mu.lambda,s.lambda))*dnorm(current.lambda,mu.lambda,s.lambda)
+ accept<- min(1,a)
+ u<-runif(1)
+ if(u[t]<=accept[t]){
+ current.lambda<- prop.lambda
+ acc.prob<- acc.prob +1
+ }
+ lambda[t,]<- current.lambda
+ }
+ lambda
+ }
mh(1000,0,0)
Error in if (u[t]<= accept[t]) { : missing value where TRUE/FALSE needed
As the error message tells you:
accept[t] is NA t > 1 since
accept<- min(1,a)
Uwe Ligges
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
