Re: [R] R2WinBUGS stopping execution

2009-01-13 Thread Gregor Gorjanc
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 R2WinBUGS just passes

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

2009-01-01 Thread Gregor Gorjanc
Anny Huang 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 that you ha

Re: [R] Automatic placement of Legends

2008-07-08 Thread Gregor Gorjanc
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. Gregor _

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

2008-05-12 Thread Gregor Gorjanc
Bert Gunter 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 http://en.wikipedia.org/wiki/Tobit_model I think y

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 post

Re: [R] R2WinBUGS sending variables as factors

2008-01-05 Thread Gregor Gorjanc
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 Syntactically > co

Re: [R] Sweave and Scientific Workplace

2007-12-23 Thread Gregor Gorjanc
Dietrich Trenkler 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 > experiences >

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

2007-12-05 Thread Gregor Gorjanc
Henrique Dallazuanna gmail.com> writes: > x[is.na(x)] <- 0 > > On 05/12/2007, Amit Patel 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) x[1, 1] <- NA NA

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

2007-11-23 Thread Gregor Gorjanc
Thompson, David (MNR 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), 'dbh.n'] <- 0

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

2007-11-10 Thread Gregor Gorjanc
Rick Bilonick 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. Neith

Re: [R] Meta-analysis mixed model

2007-10-27 Thread Gregor Gorjanc
Diane Srivastava 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 explanato