similar to: Calculation of lags

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

Hi, I install rJava in linux. The installation, from source, is ok. But I have this error on load package. > library(rJava) Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to load shared library '/usr/local/lib/R/site-library/rJava/libs/rJava.so': libjvm.so: cannot open shared object file: No such file or directory Error: .onLoad failed in

Hi, I have R 2.3.1 installed by debian p?ckage. I install only the base and recommended R packages from Debian source, all others packages I install from the R source at CRAN. But, I try to use library MASS, but I received this message: > library(MASS) Error: This is R 2.3.1, package 'MASS' needs >= 2.4.0 What is the problem? Inte Ronaldo -- "It ain't over until

Hi, I write several xyplot graphics on a source file. When I try to use source("graphics.R") the source don't work, but if I use source("graphics.R",echo=T) it work. Generally some commands work without echo=T, but xyplot dont work. Why it dont work without echo=T? It is possible to write a source that dont need echo=T on source command? Thanks Ronaldo -- Aluno de

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 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 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 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, 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, 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, is possible to make a R script to run under a console without open the R environment? Something like this example.R #!/usr/bin/R function(name="Put here your name") { print(name) } In a console I make ./example.R name="Ronaldo Reis J?nior" then program print my name. It is possible? Thanks Ronaldo -- A jury consists of twelve persons chosen to decide who has the

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 made a survival analysis using survreg function. Its OK. But, I need to show in the graph 1-weibull curves. Ploting weibull curve a add the estimated points following these steps: fit <- survfit(Surv(Time,Censor)~Group) points(fit[1]) points(fit[2]) points(fit[3]) I have 3 groups. But I dont know how to plot the correct points on a inverse weilbull graph. I try to use

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