similar to: is possible to use update with lmer models

Hi, I have these vectors: > WEEK <- rep(c(1:52),2) > YEAR <- rep(c(2000,1999),c(52,52)) How to make a vector of Date with weeks in years? I try as.date from survival package, but it dont work with weeks, just only with days, months etcs. Thanks Ronaldo -- "Realmente minha cidade e muito facultativa" --Elivelton, ao rep?rter da Jovem Pan que falava das muitas

Hi, I have the follow function: function() { ## Init of function ... for(i in test) { ... while(j <= test2) { ... } } } The problem is that sometimes, naturally, the while is not possible to be resolved, and so the program abort. In this case I need that program return to the init of function and run again. How I can make this? Abort the while, abort the for

Hi, I have some data that have this behaviour: | |******* | * | * | * | * |---------------- What is the best and simpler way to fit this in R? Thanks Ronaldo -- Ela pilotava um Continenal 2001 com igni????o autom??tica Magiclic... -- |> // | \\ [***********************************] | ( ?? ?? ) [Ronaldo Reis J??nior ] |>

Hi, I have this: lme(y~x1+x2,random=~1|x1/x2) How to make this random effect using lmer? I try this: lmer(y~x1+x2+(1|x1/x2) But it dont work. Any idea? Thanks Ronaldo -- System halted! -- |> // | \\ [***********************************] | ( ?? ?? ) [Ronaldo Reis J??nior ] |> V [UFV/DBA-Entomologia ] | / \ [36570-000 Vi??osa -

Hi, look this: > summary(fam??lia) Anacardiaceae Annonaceae Bombacaceae Cecropiaceae 2 4 1 3 Chrysobalanaceae Clusiaceae Euphorbiacaea Fabacea 1 1 4 3 Fabaceae Flacourtiaceae Humiriaceae indeterminada

Hi, I try to run the example of Crawley's Book on the page 661, but it fail, look > repmeasures <- read.table("../Packages/Crawley/data/repmeasures.txt",header=T) > attach(repmeasures) > rep <- as.factor(rep) > library(nlme) > model <- lme(height~seed,random=~time|rep/seed) Erro em lme.formula(height ~ seed, random = ~time | rep/seed) : iteration limit

Hi, I try to make an analyses to discover what is the time that an area begin to have spacial autocorrelation. And after, what is the number of individuals responsible for this autocorrelation. The main idea is to discover if exist a contamination of a quadrat from others quadrats and how is the population needed to make this contamination. This is very common to use automata to simulate

Hi, I make a upgrade to R 2.1.0 and in some analysis I give an error: anova(model,test="F") Analysis of Deviance Table Model: binomial, link: logit Response: landing/total Terms added sequentially (first to last) Df Deviance Resid. Df Resid. Dev F Pr(>F) NULL 16 105.079 trat 1 93.149 15 11.930 93.15 Warning message:

Hi, Anybody can suggest a good Ecological Modelling's Book. Thanks Ronaldo -- O c?u sueco. -- pal?ndromo -- |> // | \\ [***********************************] | ( ? ? ) [Ronaldo Reis J?nior ] |> V [UFV/DBA-Entomologia ] | / \ [36571-000 Vi?osa - MG ] |> /(.''`.)\ [Fone: 31-3899-2532

Its just a simple question I guess: I have a vector with missing information like x<-c(0,1,31,131,NA,133,NA,310,NA,112,3,1,2,93) How can I make a vector like this no missing in it. I used the x[x<0] commabd and tried some more, with no success. thx Korbinian

Hi, anybody have the dataset used in Gamma Errors chapter of the Crawley's books (An Introduction to Data Analysis using S-Plus). specifically the functionalresponse and the Density datasets. Thanks Ronaldo -- For every problem there is one solution which is simple, neat, and wrong. -- H. L. Mencken -- |> // | \\ [***********************************] | ( ? ? ) [Ronaldo Reis

Hi, Please, change the brazilian mirror http://www.termix.ufv.br/CRAN/ to http://www.insecta.ufv.br/CRAN/ in R homepage. Thanks ROnaldo -- Se dois homens no mesmo trabalho concordam o tempo todo, um deles ?? demais. Se discordam sem parar, ent??o os dois s??o dispens??veis --Darryl F. Zanuck -- |> // | \\ [***********************************] | ( ?? ?? ) [Ronaldo Reis J??nior

Is possible to integrate this diferential equation: dN/dt = Nr(1-(N/K)) in R using the integrate() function? Or any other diferential equation? If yes, how? If no, anybody know any software on linux that make this? Inte Ronaldo -- O problema de ter os dois p??s bem firmes no ch??o ?? que voc?? n??o vai conseguir tirar as cal??as. --An??nimo -- |> // | \\

Hi, exist in R or any other program for linux that estimate the best fit for data using severals functions? Somethink like tablecurve Thanks Ronaldo -- Um menino-prodigio e uma crianca cujos pais tem muita imaginacao. -- Jean Costeau -- |> // | \\ [***********************************] | ( ?? ?? ) [Ronaldo Reis J??nior ] |> V [UFV/DBA-Entomologia

Hi, I try to put a legend in a xyplot graphic. xyplot(y~x|g,ylim=c(0,80),xlim=c(10,40),as.table=T,layout=c(2,3), ylab="N??mero de machos capturados",xlab=expression(paste("Temperatura (",degree,"C)")), key=list(corner=c(0,0),x=0, y=0, text=list(legenda),lines=list(col=cor, lwd = espessura, lty=linha),columns=7,between=0.5,betwen.columns=0.5,cex=0.8)) The problem

Hi, At this moment, I make my R output using comments and latex tags in R. Something like this: ### make a plot plot(y~x) dev.copy2eps(file="plot.eps") # @\includegraphics[width=\linewidth]{plot}@ ... After I include the file saved with this code in a latex document using the listings package. It work. I reading about Sweave and it make a good output. But all example is made with

Hi, I make an analysis and depending of the order of the variables, the significance change, look. m1 <- glm((infec/ntot)~idade+sexo+peso,family=binomial,weights=ntot) > anova(m1,test="F") Analysis of Deviance Table Model: binomial, link: logit Response: (infec/ntot) Terms added sequentially (first to last) Df Deviance Resid. Df Resid. Dev F Pr(>F)

Hi, exist in R any glmm function that have any tools for test for overall goodness of fit? Thanks Ronaldo -- O papel da impressora ? sempre mais forte na parte picotada. -- | // | \\ [***********************************] |> ( ? ? ) [Ronaldo Reis J?nior ] | V [UFV/DBA-Entomologia ] |> / \ [36571-000 Vi?osa - MG ] |

Hi, I make an analysis sequence on R. In some cases the function glm.nb fail to ajust the model. Its Ok. The problem is that this error stop the program. I need treat this error and not stop the program. Something like this: ... model <- glm.nb(y~x,maxit=1000) if(glm.nb fail) { teste[i] <- 0 } else { teste[i] <- anova(modelo)$"P(>|Chi|)"[2] } ... I try this: ...

Hi, I'm try to use the SJava package. The install is OK. In R I have this error message: -------------------------- > library(SJava) Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to load shared library "/opt/lib/R/site-library/SJava/libs/SJava.so": libRSNativeJava.so: cannot open shared object file: No such file or directory Error in