Re: [R] R2WinBUGS stopping execution

2009-01-13 Thread Gregor Gorjanc
Richard.Cotton at hsl.gov.uk writes: I'm running OpenBUGS model via the R2WinBUGS package interface, under Windows. Is it possible to terminate running models, short of using the Windows Task Manager to forcibly exit the program? If you use OpenBUGS, then I guess you can not since

Re: [R] WinBUGS posterior samples (via R2WinBUGS)?

2009-01-01 Thread Gregor Gorjanc
Anny Huang annylhuang at gmail.com writes: I did some analysis using package R2WinBUGS to call WinBUGS. I set the iterations to 5 (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. This indicates

Re: [R] Automatic placement of Legends

2008-07-08 Thread Gregor Gorjanc
tolga.i.uzuner at jpmorgan.com writes: I am looking for a way to get legends placed automagically in an empty spot on a graph. Additional complication comes through my useage of multiple graphs on the same plot through mfrow. Take a look in Hmisc package. There is function for this task.

Re: [R] Left censored responses in mixed effects models

2008-05-12 Thread Gregor Gorjanc
Bert Gunter gunter.berton at gene.com writes: Dear R Fellow-Travellers: What is your recommended way of dealing with a left-censored response (non-detects) in (linear Gaussian) mixed effects models? Your description of the data calls for a tobit model

Re: [R] rbugs on linux and wine

2008-04-21 Thread Gregor Gorjanc
Hi Alexander! You are mixing WinBUGS and OpenBUGS. R package Rbugs works with OpenBUGS, but the later does not work with Rbugs under Linux! Gregor __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

Re: [R] R2WinBUGS sending variables as factors

2008-01-05 Thread Gregor Gorjanc
James.Dell at csiro.au writes: ... 1) I can't seem to send variables classed as factors (Month), is there a way do this? You can not use factors per se in BUGS. You have to convert them to numeric (integer) variables before. 2) Checking the Log in WinBUGS I can see that the model is

Re: [R] Sweave and Scientific Workplace

2007-12-23 Thread Gregor Gorjanc
Dietrich Trenkler Dietrich.Trenkler at uni-osnabrueck.de writes: Dear HelpeRs, a colleague of mine uses Scientific Workplace to write his LaTeX documents. I made his mouth water mentioning the advantages of using Sweave. Not using SW myself I wonder if anyone out there has gathered some

Re: [R] Dealing with NA's in a data matrix

2007-12-05 Thread Gregor Gorjanc
Henrique Dallazuanna wwwhsd at gmail.com writes: x[is.na(x)] - 0 On 05/12/2007, Amit Patel amitpatel_ak at yahoo.co.uk wrote: Hi I have a matrix with NA value that I would like to convert these to a value of 0. any suggestions also library(gdata) x - matrix(rnorm(16), nrow=4, ncol=4)

Re: [R] Summary: Process multiple columns of data.frame

2007-11-23 Thread Gregor Gorjanc
Thompson, David (MNR David.John.Thompson at ontario.ca writes: Thank you Jim Holtman and Mark Leeds for your help. Original question: How do I do the following more concisely? Bout[is.na(Bout$bd.n), 'bd.n'] - 0 Bout[is.na(Bout$ht.n), 'ht.n'] - 0 Bout[is.na(Bout$dbh.n),

Re: [R] Confidence Intervals for Random Effect BLUP's

2007-11-10 Thread Gregor Gorjanc
Rick Bilonick rab at nauticom.net writes: ... I think prediction interval is what is usually used. Regardless, I'm not sure how predict.lm will be of much help because I asked specifically about BLUP's for random effects and the last time I checked lm did not handle mixed effects models.

Re: [R] Meta-analysis mixed model

2007-10-27 Thread Gregor Gorjanc
Diane Srivastava srivast at zoology.ubc.ca writes: I have a meta-analysis dataset which I would like to analyze as a mixed model, where the y-variable is a measure of effect size, the random effect is the study from which the effect size was extracted, and the fixed effect is a categorical