search for: pausa

...t;wide", timevar="tr", idvar=c("code","sp" )) I'm getting very close. The only difference is in the 3rd row (**), that is when sp and code are the same I only get one record. Is there a way to get all records? Any idea? Thank you very much for any help Juli Pausas -- http://www.ceam.es/pausas

...ine each time. It is been running for 25 hours and I think it has done less than half of the file (Yes, probably I do not have a very good computer and I'm working under Windows ...). So my question is: do you know any other faster way to do this? Thanks in advance Juli -- http://www.ceam.es/pausas [[alternative HTML version deleted]]

...dom factor (C, the plot), using the aov(). So my first idea would be something like: aov(H ~ B * A %in% B * Time + Error(id) ) where id is the factor coded for the repeated subjects. But my question is how to include the random factor C ? Any help would be appreciated. Thanks Juli -- Juli G. Pausas CEAM & UA http://www.ceam.es/lass/pausas.htm

...d chisq.test() They may different results. E.g.: aa <- gl(2, 10) bb <- as.factor(c(1,2,2,2,1,2,1,2,2,2,1,2,2,2,1,1,1,2,1,1)) aa <- c(aa, aa) bb <- c(bb, bb) table(aa, bb) summary(table(aa, bb)) chisq.test(aa, bb) Could somebody give me some clues? Thanks in advance Juli -- Juli G. Pausas Centro de Estudios Ambientales del Mediterraneo (CEAM) C/ C.R. Darwin 14, Parc Tecnologic, 46980 Paterna, Valencia, SPAIN Tel: (+ 34) 96 131 8227; Fax: (+ 34) 96 131 8190 mailto:juli at ceam.es http://www.gva.es/ceam GCTE Fire Network - http://www.gva.es/ceam/FireNetwork -.-.-.-.-.-.-.-.-.-.-.-...

...but with an specific order of columns. Is there a direct way to do it? or I have to generate a new dataframe as follows: t <- data.frame(c=1:10, b=11:20, a=letters[1:10]) t2 <- data.frame(a=t$a, b=t$b, c=t$c) write.table(t2, row.names=F) Thanks for any comment Juli -- http://www.ceam.es/pausas

...end of strings in a vector? Something like the =TRIM functions of the OpenOffice spreadsheet. E.g., a <- c("hola ", "Yes ", "hello ") # I'd like to get: c("hola", "Yes", "hello") Thanks Juli -- http://www.ceam.es/pausas

Dear all, How could I use variables in a loop that their names are in a vector? For example: aaa <- 1:10 bbb <- aaa*2 ccc <- aaa+bbb varn <- c("aaa", "bbb", "ccc") m <- rep(NA, 3) for (i in 1:length(varn)) m[i] <- mean(varn[i]) # wrong thanks in advance Juli -- "Wars do not solve problems, wars generate even more problems"

...unused argument(s) (nsmall ...) > anova.negbin(m.nb) Error in get(x, envir, mode, inherits) : invalid first argument In addition: Warning message: tests made without re-estimating theta in: anova.negbin(m.nb) > Any suggestion of why I'm getting this error? Many thanks Juli -- Juli G. Pausas Centro de Estudios Ambientales del Mediterraneo (CEAM) C/ C.R. Darwin 14, Parc Tecnologic, 46980 Paterna, Valencia, SPAIN Tel: (+ 34) 96 131 8227; Fax: (+ 34) 96 131 8190 mailto:juli at ceam.es http://www.gva.es/ceam GCTE Fire Network - http://www.gva.es/ceam/FireNetwork -------------- next part...

Dear R users, - Is it possible to produce a figure with 2 plots that they share one of the axis, e.g., the y-axis? I did not succeed by setting mai[4] <- 0. Is there a simple way? - How could I convert 2 factor variables in a single factor variable which is the combination of the other 2. Example: lith: a factor 2 levels "ca", "ma" sp: a factor with 2 levels,

...s (and viceversa), than my function below (it works). Any help? m2coord <- function(m) { k <- nrow(m)*ncol(m) aa <- data.frame(r=1:k, c=1:k, v=1:k) k <- 0 for (i in 1:nrow(m)) for (j in 1:ncol(m)) { k <- k+1 aa$f[k]=i; aa$c[k]=j; aa$v[k]=m[i,j] } aa } Juli -- Juli G. Pausas Centro de Estudios Ambientales del Mediterraneo (CEAM) C/ C.R. Darwin 14, Parc Tecnologic, 46980 Paterna, Valencia, SPAIN Tel: (+ 34) 96 131 8227; Fax: (+ 34) 96 131 8190 mailto:juli at ceam.es http://www.gva.es/ceam GCTE Fire Network - http://www.gva.es/ceam/FireNetwork -- "Wars do not s...

...,28,32,37,34,35) tr <- gl(2, 4) sf <- cbind(s,f) # numbers of successes and failures summary(glm(sf ~ tr, family="binomial")) # tr ns sf[8,1] <- 1 summary(glm(sf ~ tr, family="binomial")) # tr significative ** Thanks for any suggestion Juli -- http://www.ceam.es/pausas

