similar to: plot survival points

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 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 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''m trying to understand the lme output and procedure. I''m using the Crawley''s book. I''m try to analyse the rats example take from Sokal and Rohlf (1995). I make a nested analysis using aov following the book. > summary(rats) Glycogen Treatment Rat Liver Min. :125.0 Min. :1 Min. :1.0 Min. :1 1st Qu.:135.8

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