similar to: Zero truncated Poisson distribution & R2WinBUGS

Dear All, I am trying to use the do.call command to avoid tedious loops in R when I need to call repetitively a function. I am experiencing a problem when the arguments of a function are given with a vector (I need to express then this way to be able later on to integrate them numerically with the adapt package). Consider the small script: remove(list=ls()) #library(adapt) # first a list of the

Dear All, I'd like to know if it is possible to create animations with R. To be specific, I attach a code I am using for my research to plot some analytical results in 3D using the lattice package. It is not necessary to go through the code. Simply, it plots some 3D density profiles at two different times selected by the user. I wonder if it is possible to use the data generated for different

I'm not aware of any routine that those the job, although I think that it could be relatively easily done by multiplication the manifest variable vector with the estimates for the specific effect. To make an example: v1; v2; v3; v4 are manifest variables that loads on one y latent variablein a data frame called "A" the code for the model should be like: model <-specifymodel( y

hi all, I want to get a cumsum according to the order of some variable. However, it doesnt' work. For example, ********************** test<-data.frame(cbind(x=c(3,5,2,6,7),y=c(8,1,4,9,0))) test[order(test$x),]$sumy<-cumsum(test[order(test$x),]$y) ********************** R complians Warning message: In `[<-.data.frame`(`*tmp*`, order(test$x), , value = list(x = c(2, : provided 3

Hi all, I am new to R. I made this script that plots the deviation of a skellam-distributed random variable with respect to the skellam distribution. I would expect to get random errors, but the plot sistematically shows a non-random pattern (first a peak and then a low). I don't know how to explain this. Can somebody enlighten me?? require(skellam) n <- 5000 lam1 <- 5.45 lam2 <-

Just while I use R2Winbugs, the following error will be presented. Please tell me how to cope with. Error in file(file, "r") : unable to open connection In addition: Warning message: cannot open file 'codaIndex.txt', reason 'No such file or directory' in: file(file, "r") -- View this message in context:

Greetings - I run a multi-pro server box (GNU/Linux), on which I've installed winBUGS under wine. Works fine, and plays nice with r2WinBUGS, provided everything is done through a graphical front-end (either by working from console, or using a virtual desktop via VNC or equivalent). However, I'm wondering if there is a way to use R2winBUGS and winBUGS such that winBUGS does not spawn

Hello ALL! I am running Linux, Fedora 15 64-bits, and R on it. I need to use WinBUGS and R2WinBUGS, but as far as I read, WinBUGS is closed project, to be continued with/as OpenBUGS. Thus, I have found R2OpenBUGS on OpenBUGS Contributed Code (http://openbugs.info/w/UserContributedCode), not on CRAN. Author(s) states that it is equivalent for R2WinBUGS. I tried briefly, and realized few minor

I have an array from an mcmc simulation and I'd like to try and use the monitor function in R2WinBUGS. The R2WinBUGS package is in the search path and when I try help(monitor) the help page appears for the function. But, when I try and use it the following happens: > monitor(sims) Error: couldn't find function "monitor" I've searched the archives but have not seen any

This weekend I noticed that my R2WinBUGS connection was no longer working on my Windows computer at work AND my Ubuntu linux computer at home. As soon as WinBUGS opens, the message reads Index Out Of Range. I have un-installed and re-installed R on both computers and the problem remains the same. At first I thought it was the new version of R, but I was able to overcome the problem (on Windows

Dear R users I have been encountering an error (in the Winbugs window: Blackbox - index out of range ..... ). So to me it seems R2WinBUGS is not linking up well with the new R version since the same model runs well under r-2.9.2. So is my assumption correct or am I doing something wrong. The data I used and the code are in the help of R2WinBUGS (?bugs). Regards, -- Luwis Diya, Leuven

Hi! I am new to BUGS and running BUGS from R. I am trying to run a regression model from R, however I have this error message: Error in file(con, "wb") : cannot open the connection In addition: Warning messages: 1: In file.create(to[okay]) : cannot create file 'c:/Program Files/WinBUGS14//System/Rsrc/Registry_Rsave.odc', reason 'Permission denied' 2: In file(con,

Dear all I want to simulate a stochastic jump variance process where N is Bernoulli with intensity lambda0 + lambda1*Vt. lambda0 is constant and lambda1 can be interpreted as a regression coefficient on the current variance level Vt. J is a scaling factor How can I rewrite this avoiding the loop structure which is very time-consuming for long simulations? for (i in 1:N){ ... N <- rbinom(n=1,

Hi R Users, I have--what I think--is a very basic question. I'm trying to use R2WinBUGS and am having difficulty setting up the data properly for repeated measures data in a hierarchical model. First, I know that it is me...not the package, as I sucessfully used R2WinBUGS to run a basic regression. Second, I know that it's not my hierarchical model as I've run it successfully in

For some purpose, I need to transfer a NAs array to WinBugs through R2WinBugs, But I constantly got an error message:"'type' must be "real" for this format". Here is my data to transfer: x = matrix(data=NA,nrow=3,ncol=3) x = as.array(x) data <- list ("x") if I add a line to above setting, then I can pass R2WinBugs: x[1,1] = 0 If I manually input the

Hi all, I did some analysis using package R2WinBUGS to call WinBUGS. I set the iterations to 50000 (fairly a large number, I think), but after the program was done, the effective posterior samples contained only 7 draws. I don't know why. By the way, I checked posterior sample size by using bugsobj$n.sims. And, for my previous practice with WinBUGS/R2WinBUGS, no such strange thing happend.

I am generating 1000 replicate data sets in R, each data set is then analyzed with WinBUGS in batch mode using R2WinBUGS. Unfortunately, occasionally some data sets lead WinBUGS to open a "trap" window; and the simulations are interrupted as result of the message. Is there any ways to set R2WinBUGS so that it would ignore the trap message and proceed with the simulations? Cheers

Hi there - I have a question for any of you who use R2WinBUGS to call WinBUGS using the useWINE option in that package. I have Codeweaver's Crossover emulator installed on my Intel Mac with WinBUGS working fine if I start it from Crossover. But there is supposed to be a way to start it directly from a command line, which I could pass to the R2WinBUGS using the bugs(...,

Hi all, I am having problems using the R2Winbugs function When I perform an analysis directly in Winbugs I can specify that the first 'n' iterations are to be done using an 'adaptive' phase. After this phase the markov chain seems to mix a lot better. I don't seem to be able to specify R2winbugs to carry out this adaptive phase, I can just specify the burnin length but

I am trying to figure out if it is possible to run winBUGS from within R, using R2winBUGS, without having winBUGS spawn any windows (basically - 'true' batch - no GUI actions at all). The reason being I have a machine which I (and several others) ssh/telnet into, and would like to run winBUGS without having to mount a virtual desktop of any kind. I've looked through the r2winBUGS