[R] Help on survival analysis

Hi, I need a help in a survival analysis using survreg function with weibull distribution from survival package. Look the data sample: ## Start of script dados - structure(list(TFD = c(20L, 34L, 1L, 2L, 3L, 3L, 50L, 26L, 1L, 50L, 21L, 3L, 13L, 11L, 22L, 50L, 50L, 1L, 50L, 9L, 50L, 1L,

[R] Help in cluster MPI and R

Hi, ImountedanMPIclusterwithdebian.I installedthepackagesnecessarytorunusingRinthiscluster,butunsuccessfully, thefunctionsofexampleswere runbutonlyin theserver processorin the cluster.I thinkitlacks an essentialconfigurationofR packagesfor him tocommunicatewith

[R] Very simple question

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

Re: [R] Very simple question

Thanks That is it really simple. Ronaldo Em 08-04-2011 07:48, Mario Valle escreveu: paste(teste, collapse=+) ciao! mario On 08-Apr-11 12:44, Ronaldo Reis Junior wrote: teste - c(A,B,C) -- 1ª lei - Suas férias começam após a defesa e entrega de sua dissertação

[R] Stress in metaMDS from vegam package

Hi, I have a newbie doubt. I try to use metaMDS in R, it is OK, but the stress that I get is 16.32, in other programs stress is showed like 0.173. I don't find a good explanation about this calculations on books, I think that R show stress em a percentage data (0-100) and other software in a

Re: [R] Problem with vegan package instalation on linux

Gavin, yes, I have GCC, Gfortran, [root@amandix ~]# gcc --version gcc (Debian 4.5.2-5) 4.5.2 [root@amandix ~]# gfortran --version GNU Fortran (Debian 4.5.2-5) 4.5.2 [root@amandix ~]# make --version GNU Make 3.81 I have another linux installed that I dont have any problem, the installations

Re: [R] Problem with vegan package instalation on linux

Hi, I try with --no-multiarch and it work now. I use the same R version in both systems: R version 2.12.2 (2011-02-25) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-pc-linux-gnu (32-bit) Now, how to make the --no-multiarch option default to R

[R] Using merge

Hi, I have two tables and I need to merge both. I use the merge command, but in this way the name must be exactly. How I can make to compare independently of upper or lower-case? Look: data1-data.frame(journal=c(Ecology,Environmental Entomology,Neotropical Biology And Conservation))

Re: [R] Using merge

Hi, ignore my e-mail, a just use tolower function. thanks and sorry Ronaldo Em 11-02-2011 11:24, Ronaldo Reis Junior escreveu: Hi, I have two tables and I need to merge both. I use the merge command, but in this way the name must be exactly. How I can make to compare independently

[R] sum with times

Hi, I have this vector: A - c(00:00:36,00:02:18) A [1] 00:00:36 00:02:18 I use as.difftime to convert this to time vector based B - as.difftime(A) B Time differences in secs [1] 36 138 attr(,tzone) [1] Now i try to make a sum sum(B) Time difference of 174 secs The sum is ok, but

[R] sum with times

