similar to: Very simple question on plot

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 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, 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 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, 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, 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. Anybody know a more easy or more generalized way to make this? Somethink like this

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 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 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 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, 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, after a debian upgrade my dev.copy2eps dont work anymore. I have this message: Error in matchFont(postscriptFonts(family)[[1]], old$encoding) : unknown font I try to change the font family in my .Rprofile: setHook(packageEvent("graphics", "onLoad"), function(...) { grDevices::ps.options(family="ComputerModern") } ) But it no fix

Hi, I need to add some different curves to a each panel in a xyplot. I have a old function to make this using panel.number, like this: panel=function(x,y,panel.number,...){ panel.xyplot(x,y,...) if(panel.number==1){ panel.curve(-655.8689+769.1589*log(5)+64.7981*log(x)-206.4475*log(5)^2) } if(panel.number==2){

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 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 try to make a model using lmer, but the weigths is not accept. m1<-lmer(ocup/total~tempo+(tempo|estacao),family=binomial,weights=total) Erro em lmer(ocup/total ~ tempo + (tempo | estacao), family = binomial, : object `weights' of incorrect type I dont understand why this error, with glm this work. the total object is a vector. Any idea? Thanks Ronaldo -- God is subtle, but

Hi, anybody have a JGR launcher for linux? Maybe a script that launch JGR directly without open R then library(JGR) and JGR(). Thanks Ronaldo -- Deflector shields just came on, Captain. -- > Prof. Ronaldo Reis J?nior | .''`. UNIMONTES/Depto. Biologia Geral/Lab. Ecologia Evolutiva | : :' : Campus Universit?rio Prof. Darcy Ribeiro, Vila Mauric?ia | `. `'` CP: 126, CEP:

Hi, I testing JGR and I like, but my ~ caracter dont work. My keyboard is Brazilian ABNT2. The key is OK, only in JGR it dont work. Anybody have any idea about this? Thanks Ronaldo -- Mais variado que baldea??o em Cacequi. -- > Prof. Ronaldo Reis J?nior | .''`. UNIMONTES/Depto. Biologia Geral/Lab. Ecologia Evolutiva | : :' : Campus Universit?rio Prof. Darcy Ribeiro, Vila

Hello all. I just switched from Windows to Mac and I am trying to find a good front end for R in Mac that is relatively user friendly and works well. I have used RWinEdt for several years in Windows, and something similar to RWinEdt for Mac would be ideal. I've looked into some of the emacs editor (e.g., Aquamacs) but I wanted to see what other people prefer, as I'm not sure I'm an