[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 me

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

[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. An

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

[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 o

[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] 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 NA

[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 d

[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

[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 | .''`. UNIMONTES/DBG/L

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

[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] 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")) d

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 >

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 a

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 in

[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 p

[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 elemen

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")

[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 thecluster.Anyoneknowanymorebasi

[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 <- as.formula(gsub(qvar,"qvarte

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] "baixomedi

Re: [R] Problem with object inside a function

Em Dom 20 Abr 2008, Duncan Murdoch escreveu: > On 20/04/2008 9:40 AM, Ronaldo Reis Junior wrote: > > 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

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 sequ

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        

[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] 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. > &

[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) o

[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 exam

[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

[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" "t

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 Ribei

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 u

[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 Biolog

[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 equatio

[R] Make playwith a default graphic device

Hi, I looking for a graphic device on Linux with more functionalities than the X11 (). I find the plotAndPlayGTK package and the playwith function. It is a good device. It is possible to make it default in my system. Now I need: > library(plotAndPlayGTK ) Loading required package: RGtk2 Loading

Re: [R] Make playwith a default graphic device

Em Dom 14 Out 2007, Prof Brian Ripley escreveu: > You seem to be under the misapprehension that playwith() is a graphics > device: it is not. > > To set a graphics device as the default, use options(device=). I am not > sure what 'more functionalities' you are looking for: the playwith() uses > th

[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 Carreg