similar to: warning in binomial analysis

Hi all, I try to understand the output of the summary(model). I make a simple model. > summary(m1anova) Call: glm(formula = peso ~ gen) Deviance Residuals: Min 1Q Median 3Q Max -4.3114 -2.3788 -0.9167 2.1581 5.7856 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 4.8660 0.9808 4.961 0.000101 *** gen2g 3.8504

Hi all, In my continuous transition of GLIM to R I try to make a glm with binomial errors. The data file have 3 vectors: h -> the factor that is ajusted (have 3 levels) d -> number of animais alive (the response) n -> total number of animals To test proportion of alive, make d/n. In GLIM: $yvar d$ $error binomial n$ $fit +h$ scale deviance = 25.730 (change = -9.138) at cycle 4

Hi all, I have experimented a segmentation fault error using RMySQL_0.4-6 library in R 1.5.1. Look the sequence: ------------------------------------------- > library(RMySQL) > m <- dbManager("MySQL") > m <MySQLManager:(1137)> > describe(m) <MySQLManager:(1137)> Driver name: MySQL Max connections: 10 Conn. processed: 0 Default records per

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 all, I have a little problem. I make an weibull survival analysis using the survival package. It,s OK, them I have the functions. I plot this funcions with curve(). I want to make a plot with the real survival points (proportion of alive x time) and them add the curves to points. I have the time to dead, the censor data and my trataments. To analysis the model is: model1 <-

Hi, I have a new (for me) situation. I have a function for one linear model: y = 9.7909-0.035*x I have the variance of system var = 13.01403 Is possible to create a random vector points of Y (simulating the real data) over this curve using the function and the variance? I need this to a discipline. I have this: | \ | \ -> y = 9.7909-0.035*x and var = 13.014 | \ | \ | \ | \

Hi, I have a little doubt about contrib packages. I make some little functions in a R package format, this functions are for my personal use (envelopes for glm models and a plot bar routine). Some friends ask me about the liberation of this functions in a web site, R contrib etc. what is the correct procedure for send a package to R contrib? The functions is not create for me, is adaptaded to R

Hi, I have a pentium III 600 with 128 ram running Linux. I have a dataset with this variables in 12000 lines. coleta = response (count data) data = continuous data latitude = continuous data longitude = continuous data variedade = categorical with 14 levels solo = categorical with 22 levels idade = continuous data app = categorical with 4 levels I want to test several models beginning by the

Hi how I must make to officialize my mirror in the main R website? My CRAN's mirror is sited in Federal University of Vi?osa in Minas Gerais State - Brazil. The address is: http://www.termix.ufv.br/CRAN It is diary updated Bie Ronaldo -- Q: How does a hacker fix a function which doesn't work for all of the elements in its domain? A: He changes the domain. -- | //|\\

Hi, I try to put an equation on a object to use in curve for example, but it don't work, it possible to make a object of an equation? ex. fx <- a + b*x for(a in 0){for(b in 1){curve(fx,...)}} Thanks you Ronaldo ps. Exist in R any functions to estimate severals curves parameters (non interactive) for a dataset? Something like table-curve 2d. -- If Karl, instead of writing a lot about

Hi, I make a graphic using xyplot, it is very good. xyplot((ocorrencia/isca)~frag|especie) But I need to plot a curve for each especie's plot. I try to use panel.abline and others panel., but I dont succeed. I need to plot the function exp(a+b*x+c*x^2)/1+exp(a+b*x+c*x^2), a b and c are different for each especie. Normally I use the curve function to add a curve in a points

Hi, I make a model and have compared with null model. anova(mnull,m1) Model theta Resid. df 2xlog-lik test df LR stat. Pr(chi) 1 ... 0.39 161 -577.9129 NA NA NA 2 ... 1.30 150 -475.6839 1 vs 2 11 102.229 1.11e10-16 anova(m1) Df Deviance Resid. Df Resid. Dev P(>|chi|) NULL 161 282.139 9... ... ...

Hi all, I have this, Weibull distribution n= 1000 > anova(pk.m0) Df Deviance Resid. Df -2*LL P(>|Chi|) NULL NA NA 999 1451.31 NA > anova(pk.m0,pk.m1) Terms Resid. Df -2*LL Test Df Deviance P(>|Chi|) 1 1 998 1451.310 NA NA NA 2 seed 996 1229.714 2 221.59557 7.605594e-49 With

Hi all, I try to test a linear slope using offset. I have: > m2 <- glm(Y~X*V) > summary(m2) Call: glm(formula = Y ~ X * V) Deviance Residuals: Min 1Q Median 3Q Max -2.01688 -0.56028 0.05224 0.53213 3.60216 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 1.3673 0.8476 1.613 0.119788 X

Hi, I have the graph below: | 422 | 1 321 | 1 | 1 |32 1 |23 |_______________________ The numbers are points in coordinates (simplified graph). I need make a cluster analysis of this for identify groups and make a graph with this groups under circle. In R exist several packages for cluster analysis, but, which the package for this analysis and

Hi all, I use a function update.packages to update all packages in my R instalation. But I need compare one by one packages to discovery new packages. My question is: Exist in R any function that compare all installed packages in my R and all packages available in CRAN and return a list of the new packages (ie packages that I dont have im my installation) and automatic download and install

Hi, I'm a new R's user. I am using "faq" and "An Introdution to R" to learn. My main problem for start the learning is: how to save my result in a output file? I has used the sink CMD or $ R > [output file], but the commands go to file and can't to be read in Linux Terminal. The workspace saved at the end of section save only commands and dont save results.

Hi, I have R installed in a system, but I need to install some packages as normal user on linux. Is possible to R use others library to get packages or one by one? Ex. I have some main package installed in /usr/lib/R/library And now I need to install my own packages in ~/R/library And I need to use the packages in main library and in my own library. Its possible? Thanks Ronaldo --

Hi, I make a glm analysis with binomial errors and its OK. Now I need plot the ajusted logist line in my points. How do I make it? | * * * * -> I need plot an ajusted | * * line in this points. | * * | * | * * | * * |* * * * |_________________________ Exist an abline like function for ajust others formulas? Exist in R

Hi, Is possible to anova function on R calculate the significance of complete glm model? I explain better. I have a glm model: > model1 <- glm(nsp~caules+mo,family=poisson) I make an anova of this: > anova(model1,test="F") Analysis of Deviance Table Model: poisson, link: log Response: nsp Terms added sequentially (first to last) Df Deviance Resid. Df Resid. Dev