Hi, Forget, the chron package work with this Thanks Ronaldo -- 8ª lei - Colete seus dados hoje como se você soubesse que seu equipamento vai quebrar amanhã. --Herman, I. P. 2007. Following the law. NATURE, Vol 445, p. 228. Prof. Ronaldo Reis Júnior | .''`.

[R] Color on trellis device

Hi, I try to make a xyplot withou colors. I try: pdf(file=test.pdf) trellis.device(color=F) Depth - equal.count(quakes$depth, number=8, overlap=.1) xyplot(lat ~ long | Depth, data = quakes) dev.off() the graphic is showed without colors in a X device but it dont work saving directly on pdf

[R] Different colour in each bar in lattice package

Hi, I try to plot bars with different colours in a barchart graphic. My idea is make that all X-Levels from trat var with different colour (grey scale). I search for a solution but dont find any. Any help? Thanks dados - structure(list(Medias = c(0.994169096209855, 0.99416342412449,

[R] Help with tapply selection

Hi, I try to create a data.frame from results of a tapply. I have this results from tapply: -- , , tapinom sp 1 1 2 3 4 5 6 7 8 9 10 1 NA NA NA NA NA NA NA NA NA NA 4 NA NA NA NA NA NA NA NA NA 4 5 NA NA NA NA NA NA NA NA NA NA 6 NA NA NA NA NA NA NA NA NA

[R] Help on rKward

Hi, I'm testing rKward and it's become a great GUI for R on linux, mainly for new linux users. I'm not a new linux user and I use Emacs for my own R's script. But I always try new GUIs or IDEs for to recommend to my students. The most dificult for new R and Linux users is: I have installed R

[R] Yet Testing rKward

Hi, I continue testing rKward. I dont know how to save the results from a script execution without use copy and paste or using a rkward output system. Now I try to understand how adapt my script do use the rKward output system. Example: I have this script: -- ## Carregar a

[R] change the panel name in xyplot

Hi, I try to change the panel name in a xyplot without success. Look this example from xyplot manual: xyplot(Murder ~ Population | state.region,data=states) The panel title are: Northeast, South, North Central, West, that are factor from state.region. I need do change some names and, for

[R] How to increase the iteration limit in lmer

Hi, I have this message: Warning message: In mer_finalize(ans, verbose) : iteration limit reached without convergence (9) How to increase this limit? I try control=list(maxit=1000) I try maxIter and it dont work. Any idea Thanks Ronaldo -- Opportunities are usually disguised as hard work,

[R] Add arrows to a filled.contour margin

Hi, I try to put an arrow to a filled.contour plot. I make this: filled.contour(volcano, color = terrain.colors,key.axes=F) arrows(0.95,0,0.95,1,lwd=2) This work, but the problem is that I find the coordinates at random try, for other graphic with other scale is other fight with coordinates.

[R] How to get the p-value from lmer on a longitudinal analysis

Hi, I have a modelo like this: Yvar - c(0, 0, 0, 0, 1, 0, 0, 0, 1, 2, 1, 1, 2, 3, 6, 6, 3, 3, 4) TIME - 4:22 ID - rep(PlotA,19) m - lmer(Yvar~TIME+(TIME|ID),family=poisson) anova(m) summary(m) How to get the p-value for this case? Thanks Ronaldo -- Just because you're paranoid doesn't mean

Re: [R] How to get the p-value from lmer on a longitudinal analysis

Em Sex 01 Ago 2008, Marc Schwartz escreveu: on 08/01/2008 10:11 AM Ronaldo Reis Junior wrote: Hi, I have a modelo like this: Yvar - c(0, 0, 0, 0, 1, 0, 0, 0, 1, 2, 1, 1, 2, 3, 6, 6, 3, 3, 4) TIME - 4:22 ID - rep(PlotA,19) m - lmer(Yvar~TIME+(TIME|ID),family=poisson) anova(m

Re: [R] Errors bar in barchart

Em Sex 02 Mai 2008, Deepayan Sarkar escreveu: On 5/2/08, Ronaldo Reis Junior [EMAIL PROTECTED] wrote: Hi, I user barplot2 to make a plot bar with errors bars. In old times I needed to use a sequence of segments commands to make this. Now I try to make the same but using lattice

[R] Errors bar in barchart

Hi, I user barplot2 to make a plot bar with errors bars. In old times I needed to use a sequence of segments commands to make this. Now I try to make the same but using lattice. Is possible to use barplot2 in barchart function? If not, what is the simplest way to put errors bar in barchart? I

Re: [R] Survival Regression with multiple records per subject

Em Qua 30 Abr 2008, Fabian Hefner escreveu: Here is an example:      itm      ID     exercise      time 1.401869    1        0             1 1.324390    1        0             2 1.324390    1        0             3 1.38    1        0             4 1.346761    1        0             5

Re: [R] R on X11 under Linux (newbie)

Em Seg 28 Abr 2008, Dumblauskas, Jerry escreveu: Thx I verified by typing in plot(1,1) at the R prompt -- sure enough an X window popped up on my desktop. Off to investigate Linux GUI front ends! Thx again Hi, at this moment the best IDE to use R on Linux in my opinion are in sequence:

[R] Problem with object inside a function

Hi all, I have a problem with an object inside a formula in a function that I make. I create an object named qvartemp. This object is OK, look the levels: print(levels(qvartemp)) [1] baixomedio alto Now I create a new object with a new formula: new.form -

Re: [R] Problem with object inside a function

Em Dom 20 Abr 2008, Ronaldo Reis Junior escreveu: Hi all, I have a problem with an object inside a formula in a function that I make. I create an object named qvartemp. This object is OK, look the levels: print(levels(qvartemp)) [1] baixomedio alto Now I create a new object with a new

[R] Script to start Rcmdr

Hi, anybody know any way to make a bash script to start Rcmdr directly widthout need to open R and execute library(Rcmdr)? I try to make this using R CMD BATCH somethink, but dont find the way. Thanks Ronaldo -- Prof. Ronaldo Reis Júnior | .''`. UNIMONTES/Depto. Biologia Geral/Lab. de

[R] update don't find an object

Hi, I have a situation here. I try this update: mmaa - update(mma,biomass~qvartemp) but I have this message: Error in eval(expr, envir, enclos) : object qvartemp not found but this object exist: [1] cont i levelsord mma qvar qvarmma [7] qvartemp test yvar

Re: [R] factors ordered by mean

Em Qua 20 Fev 2008, Henrique Dallazuanna escreveu: Try this: factor(pH,levels=names(sort(tapply(Riqueza,pH,mean Thanks it work Ronaldo -- Prof. Ronaldo Reis Júnior | .''`. UNIMONTES/Depto. Biologia Geral/Lab. de Biologia Computacional | : :' : Campus Universitário Prof. Darcy Ribeiro,

Re: [R] update don't find an object

Em Qui 21 Fev 2008, David Winsemius escreveu: Ronaldo Reis Junior [EMAIL PROTECTED] wrote in news:[EMAIL PROTECTED]: I try this update: mmaa - update(mma,biomass~qvartemp) That does not look like it has proper arguments to update.formula. As the update() help page suggests

[R] factors ordered by mean

Hi, How to order the levels os factor not by alphabetic order but by mean of Y. Somethink like this: I have this alphabetic order: levels(pH) [1] alto baixo medio the order by mean os yvar is: sort(tapply(Riqueza,pH,mean)) baixomedio alto 11.56667 20.0 26.8 How to

[R] Very simple question on plot

Hi, I try to make a plot like this: Y | | | o | oo | o o | o o |o o o | 0 10 20 30 40 50 60 70 80 90 100 (A) 100 90 80 70 60 50 40 30 20 10 0 (B)

[R] Package to make stepwise model selection using F or Chisq test

Hi, I looking for a method that use F or Chisq test instead of AIC in a stepwise modelo selection. I try the grasp package using the grasp.step.anova, but It dont work. library(grasp) Carregando pacotes exigidos: gam Carregando pacotes exigidos: splines Carregando pacotes exigidos: mda

[R] help in substitute function

Hi, I try to use a substitute function to generalise a equation. I have this: expression(1+2*pred+3*lat) I need to define the equation in function of x1 and x2 variables. Also I try to define who is x1 and x2. I try this: X1 - pred X2 - lat Now I need substitute pred and lat in equation by