similar to: Help on a simple function.

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible URL : <https://stat.ethz.ch/pipermail/r-help/attachments/20080725/bcafe250/attachment.pl>

Hello, I am using Hmisc summary.formula, latex and Sweave to produce tables for publication. Is it possible to change the formats for binary and continuous variables? I would prefer to show 35 (10%) and 1.5 (1.2-1.8) rather than 10% (35) and 1.2 / 1.5 / 1.8. Here is a simple example: sex <- factor(sample(c("m","f"), 500, rep=TRUE)) age <- rnorm(500, 50, 5) treatment

Hi, I have a very simple doubt. Look: > teste <- c("A","B","C") > teste2 <- paste(teste[1],teste[2],teste[3],sep="+") > teste2 [1] "A+B+C" > How to make it automatic, like I try to use paste(teste,sep="+") but the paste dont get the teste elements separately to join again in a unique element. Exist any

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: ...

Hola colegas: He estado usando R Studo para obtener unos reportes que espero me ayuden a comprender un ejemplo que estoy tomando de la red. Usando sink("Tabla.txt") , al pedir summay(n), obtengo: sr tl 1 : 1 1: 4 10 : 1 2:13 18 : 1 3:10 30 : 1 4: 4 47 : 1 7: 9 67 : 1 (Other):34 d Min. : 1.000 1st Qu.: 2.000 Median : 5.000

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,

Hello all, I'm trying to use confint from the MASS package to compute confidence intervals for an nls object. When I plot the results, however, they don't make sense - lines cross over the fitted model or just don't match the data. Code is : Thanks for help dat<-data.frame(a,b) with(dat, plot(a,b)) model<-(nls(b~(1/exp(a*x))*n, data=dat, start=list(x=.001,n=20),

Hi, I have problem with the restore function in some .RData using R 1.9.0 Look the error: [ronaldo at zeus RAnalise]$ R ... Error: object 'family' not found whilst loading namespace 'MASS' Fatal error: unable to restore saved data in .RData But if I load this .RData with the load() function all objects are recovered. [ronaldo at zeus RAnalise]$ R --no-restore-data ... > ls()

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, I make a upgrade to R 2.1.0 and in some analysis I give an error: anova(model,test="F") Analysis of Deviance Table Model: binomial, link: logit Response: landing/total Terms added sequentially (first to last) Df Deviance Resid. Df Resid. Dev F Pr(>F) NULL 16 105.079 trat 1 93.149 15 11.930 93.15 Warning message:

Is possible to integrate this diferential equation: dN/dt = Nr(1-(N/K)) in R using the integrate() function? Or any other diferential equation? If yes, how? If no, anybody know any software on linux that make this? Inte Ronaldo -- O problema de ter os dois p??s bem firmes no ch??o ?? que voc?? n??o vai conseguir tirar as cal??as. --An??nimo -- |> // | \\

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