similar to: integrate function

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 these vectors: > WEEK <- rep(c(1:52),2) > YEAR <- rep(c(2000,1999),c(52,52)) How to make a vector of Date with weeks in years? I try as.date from survival package, but it dont work with weeks, just only with days, months etcs. Thanks Ronaldo -- "Realmente minha cidade e muito facultativa" --Elivelton, ao rep?rter da Jovem Pan que falava das muitas

Hi, Please, change the brazilian mirror http://www.termix.ufv.br/CRAN/ to http://www.insecta.ufv.br/CRAN/ in R homepage. Thanks ROnaldo -- Se dois homens no mesmo trabalho concordam o tempo todo, um deles ?? demais. Se discordam sem parar, ent??o os dois s??o dispens??veis --Darryl F. Zanuck -- |> // | \\ [***********************************] | ( ?? ?? ) [Ronaldo Reis J??nior

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 have the follow function: function() { ## Init of function ... for(i in test) { ... while(j <= test2) { ... } } } The problem is that sometimes, naturally, the while is not possible to be resolved, and so the program abort. In this case I need that program return to the init of function and run again. How I can make this? Abort the while, abort the for

Hi, I have some data that have this behaviour: | |******* | * | * | * | * |---------------- What is the best and simpler way to fit this in R? Thanks Ronaldo -- Ela pilotava um Continenal 2001 com igni????o autom??tica Magiclic... -- |> // | \\ [***********************************] | ( ?? ?? ) [Ronaldo Reis J??nior ] |>

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

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''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 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 have this: lme(y~x1+x2,random=~1|x1/x2) How to make this random effect using lmer? I try this: lmer(y~x1+x2+(1|x1/x2) But it dont work. Any idea? Thanks Ronaldo -- System halted! -- |> // | \\ [***********************************] | ( ?? ?? ) [Ronaldo Reis J??nior ] |> V [UFV/DBA-Entomologia ] | / \ [36570-000 Vi??osa -

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

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, look this: > summary(fam??lia) Anacardiaceae Annonaceae Bombacaceae Cecropiaceae 2 4 1 3 Chrysobalanaceae Clusiaceae Euphorbiacaea Fabacea 1 1 4 3 Fabaceae Flacourtiaceae Humiriaceae indeterminada

Hi, I try to run the example of Crawley's Book on the page 661, but it fail, look > repmeasures <- read.table("../Packages/Crawley/data/repmeasures.txt",header=T) > attach(repmeasures) > rep <- as.factor(rep) > library(nlme) > model <- lme(height~seed,random=~time|rep/seed) Erro em lme.formula(height ~ seed, random = ~time | rep/seed) : iteration limit

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

Its just a simple question I guess: I have a vector with missing information like x<-c(0,1,31,131,NA,133,NA,310,NA,112,3,1,2,93) How can I make a vector like this no missing in it. I used the x[x<0] commabd and tried some more, with no success. thx Korbinian

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, 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 -- | // | \\ [*****************************][*******************] ||