similar to: help with time series

Hi, I make this model: > moths.m6 <- glm(nind~metros+especie+habitat+metros:habitat+habitat:especie,family=poisson) > anova(moths.m6,test="F") Analysis of Deviance Table Model: poisson, link: log Response: nind Terms added sequentially (first to last) Df Deviance Resid. Df Resid. Dev F Pr(>F) NULL 81 488.32

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

Hi all, I have a little doubt. Is correct to make a regression when the independent variate is the time? If not, why and how to analyse my data. I dont want to know about time series or survival analysis. For exemple: Y(growing) X(days) 1 1 3 2 5 3 6 4 7 5 Thanks Ronaldo -- Mind your own business, then you don't mind mine. -- | //|\\

Hi all, I try to make a residuals plot in a model, but the warnings occur. > anova(m2) Analysis of Deviance Table Model: Negative Binomial(5.6937), link: log Response: nind Terms added sequentially (first to last) Df Deviance Resid. Df Resid. Dev P(>|Chi|) NULL 81 270.242 metros 1 0.005 80 270.237 0.943 habitat

Hi, I make a CRAN mirror in Brazil at Federal University of Vi?osa. The address is: http://termix.ufv.br/CRAN Every problem please send report for me. Thanks Ronaldo -- Hug me now, you mad, impetuous fool!! Oh wait... I'm a computer, and you're a person. It would never work out. Never mind. -- | //|\\ [*****************************][*******************] || ( ? ? ) [Ronaldo

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 all, I have some little questions. 1) I make a function, with list's help, for make a map, its work OK, but it show an extra output of this: lapply(split(dataset,dataset$trajectory),lfun) Is possible to eliminate any text output from an function? 2) 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?

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 all, Firstly excuse-me by this off-topic question. I'm very confused to decide what is the better analysis to use. Whem the experiments is a tradicional design of manipulative experiments, it is easy to decide. But in natural observations is too complicated. I have 6 farms, each farm cultives some sugarcane's varieties (maybe or not the same varieties by farm). For each farm I make

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 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 have make an analise with presence and absence, y=(1 e 0). I have a area continuous data and a sp data with 25 levels. I have 300 points. When I make glm((presenca/peso)~area,weights=peso,family=binomial,maxit=1000) where presenca is 0 or 1. peso is the unit = 1. area is the continuous data. The analysis is OK. When I put the sp and interactions in analysis this warning appear.

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