[Rcpp-devel] Help for call R function
bbonit at tin.it
bbonit at tin.it
Fri Jun 15 22:04:01 CEST 2012
Hello to everyboby Im Gianluca ....and i m a statistic student.
I have one problem that i think Dirk or some one "good like dirk" could help to resolve :
iI would like to implement this following function in R code into inline code with Rccp plugin :
MHWG<-function (logpost, start, n.iter, burn, thin, scale, ...)
{
pb <- txtProgressBar(min = 0, max = n.iter, style = 3)
p = length(start)
vth = array(0, dim = c(n.iter/thin, p))
f0 = logpost(start, ...)
arate = array(0, dim = c(1, p))
th0 = start
for (i in (-burn):n.iter) {
setTxtProgressBar(pb, i)
for (j in 1:p) {
th1 = th0
th1[j] = th0[j] + rnorm(1) * scale[j]
f1 = logpost(th1, ...)
u = runif(1) <= min(exp(f1 - f0),1)
th0[j] = th1[j] * (u == 1) + th0[j] * (u == 0)
f0 = f1 * (u == 1) + f0 * (u == 0)
vth[floor(i/thin), j] = th0[j]
arate[j] = arate[j] + u
}
}
arate = arate/n.iter
stuff = list(par = vth, accept = arate)
return(stuff)
}
So one argument of MHWG is logpost and is a R function with argument par where par iso ne vector .... for example logpost(par) Logpost<-function(par) { ...................} (in r code).
How can i call teh logpost function into the inline rccp code . In the MHWG there is a difference f1-f0 where fi=logpost(par1) e f0=logpost(par 2) for exemple.
help me please to implemet my problem ...and also i m "open" about suggestion for implemeting the function
Sorry for all noise and thank also for the past mail that help me to learn more about programming
Best regards... sorrry fo rmy bad English
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/rcpp-devel/attachments/20120615/7892c407/attachment.html>
More information about the Rcpp-devel
mailing list