similar to: Add arrows to a filled.contour margin

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 they AREN'T after you. -- > Prof. Ronaldo Reis J?nior |

Hi, this is a simple question I have this data.frame: > test <- data.frame(var1=c(1,1,1,1,1,1),var2=c("a","a","b","c","d","e"),var3=c("a1","a1","b1","a1","c1","d1")) > test var1 var2 var3 1 1 a a1 2 1 a a1 3 1 b b1 4 1 c a1 5 1

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 example, put some of these in italic. I dont find how change this. I looking for this in Deepayan

Hi, I have a table like that: > datatest var1 var2 var3 1 1 1 1 2 3 1 2 3 8 1 3 4 6 1 4 5 10 1 5 6 2 2 1 7 4 2 2 8 6 2 3 9 8 2 4 10 10 2 5 I need to create another table based on that with the rules: take a random sample by var2==1 (2 sample rows for example): var1 var2 var3 1 1

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 x1 and x2 defined on X1 and X2 above, I try this way: >

Hi, I try to make a plot like this: Y | | | o | o o | 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) or B 100 90 80 70 60 50 40 30 20 10 0 Y |-------------------------------- | | o

Hi, I use glht to make multcomp, using Tukey, from a glm model. It is possible to get a more simplified output of result? Somethink like ordering by letters. Thanks Ronaldo -- Human kind cannot bear very much reality. -- T. S. Eliot, "Four Quartets: Burnt Norton" -- > Prof. Ronaldo Reis J?nior | .''`. UNIMONTES/Depto. Biologia Geral/Lab. de Ecologia | : :' :

Hi, I dont know if this a coding rule or a program bug. When I write this function below it work: a <- c(1,2) for(i in a){ print(i) } [1] 1 [1] 2 But if I write this function with line space it dont work. a <- c(1,2) for(i in a){ print(i) } Error: unexpected '}' in "}" I remember that in early R version this limitation dont exist. This is a bus, a rule or a

Hi, It is possible to use a shell command inside a R script? I'm write a R script and I like to put somes shell commands inside to R. Somethink like: convert fig01.png fig01.xpm or sed ..., etc. It is possible? How? Thanks Ronaldo -- Acalme-se, s??o somente 0's e 1's -- > Prof. Ronaldo Reis J?nior | .''`. UNIMONTES/Depto. Biologia Geral/Lab. de Ecologia | : :' :

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 required package: cairoDevice > > playwith(plot(1,1)) Error in plot.new() : figure margins

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, so most people don't recognize them. -- > Prof. Ronaldo Reis J?nior | .''`.

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 Biologia Computacional | : :' : Campus Universit?rio Prof. Darcy Ribeiro, Vila Mauric?ia |

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)) baixo medio alto 11.56667 20.00000 26.80000 How to make the levels of pH ordered by this mean to the result to

Hi, I have a dummy question. Suppose that I have two explanatory variable, T1 (A, B) and T2 (C, D) and one response variable. > attach(dados) > tapply(Y,list(T1,T2),mean) C D A 2.200000 10.20000 B 2.223333 20.26667 In this case, "A" and "B" inside "C" have no difference, but have differences inside "D" I make this model: > m

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

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,"qvartemp",as.expression(mma$formula))) Look this new object: >

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

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 have a R Package. It is working, but I can't that it install the dependencies automatically. When I try to install I receive this message: > install.packages("./RT4Bio_1.0.tar.gz",dependencies=TRUE) Installing package(s) into '/usr/local/lib/R/site-library' (as 'lib' is unspecified) inferring 'repos = NULL' from the file name ERROR: dependency

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 device. How I made to put color off in trellis device? Thanks Ronaldo -- 2ª lei - Na