similar to: R and default library

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 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 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 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 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, 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 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 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, I have a simple question. I have a variable with 3 levels. Variable A A A A A A B B B B B B C C C C C C I need to make a new vector renaming de levels B and C for D. NewVariable A A A A A A D D D D D D D D D D D D I dont find how to make it in manuals or in this list. Can anyone help-me??? Thank you Ronaldo -- Agnes' Law: Almost everything in life is easier to get into than out of.

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