[R] foreach loop - rejection method

[dunno87]

Hi Everybody,

Thanks in advance for your help.
This is my first time using the foreach statement and I cant get it to work
properly so here is what i have

test<-function(){
repeat {
cand2[l-1]<-rinvgamma(1,phi,lambda[l-1])
q2<-dinvgamma(cand2[l-1],phi,lambda[l-1])
p2<-cand2[l-1]^-1.5*exp(-y[l]^2/(2*cand2[l-1]))*exp(-((log(cand2[l-1])-mu_t_cand[l-1])^2)/2*sigmasq)
ratio<-p2/(c[l-1]*q2)
if (runif(1)< ratio) {break}
}
f2[l-1]<-min(p2,c[l-1]*q2)
}}
foreach(l=2:(n-1),.combine=c,.packages='MCMCpack',.options.nws=li
st(chunkSize=250)) %dopar% test()

even though i've created the array earlier called cand2 when i run this code
i get the error object "cand2" not found.

What I'm trying to do is the rejection method for my cand2. i need cand2 and
f2 as vectors as output both of length l-2 in the correct order.

if i do a for loop as normal it works perfectly. Any ideas?

Ben
-- 
View this message in context: http://www.nabble.com/foreach-loop---rejection-method-tp25796927p25796927.html
Sent from the R help mailing list archive at Nabble.com.