similar to: WinBUGS on survival, simple but confusing question

You can make this much more readable with apply functions. result <- apply( all_combine1, 1, function(x){ p.value <- sapply( seq_len(nSims), function(sim){ gamma1 <- rgamma(x["m"], x["sp(skewness1.5)"], x["scp1"]) gamma2 <- rgamma(x["n"], x["scp1"], 1) gamma1 <- gamma1 -

Always keep the mailing list in cc. The code runs for each row in the data. However I get the feeling that there is a mismatch between what you think that is in the data and the actual data. ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance Kliniekstraat 25 1070

Hi Jeem, First, please send questions like this to the help list, not me. I assume that you are in a similar position to sjtan who has been sending almost exactly the same questions. The problem is not in the loops (which look rather familiar to me) but in your initial assignments at the top. For instance: scale parameter=(1,1.5,2,2.5,3) produces an error which has nothing to do with the

Hi Si Jie, Again, please send questions to the list, not me. Okay, I may have worked out what you are doing. The program runs and produces what I would expect in the rightmost columns of the result "g". You are storing the number of each test for which the p value is less than 0.05. It looks to me as though the objects storing the results should be vectors as you are only storing 100 p

Hi, i am sorry, the output should be values between 0 and 0.1 and not supposed to be 1.00, it is because they are type 1 error rate. And now i get output 1.00 for several samples,rhis is no correct. The loop do not run for every row. i do not know where is my mistake. As i use the same concept on normal distribution setup, i get the result. Sent from my phone On Thierry Onkelinx

i have combined all the variables in a matrix, and i wish to conduct a simulation row by row. But i found out the code only works for the every first row after a cycle of nine samples. But after check out the code, i don know where is my mistake... can anyone pls help .... #For gamma disribution with equal skewness 1.5 #to evaluate the same R function on many different sets of data

Dear anonymous, The big mistake in the output might be obvious to you but not to others. Please make clear what the correct output should be or at least what is wrong with the current output. And please DO read the posting guide which asks you not to post in HTML. ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest team Biometrie &

Hi, I have a bivariate longitudinal dataset. As an example say, i have the data frame with column names var1 var2 Unit time trt (trt represents the treatment) Now suppose I want to fit a joint model of the form for the *i* th unit var1jk = alpha1 + beta1*timejk + gamma1* trtjk + delta1* timejk:trtjk + error1jk var2 = alpha2 + beta2*timejk + gamma2* trtjk + delta2* timejk:trtjk +

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

Workshop on Bayesian methods and WinBUGS *************************************** A two-day workshop on Bayesian methods is being held on Friday 3 - Saturday 4 February 2012 at the University of Sydney. This course is suitable for graduate students, academics, researchers and professionals who wish to introduce themselves in the application of Bayesian methods and the use of WinBUGS software.

Hi, I'm having trouble running a script within winbugs when calling winbugs from R. Here are a few lines from the winbugs log: display(log) check(C:\DOCUME~1\michael\LOCALS~1\Temp\RtmpjKnHMu\mm_TO_bugs9bc21e.txt) cannot open C:\DOCUME~1\michael\LOCALS~1\Temp\RtmpjKnHMu\mm_TO_bugs9bc21e.txt data(C:/DOCUME~1/michael/LOCALS~1/Temp/RtmpjKnHMu/data.txt) command #Bugs:data cannot be executed (is

Hello everyone, I need some help regarding calling WinBUGS from R. I have a model for WinBUGS and an R code which calls WinBUGS. On running the code I am being shown the 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

Dear list members, Apologies for cross-posting. Please, find below the information of two statistical courses with R: 1) Statistical Analysis with R 2) Bayesian Data Analysis with R and WinBUGS If you have any question don't hesitate to contact me. Best regards, Pablo ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ *Two days course in: Statistical Analysis with R *Where:

Dear list members, Best greetings and apologies for cross-posting. There are available places for the course: "Bayesian Data Analysis with R and WinBUGS", please find the description of the course below. If you have any question don't hesitate to contact me. Best regards, Pablo ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ *Three days course in: Bayesian Data

Dear all, I want to run two winBugs process through R (by R2WinBUGS package). The first WinBUGS process can successfully perform under R. The error message appears while I put the code to run for the second WinBUGS. The error message: Error in file(file, ifelse(append, "a", "w")) : cannot open the connection In addition: Warning message: In file(file, ifelse(append,

http://biomet.oxfordjournals.org/cgi/reprint/86/3/677 biometrika1999 http://biomet.oxfordjournals.org/cgi/reprint/94/4/1006 biometrika2000 Hi All: I just want to try some luck. I am currenly working on my project,one part of my project is to reanalysis the kenward cattle data by using the method in Mohsen's paper,but I found I really can get the same or close output as he did,so,any

Hi, I am very new to R and Winbugs. I have been provided with a code that runs Winbugs from R. The code runs 60 different models and I would like to either save the Winbugs log file for each of these or reproduce the WinBugs plots which show convergence over the iterations.I can do this by using Debug=true and saving manually each time but obviously this is very time consuming. Many thanks Ruth

Hi, I am trying to use WinBUGS using R, but i am having some problems using bugs function: out <- bugs(data = win.data, inits = inits, parameters.to.save = params, model.file = "model.txt", n.thin = nt, n.chains = nc, n.burnin = nb, n.iter = ni, debug = TRUE, DIC = TRUE, working.directory = getwd()) ls() WinBUGS file is in working directory=D:/Line/documents but R looks for

Apologies for cross-posting. Places are available in the following course: "Introduction to Bayesian Data Analysis using WinBUGS and R" Best Regards, Pablo ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ Course: Introduction to Bayesian Data Analysis using WinBUGS and R Where: Linux Hotel, Essen-Horst, Germany

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