similar to: glm.nb

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 All, This question was brought up some time ago but I never saw a reply so I'd like to bring it up again. When using ROracle package (version 0.5-5), I am unable to run any queries that are greater than 4000 characters in length. If I do, I get the following message: Error in oraPrepareStatement(con, statement, bind=NULL) : RS-DBI driver: (too long a statement -- it must has less than

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

I think you've badly misinterpreted the purpose of the R listserv with this request: https://www.stat.math.ethz.ch/mailman/listinfo/r-help says "The `main' R mailing list, for announcements about the development of R and the availability of new code, questions and answers about problems and solutions using R, enhancements and patches to the source code and documentation of R,

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, 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

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, I have a R mirror on www.termix.ufv.br/CRAN I try to keep it updated using rsync and cron. rsync -rzuav --progress --delete cran.r-project.org::CRAN /var/www/termix/mirrors/linux/CRAN/ & But all times the rsync abort with error: 95.99kB/s 0:00:02 548960 77% 110.50kB/s 0:00:01rsync: read error: Connection reset by peer (104) rsync error: error in rsync protocol data

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

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, I try to make a key alias to XEmacs+ESS function Complete File Name. But I dont know the name of this funcion. Anybody know how this funcion is call? (global-set-key "\C-tab" 'name of the function Complete File Name) Thanks for all Inte Ronaldo -- To beer or not to beer. -- ShakesBeer -- | // | \\ [*****************************][*******************] ||

Hi, I have this non-linear function: y=115.251 - 118.69 * exp(-0.123517*x) I try to discovery the x value for the assyntoptic value of y and the x value where the behavior of curve change, the inflection point. How to make this? Is poss?ble to make this on R? Thanks Ronaldo -- O Flamengo e o Bangu terminaram o jogo em 0 x 0. Quem fez o gol? A Volkswagen -- | // | \\

Hi, I try to make a vector in a for for loop, but it dont work. Look: > a <- 0;for(i in c(1:2)) { for(j in c(1:2)) { a <- i+j; print(a)}} [1] 2 [1] 3 [1] 3 [1] 4 I try to make this a vector, like this: [1] 2 3 3 4 > a <- 0;for(i in c(1:2)) { for(j in c(1:2)) { a[j] <- i+j}}; print(a) [1] 3 4 > a <- 0;for(i in c(1:2)) { for(j in c(1:2)) { a[i] <- i+j}}; print(a) [1]

Hi, I know that is not possible make a stepwise procedure using REML in R, I can use ML for this. For nested design it may be very dangerous due the difference in variance structure, mainly in a splitplot design. ML make significative variables that REML dont make. I read an article that is made a stepwise procedure using GENSTAT. from article: "Terms were dropped from a model in a

Ronaldo Your mathematical question is of interest to me too. But I am particularly interested in a translation into American English of the pun attached to you email. Could you try to explain O Flamengo e o Bangu terminaram o jogo em 0 x 0. Quem fez o gol? A Volkswagen My rough translation into English is: Two soccer teams ended play 0 to 0. Who got the goal? A Volkswagen. I need all the

Hi, I have problem with the restore function in some .RData using R 1.9.0 Look the error: [ronaldo at zeus RAnalise]$ R ... Error: object 'family' not found whilst loading namespace 'MASS' Fatal error: unable to restore saved data in .RData But if I load this .RData with the load() function all objects are recovered. [ronaldo at zeus RAnalise]$ R --no-restore-data ... > ls()

Hi all, I have problems in a dataframe variables types. Look: from a loop function: for(...){ ... dados.fin <- rbind(dados.fin, c(L=j, A=j^2, Nsp=nsps, N=length(amosfin$SP), AmT="am",NAm=nam, AMST=amst)) dados.fin <- rbind(dados.fin, c(L=j, A=j^2,

Hi all, I'm here again with my newbies questions :( I have a simple example: count of slugs in two fields. I need to make a barplot with mean and SE of mean. So I have: The mean: > tapply(slugs,field,mean) Nursery Rookery 1.275 2.275 The SE: > tapply(slugs,field,sd)/sqrt(tapply(slugs,field,length)) Nursery Rookery 0.3651264 0.3508004 If the data has been normally

Hi, I have an data like this x <- c(1:11) y <- seq(245,445,20) I make a plot plot(x,y) The first y value in y axis is 250 and the last y valeu is 450, the spacing between values are 50, so the y values showed are: 250, 300, 350, 400, 450 I need that values in y axis are: 240, 280, 320, 360, 400, 440, 480. I try plot(x,y,yaxp=c(240,480,6)) dont work I try

Hi, I try to make a binomial analysis using GLMM in a longitudinal data file. Is correct to use anova(model) to access the significance of the fixed terms? Thanks Ronaldo -- Todos somos iguais perante a lei, mas nao perante os encarregados de faze-las cumprir. -- S. Jerzy Lec -- |> // | \\ [***********************************] | ( ? ? ) [Ronaldo Reis J?nior ]