search for: mckann

I seem to come to this problem alot, and I can find my way out of it with a loop, but I wish, and wonder if there is a better way. Here's an example (lmer1-5 are a series of lmer objects): bs=data.frame(bic=BIC(lmer1,lmer2,lmer3,lmer4,lmer5)$BIC) rownames(bs)=c('lmer1','lmer2','lmer3','lmer4','lmer5') best=rownames(bs)[bs==min(bs)] > best [1]

Hello all. Why doesn't this work? d=data.frame(y=rpois(10,1),x=rnorm(10),z=rnorm(10),grp=rep(c('a','b'),each=5)) library(lme4) model=lmer(y~x+z+(1|grp),family=poisson,data=d) update(model,~.-z)###works, removes z var='z' update(model,~.-get(var))##doesn't remove z update(model,~. -get(var,pos=d))###doesn't remove z I am trying to remove z from the model in

Hello all, I am trying to write a program in R in which I call a function multiple times within a loop. The problem is that sometimes the function breaks down while calling another function, and produces an error message that breaks my loop and the program stops. I would like to keep the loop running when this function breaks down, and just move on to the next iteration in the loop. Is there

Hello all, I don't understand why this won't work. I have entered: WriteXLS(alldata,'test.xls') and I get this error message: Error in get(x, envir = envir) : variable names are limited to 256 bytes. My variable names are not very long, and are accepted by write.csv. alldata is a list containing 4 dataframes, with each dataframe having the the same variable names, which are:

...Winsemius) 66. Re: Encoding problem - I fails to read Hebrew text from online (Tal Galili) 67. spatial clusters (dorina.lazar) 68. Re: Time Series Row Label (Gabor Grothendieck) 69. WriteXLS error:Error in get(x, envir = envir) : variable names are limited to 256 bytes (Patrick McKann) 70. How to print colorful R output?? (casperyc) 71. Re: WriteXLS error:Error in get(x, envir = envir) : variable names are limited to 256 bytes (David Winsemius) 72. Re: WriteXLS error:Error in get(x, envir = envir) : variable names are limited to 256 bytes (David Winsemius) 73...

I tried calculating pD as 1/2 variance of the deviance, but I got hugely inflated numbers. Does anybody know how to calculate pD from the coda files output from R2WinBUGS? By the way, 'set DIC' is greyed out for some reason within WinBUGS, so I can't monitor DIC. [[alternative HTML version deleted]]

Hello. I have made many pngs in the past, but started having trouble with them yesterday. I just downloaded R 2.13.0 this morning, and was able to make some pngs earlier this morning. Now I am trying to remake them, but get the same error message I got yesterday: > u='in' > w=10.5 > h=8 > png('winbugs changepoint graphs ind hibs gt 400

Hi all. This is more of a stats question, I suppose. Let's say I have two separate simple regressions of weight on year from two different datasets. I want to combine the regressions so that I can come up with a single equation for the total weight regressed on year. In reality, there is missing data, so I can't just sum the data across datasets and come up with a regression on the