Re: [R] Abundance data ordination in R

2007-04-01 Thread José Rafael Ferrer Paris
There are many ways to do this, really. For example if you use constrained (~ canonical) correspondence analysis the distance measure between sites is Chi-square and absences are not informative to the analysis. Or you can use an ecological distance measure (similarity indices like Soerensen,

Re: [R] cex in xlab, ylab and zlab of persp

2007-03-15 Thread José Rafael Ferrer Paris
I had similar problems, actually it is very difficult to customize persp graphics, you should try wireframe in lattice instead. This reference might help to customize the wireframe plots: http://www.polisci.ohio-state.edu/faculty/lkeele/3dinR.pdf JR El mié, 14-03-2007 a las 20:40 -0700, Joseph

[R] segfault with correlation structures in nlme

2007-03-13 Thread José Rafael Ferrer Paris
Hi out there, I am trying to fit a species accumulation curve (increase in number of species known vs. sampling effort) for multiple regions and several bootstrap samples. The bootstrap samples represent different arrangements of the actual sample sequence. I fitted a series of nlme-models and

Re: [R] How to open more windows to make more graphs at once!

2007-03-07 Thread José Rafael Ferrer Paris
Dear Monireh, try using lattice: library(lattice) set.seed(1234) dat - data.frame(months=rep(1:10,80),upper = rnorm(800)+1, lower = rnorm(800)-1, observed = rnorm(800), best.sim = rnorm(800), stname = factor(gl(80, 10))) jpeg(filename =

Re: [R] anova applied to a lme object

2007-03-07 Thread José Rafael Ferrer Paris
The variances of the random effects and the residual variances are given by the summary function. Maybe VarCorr or varcomp gives you the answer you are looking for: library(nlme) library(ape) ?VarCorr ?ape JR El mié, 07-03-2007 a las 13:09 +0100, Berta escribió: Hi R-users, when carrying out

Re: [R] No years() function?

2007-03-07 Thread José Rafael Ferrer Paris
From the help of weekdays: Note: Other components such as the day of the month or the year are very easy to compute: just use 'as.POSIXlt' and extract the relevant component. Yet another option: help(package=chron) JR El mié, 07-03-2007 a las 15:35 +, Sérgio Nunes

Re: [R] Is there a quick way to count the number of times each element in a vector appears?

2007-03-06 Thread José Rafael Ferrer Paris
El lun, 05-03-2007 a las 22:16 -0800, Dylan Arena escribió: So here is my question in a nutshell: Does anyone have ideas for how I might efficiently process a matrix like that returned by a call to combinations(n, r, rep=TRUE) to determine the number of repetitions of each element in each row

Re: [R] Identifying points in a plot that have duplicate values

2007-03-06 Thread José Rafael Ferrer Paris
A very simple solution is given in: help(sunflowerplot,package=graphics) ##If you want to see it in action: example(sunflowerplot) El mar, 06-03-2007 a las 19:55 +1100, Jim Lemon escribió: David Lloyd wrote: I have code like this: -

Re: [R] The plot of qqmath

2007-03-06 Thread José Rafael Ferrer Paris
qqmath is in the lattice package, which is not compatible with the conventional graphics parameter. Why not using qqnorm and qqline instead? If you want to combine qqmath with other lattice plots you should look at the documentation of the lattice and grid packages. If you read carefully in

[R] different random effects for each level of a factor in lme

2007-03-05 Thread José Rafael Ferrer Paris
I have an interesting lme - problem. The data is part of the Master Thesis of my friend, and she had some problems analysing this data, until one of her Jurors proposed to use linear mixed-effect models. I'm trying to help her since she has no experience with R. I'm very used to R but have very