...plot do not fall in the midpoint of the bars in the first. Any trick for setting the 2 plots at the same scale? I have unsuccessfully tried: plot(t, data, type="b", xlim=c(0,7)) plot(t, data, type="b", xlim=c(min(t),max(t))) (R 1.8.1, for Windows) Thanks Juli -- Juli G. Pausas Centro de Estudios Ambientales del Mediterraneo (CEAM) C/ Charles R. Darwin 14, Parc Tecnologic, 46980 Paterna, Valencia, SPAIN Tel: (+ 34) 96 131 8227; Fax: (+ 34) 96 131 8190 mailto:juli at ceam.es http://www.gva.es/ceam GCTE Fire Network - http://www.gva.es/ceam/FireNetwork

Dear colleges, I do not understand the following behaviour: > aa <- data.frame(a1= 1:10, a2= c(rep(NA, 5), 1:5) ) > aa a1 a2 1 1 NA 2 2 NA 3 3 NA 4 4 NA 5 5 NA 6 6 1 7 7 2 8 8 3 9 9 4 10 10 5 > aa[!aa$a2==1, ] # removing rows with a2==1 a1 a2 NA NA NA NA.1 NA NA NA.2 NA NA NA.3 NA NA NA.4 NA NA 7 7 2 8 8 3 9 9 4 10 10 5 I didn't

...# Cuadrado Original quads3d(v[,1],v[,2],v[,3],col="green") } # Rotacion v[1,] <- rotate3d(v[1,], ang, 0, 1, 0) v[2,] <- rotate3d(v[2,], ang, 0, 1, 0) v[3,] <- rotate3d(v[3,], ang, 0, 1, 0) v[4,] <- rotate3d(v[4,], ang, 0, 1, 0) # Pausa Sys.sleep(0.3) # Limpiar if(rgl.cur() !=0) rgl.clear() # Nuevo cuadrado quads3d(v[,1],v[,2],v[,3],col="green") # cont <- cont + 1 } ¿Hay alguna forma de dar movimiento a una figura con rgl sin que aparezcan esos flashes? Muchas gracias por adelantado...

...= Y, : inner loop 1; can't correct step size In addition: Warning message: Step size truncated due to divergence in: (if (is.empty.model(mt)) glm.fit.null else glm.fit)(x = X, y = Y, (this error message doesn't appear with family= gaussian) Juli ******************** Juli G. Pausas Centro de Estudios Ambientales del Mediterraneo (CEAM) C/ C.R. Darwin 14 Parc Tecnologic, 46980 Paterna, Valencia, SPAIN Tel: (+ 34) 96 131 8227 Fax: (+ 34) 96 131 8190 E-mail: juli at ceam.es http://www.gva.es/ceam -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r...

Hi all,, Is it possible to set a thicker line for the box around the plots? (i.e. for the four axes) Something like lwd (lines) but for the box. Thanks juli -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the

...9; 0.05 `.' 0.1 ` ' 1 (Dispersion parameter for poisson family taken to be 1) Null deviance: 40605 on 283 degrees of freedom Residual deviance: 32881 on 282 degrees of freedom AIC: 35207 Number of Fisher Scoring iterations: 3 Thank you very much for your help Juli -- Juli G. Pausas Centro de Estudios Ambientales del Mediterraneo (CEAM) C/ C.R. Darwin 14, Parc Tecnologic, 46980 Paterna, Valencia, SPAIN Tel: (+ 34) 96 131 8227; Fax: (+ 34) 96 131 8190 mailto:juli at ceam.es http://www.gva.es/ceam GCTE Fire Network - http://www.gva.es/ceam/FireNetwork -------------- next part...

...9; 0.05 `.' 0.1 ` ' 1 (Dispersion parameter for poisson family taken to be 1) Null deviance: 40605 on 283 degrees of freedom Residual deviance: 32881 on 282 degrees of freedom AIC: 35207 Number of Fisher Scoring iterations: 3 Thank you very much for your help Juli -- Juli G. Pausas Centro de Estudios Ambientales del Mediterraneo (CEAM) C/ C.R. Darwin 14, Parc Tecnologic, 46980 Paterna, Valencia, SPAIN Tel: (+ 34) 96 131 8227; Fax: (+ 34) 96 131 8190 mailto:juli at ceam.es http://www.gva.es/ceam GCTE Fire Network - http://www.gva.es/ceam/FireNetwork -------------- next part...

Hi, This is a very basic question, but I would like to undestand hist(). I thought that the hist( , freq=FALSE) should provide the relative frequencies (probabilities), and so they should sum 1, however: set.seed(2) ah <- hist(rnorm(100), freq=F) sum(ah$intensities) [1] 2 set.seed(2) bh <- hist(rlnorm(100), freq=F) sum(bh$intensities) [1] 0.4999996 I'm getting similar figures with

Dear all, I'd like to subset a df within a function, and use select for choosing the variable. Something like (simplified example): mydf <- data.frame(a= 0:9, b= 10:19) ttt <- function(vv) { tmpdf <- subset(mydf, select= vv) mean(tmpdf$vv) } ttt(mydf$b) But this is not the correct way. Any help? Thanks in advance Juli