Re: [R] Question on CCA and RDA analysis

2015-04-22 Thread Gavin Simpson
[[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guid

Re: [R] package mgcv - predict with bam: Error in X[ind, ] :, subscript out of bounds

2014-02-03 Thread Gavin Simpson
roject.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. -- Gavin Simpson, PhD __ R-help@r-

Re: [R] What computer power for GAMM models

2013-10-25 Thread Gavin Simpson
> 56 150 BAUD > > FRANCE > tél : 02 97 28 80 38 > fax: 02 97 28 81 10 > [[alternative HTML version deleted]] > > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do

Re: [R] Assessing temporal correlation in GAM with irregular time steps

2013-09-05 Thread Gavin Simpson
That is not that useful as you need to give the software something about when the samples occur in time, otherwise it doesn't have the information needed to properly model the decay in correlation with time. You need to give it the observation time, however you measured it. HTH G

Re: [R] Setting up 3D tensor product interactions in mgcv

2013-09-03 Thread Gavin Simpson
? > > Does anyone have experience with this, and can explain how I should do > it properly? > > Mark > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http:

Re: [R] Assessing temporal correlation in GAM with irregular time steps

2013-09-03 Thread Gavin Simpson
he 'distance' argument > > Best wishes > Tom > > > > [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting gui

Re: [R] Multiple interaction terms in GAMM-model

2013-07-26 Thread Gavin Simpson
m the R help mailing list archive at Nabble.com. > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide co

Re: [R] MGCV: Degrees of freedom of smooth terms

2013-07-25 Thread Gavin Simpson
up on this) for a (far) more authoritative answer on this part of your question. HTH G > > Thank you very much for answering! > > Best wishes, > Christoph > > -- Gavin Simpson, PhD [t] +1 306 337 8863 Adjunct Professor, Department of Biol

Re: [R] Error with metaMDS

2013-06-24 Thread Gavin Simpson
> X.dist1 <- metaMDSdist(X, method="bray",zerodist = "ignore") > X.dist2 <- metaMDSdist(X, method="bray",zerodist = "add") > > Can Please help me with this. > Thanks, > Mitra > > [[alternative HTML version deleted]] > > __

Re: [R] gamm in mgcv random effect significance

2013-06-11 Thread Gavin Simpson
> not support a direct model comparison to test the significance of the > random effects. I'm guessing the best strategy based on the summary.gam > is probably just to compare fit indices like Log Likelihoods. > > If anyone has any other suggestions, though, please do let me know.

Re: [R] gamm in mgcv random effect significance

2013-06-07 Thread Gavin Simpson
t; Formula: ~1 | g/xc > Parameter estimate(s): > Phi > 0.1809409 > Variance function: > Structure: fixed weights > Formula: ~invwt > Number of Observations: 264 > Number of Groups: > g xc %in% g > 134 > > $gam > > F

Re: [R] Vegan package treatment of zeros

2013-05-15 Thread Gavin Simpson
logy/attachments/20130419/114d1e5f/attachment.pl HTH G -- Gavin Simpson, PhD [t] +1 306 337 8863 Adjunct Professor, Department of Biology[f] +1 306 337 2410 Institute of Environmental Change & Society [e] gavin.simp...@uregina.ca 523 Research and Innovation Centre

Re: [R] NMDS using Vegan

2013-05-15 Thread Gavin Simpson
** performed an NMDS. Not sure why Jari labelled these as "MDSx". If this bothers you so, add your own labels: require("vegan") data(dune) sol <- metaMDS(dune) plot(sol, xlab = "NMDS1", ylab = "NMDS2") HTH G -- Gavin Simpson, PhD

Re: [R] NMDS in Vegan: problems in stressplot, best solution

2013-04-26 Thread Gavin Simpson
ations" are random starts, each of which has iterations of the algorithm. HTH G > parth.mds <- metaMDS(WorldPRSenv, distance = "bray", k = 2, trymax = 100, > engine = c("monoMDS", "isoMDS"), > autotransform =TRUE, wascores = TRUE, expand = TRUE, tra

Re: [R] RDA permutest envfit

2013-04-24 Thread Gavin Simpson
c à Rimouski > > remilesmeri...@yahoo.ca > > > > [[alternative HTML version deleted]] > > > > __ > > R-help@r-project.org mailing list > > https://stat.ethz.ch/mailman/listinfo/r-help > > PLEASE do

Re: [R] RDA permutest envfit

2013-04-24 Thread Gavin Simpson
..@yahoo.ca > > [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and pro

Re: [R] GAM Penalised Splines - Intercept

2013-04-23 Thread Gavin Simpson
arametric factor terms what the intercept is depends on the contrasts used - the default in R would represent the mean of the response for the reference level(s) of factor terms in the model. The intercept is separates from the penalised spline terms in the model. HTH G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~

Re: [R] CCA report "species environment correlation"

2013-04-13 Thread Gavin Simpson
> and provide commented, minimal, self-contained, reproducible code. -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Simpson [t] +44 (0)20 7679 0522 ECRC, UCL Geography, [f] +44 (0)20 7679 0565 Pearson Building, [e] gavin.simpsonATNOSPAMucl.ac.uk Gower Street, London

Re: [R] Error in plot.envfit(ef, p.max = 0.1) : (subscript) logical subscript too long

2012-12-31 Thread Gavin Simpson
t; > > > > > species<-metaMDS(species_matrix)ef<-envfit(species,environmentaldata_file,permu=999,na.rm=TRUE)efplot(species, > dis="sites")plot(ef,p.max=0.1) > > > > > Error in plot.envfit(ef, p.max = 0.1) : (subscript) logical > > > subscript too long > > > > > Many thanks! > > > > > Laura > > &g

Re: [R] Error in plot.envfit(ef, p.max = 0.1) : (subscript) logical subscript too long

2012-12-30 Thread Gavin Simpson
nger than the object being subscripted. I promise to delete the data once I have tracked the problem down. Cheers, G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Simpson [t] +44 (0)20 7679 0522 ECRC, UCL Geograph

Re: [R] CCA plot

2012-11-28 Thread Gavin Simpson
rs[,2] * mul, length = 0.05, col = "blue") text(scrs * mul * 1.1, labels = rownames(scrs), cex = 0.9, col = "blue") Not that I recommend this "pruning" a CCA... HTH G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Si

Re: [R] CCA with Vegan - Plot problem

2012-11-09 Thread Gavin Simpson
abble.com. > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, r

Re: [R] CCA with Vegan - Plot problem

2012-11-09 Thread Gavin Simpson
R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. G

Re: [R] labeling loading vectors in vegan

2012-07-31 Thread Gavin Simpson
version deleted]] > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, s

Re: [R] reception of (Vegan) envfit analysis by manuscript reviewers

2012-05-10 Thread Gavin Simpson
tions! > > > Matt > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, se

Re: [R] conducting GAM-GEE within gamm4?

2012-05-08 Thread Gavin Simpson
g the following warning message > after model completion for larger models (with many more variables > than 2): > > Warning message: > In mer_finalize(ans) : false convergence (8) > > > > > Thank you for your time, effort, and patience dealing with a novice. > > Sin

Re: [R] Clustering analysis with ordination plots

2012-05-02 Thread Gavin Simpson
ve at Nabble.com. > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-containe

Re: [R] Long command in Sweave

2012-04-13 Thread Gavin Simpson
suitable editor which can break the > command properly. > > All the best > > On 13 April 2012 19:33, Gavin Simpson wrote: > > On Fri, 2012-04-13 at 17:46 +0800, Wincent wrote: > >> Dear useRs, > >> > >> I am writing a vignette for a package,

Re: [R] Long command in Sweave

2012-04-13 Thread Gavin Simpson
ke lines where it could. In short get out of the habit of writing long lines of R code; you'll be better in the long run laying your code out logically. HTH G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Simpson [t] +

Re: [R] customizing plot, Vegan CCA

2012-04-12 Thread Gavin Simpson
ted, minimal, self-contained, reproducible code. -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Simpson [t] +44 (0)20 7679 0522 ECRC, UCL Geography, [f] +44 (0)20 7679 0565 Pearson Building, [e] gavin.simpsonATNOSPAMucl.ac.uk Go

Re: [R] Gamm and post comparison

2012-02-29 Thread Gavin Simpson
ect.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~

Re: [R] z-transform each column of a data.frame

2012-01-20 Thread Gavin Simpson
project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%

Re: [R] any DCCA function in R?

2011-12-20 Thread Gavin Simpson
r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Simpson [t] +44 (0)20 7679 0522 ECRC

Re: [R] I cannot get species scores to plot with site scores in MDS when I use a distance matrix as input. Problems with NA's?

2011-11-29 Thread Gavin Simpson
tested as there wasn't reproducible code, but hopefully you'll be able to work it out from the above. HTH G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Simpson [t] +44 (0)20 7679 0522 ECRC, UCL Geography, [f] +44 (0)20 7679

Re: [R] I cannot get species scores to plot with site scores in MDS when I use a distance matrix as input. Problems with NA's?

2011-11-29 Thread Gavin Simpson
contained, reproducible code. > > > > > -- > View this message in context: > http://r.789695.n4.nabble.com/I-cannot-get-species-scores-to-plot-with-site-scores-in-MDS-when-I-use-a-distance-matrix-as-input-Pr-tp4103699p4104406.html > Sent from the R help mailing list archive at Na

Re: [R] Vegan: Diversity Plot, label points

2011-11-29 Thread Gavin Simpson
ww.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Simpson [t] +44 (0)20 7679 0522 ECRC, UCL Geography, [f] +44 (0)20 7679 0565 Pearson Building, [e] g

Re: [R] Permutations

2011-11-18 Thread Gavin Simpson
hoose(150, 10) [1] 1.169554e+15 ?? G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Simpson [t] +44 (0)20 7679 0522 ECRC, UCL Geography, [f] +44 (0)20 7679 0565 Pearson Building, [e] gavin.simpsonATNOSPAMucl.ac.uk Gower Stre

Re: [R] cca with repeated measures

2011-11-18 Thread Gavin Simpson
he code yourself using functions from permute). HTH G > thanks, > René > > -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Simpson [t] +44 (0)20 7679 0522 ECRC, UCL Geography, [f] +44 (0)20 7679 0565 Pearson Building, [e] gavin.simpsonATNOSPAMucl.

Re: [R] cca with repeated measures

2011-11-18 Thread Gavin Simpson
in advance, > René > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, re

Re: [R] permutation within rows of a matrix

2011-11-17 Thread Gavin Simpson
cable any that don't apply. ;-) > Peter Ehlers > > On 2011-11-16 12:12, Gavin Simpson wrote: > > On Wed, 2011-11-16 at 14:29 -0500, R. Michael Weylandt wrote: > >> Suppose your matrix is called X. > >> > >> ? sample > >> X[sample(nrow(X))

Re: [R] permutation within rows of a matrix

2011-11-16 Thread Gavin Simpson
SE do read the posting guide http://www.R-project.org/posting-guide.html > > and provide commented, minimal, self-contained, reproducible code. > > > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-h

Re: [R] ordination in vegan: what does downweight() do?

2011-11-08 Thread Gavin Simpson
As a general matter of routine practice, I don't think this is a very good way of working, especially as we have no good ecological grounds for doing so and who knows what information these species could be telling us if we just listened to them instead of discarding them. HTH G -- %~%~%~

Re: [R] What does \Sexpr[results=rd]{} exactly mean in Rd?

2011-10-17 Thread Gavin Simpson
.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. G

Re: [R] binomial GLM quasi separation

2011-10-17 Thread Gavin Simpson
-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~

Re: [R] plotting issues with PCA

2011-10-17 Thread Gavin Simpson
t;HAL7","DOH1","DOH2","DOH3","DOH4","DOH5") > > fish.pca <-rda(fish.sqrt.h) > > fish.site <- scores(fish.pca,display="sites",scaling=3) > > fish.spp <- > scores(fish.pca,display="species",scali

Re: [R] pdIdent in smoothing regression model

2011-10-17 Thread Gavin Simpson
__ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~

Re: [R] pdIdent in smoothing regression model

2011-10-17 Thread Gavin Simpson
__ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. -- %~%~%~%~%

Re: [R] binomial GLM quasi separation

2011-10-14 Thread Gavin Simpson
in this. > > -- > View this message in context: > http://r.789695.n4.nabble.com/binomial-GLM-quasi-separation-tp3901687p3904261.html > Sent from the R help mailing list archive at Nabble.com. > > __ > R-help@r-project.org mailin

Re: [R] Question about GAMs

2011-10-13 Thread Gavin Simpson
about-GAMs-tp3900848p3901842.html > Sent from the R help mailing list archive at Nabble.com. > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.or

Re: [R] Vegan: Anova.CCA accessing original data using option by="margin"

2011-10-11 Thread Gavin Simpson
chive at Nabble.com. > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, se

Re: [R] About stepwise regression problem

2011-10-07 Thread Gavin Simpson
hz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Simpson [t] +44 (0)20

Re: [R] Adonis and nmds help and questions for a novice.

2011-10-04 Thread Gavin Simpson
___ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%

Re: [R] Adonis and nmds help and questions for a novice.

2011-10-04 Thread Gavin Simpson
ne who can help me as I have asked people who > use R near me and no-one is sure and uses Adonis.. > > Ash > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide h

Re: [R] Interaction plot type=o

2011-10-03 Thread Gavin Simpson
symbols of > >>>> the points and the line. > >>>> > >>>> > >>> If it's OK to have the lines going right though the symbols, then > >>> go ahead, > >>> hack the code. All you need to do is add ,"o" to the type > >>> argum

Re: [R] GAMs in R : How to put the new data into the model?

2011-09-29 Thread Gavin Simpson
ing guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Simpson [t] +44 (0)20 7679 0522 ECRC, UCL Geography, [f] +44 (0)20 7679 0565 Pearson Building, [e] gavin.simpso

Re: [R] metaMDS

2011-09-28 Thread Gavin Simpson
ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Simpson [t] +44 (0

Re: [R] problem with predict.gam in package gam

2011-09-15 Thread Gavin Simpson
; I am using R 2.13.1 on a PC inside RStudio. > > I realize I can do this in package mgcv, but just curious as to > whether this is something known and fixable? > > Thanks > -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Simpson

Re: [R] envfit vector labels with ordiplot3d

2011-09-13 Thread Gavin Simpson
oices = c(1, 2)..." which seems to suggest the function might not > allow a third axis at the moment? > Cheers, > Briony > > > On Tue, Sep 13, 2011 at 2:53 AM, Gavin Simpson [via R] < > ml-node+s789695n3807984...@n4.nabble.com> wrote: > > > On Mon, 2011-09-12

Re: [R] On-line machine learning packages?

2011-09-13 Thread Gavin Simpson
_ > > >>> r-h...@r-project.org mailing list > > >>>https://stat.ethz.ch/mailman/listinfo/r-help > > >>> PLEASE do read the posting > > >>> guidehttp://www.R-project.org/posting-guide.html > > >>

Re: [R] mvpart analyses with covariables

2011-09-13 Thread Gavin Simpson
ative HTML version deleted]] > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-co

Re: [R] envfit vector labels with ordiplot3d

2011-09-12 Thread Gavin Simpson
g guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > > -- > View this message in context: > http://r.789695.n4.nabble.com/envfit-vector-labels-with-ordiplot3d-tp3800669p3807015.html > Sent from the

Re: [R] Symbol Font Baseline, Cairo, Card Symbols

2011-08-17 Thread Gavin Simpson
ever. G > /iaw > > > Ivo Welch (ivo.we...@gmail.com) > > > > > > On Wed, Aug 17, 2011 at 5:18 AM, Gavin Simpson > wrote: > > On Tue, 2011-08-16 at 17:46 -0700, ivo welch wrote: > >> I think I found a bug in the Cairo library, plus weird beha

Re: [R] Symbol Font Baseline, Cairo, Card Symbols

2011-08-17 Thread Gavin Simpson
___ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Simpson

Re: [R] packfor

2011-08-17 Thread Gavin Simpson
could find it and if I need a special R > version to do it. packfor is part of the *sedar* project on R-forge. You can get a binary or sources from: http://r-forge.r-project.org/R/?group_id=195 HTH G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Si

Re: [R] reflecting a PCA biplot

2011-08-10 Thread Gavin Simpson
m = ylim, xlim = xlim, ylab = "PC2", xlab = "PC1", cex = 0.7, asp = 1, type = "n") ## no plotting this time first abline(h = 0, lty = "dotted") abline(v = 0, lty = "dotted") text(si.scrs[,1], si.scrs[,2], labels = rownames(si.scrs), cex = 0.7) ## add site and species scores as labels text(sp

Re: [R] Running R in a sandbox

2011-08-03 Thread Gavin Simpson
make sure you turn on the email lists and off the functions/vignettes. > Thank you. G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Simpson [t] +44 (0)20 7679 0522 ECRC, UCL Geography, [f] +44 (0)20 7679 0565 Pearson Building,

Re: [R] Help with R

2011-08-02 Thread Gavin Simpson
t does predict() not do that you are expecting? Indeed. G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Simpson [t] +44 (0)20 7679 0522 ECRC, UCL Geography, [f] +44 (0)20 7679 0565 Pearson Building, [e] gavin.simpsonATNOSPAMucl.ac.uk Gower Street, Londo

Re: [R] Elegant way to subtract matrix from array

2011-07-27 Thread Gavin Simpson
y(1:100,dim=c(2,5,10)) > > Q <- matrix(1:10,nrow=2,ncol =5) For this one, consider the often overlooked function `sweep()`: sweep(Z, c(1,2), Q, "-") does what you wanted. c(1,2) is the `MARGIN` argument over the dimensions that Q will be swept from. HTH G -- %~%~%~%~%~%~%~%~%

Re: [R] function lm, get back the coefficient

2011-07-26 Thread Gavin Simpson
t;) > > unname(z) > [1] 2 > > z > a > 2 > > zz <- unname(z) > > zz > [1] 2 > > names(z) <- NULL > > z > [1] 2 > > Cheers, > Bert > > > On Tue, Jul 26, 2011 at 1:18 PM, Gavin Simpson > wrote: > > On Tue

Re: [R] function lm, get back the coefficient

2011-07-26 Thread Gavin Simpson
ry names(result)=NULL to get rid of the name. > > Barry `unname()` exists for this purpose: > foo <- 1 > names(foo) <- "bar" > foo bar 1 > unname(foo) [1] 1 G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Simpson

Re: [R] multiple plots in single frame: 2 upper, 1 lower

2011-07-21 Thread Gavin Simpson
t.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~

Re: [R] FW: lasso regression

2011-07-12 Thread Gavin Simpson
.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Simpson [t] +44 (0)20 7679 0522 ECRC,

Re: [R] PERMANOVA+ and adonis in vegan package

2011-07-11 Thread Gavin Simpson
> R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code.

Re: [R] relative euclidean distance

2011-07-07 Thread Gavin Simpson
Pass the returned object to dist() to get the distances you want. HTH G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Simpson [t] +44 (0)20 7679 0522 ECRC, UCL Geography, [f] +44 (0)20 7679 0565 Pearson Building, [e] gavin.

Re: [R] permil symbol linux

2011-07-06 Thread Gavin Simpson
f Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ > > >> University of Oxford, Tel: +44 1865 272861 (self) > > >> 1 South Parks Road, +44 1865 272866 (PA) > > >> Oxford OX1 3TG, UKFax: +44 1865 272595 > > > > > > > > &g

Re: [R] Standards for delivery of GPL software in CRAN packages

2011-06-27 Thread Gavin Simpson
> .tar.gz files from those links should have full R, C and Fortran source code. > > I think we need counter-examples... > > Barry > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do

Re: [R] Setting up a State Space Model in dlm

2011-06-10 Thread Gavin Simpson
rm=c(0,0,0,0) , build=buildNAIRU, > hessian=TRUE, control=list(maxit=500))) > (dlmNAIRU <- buildNAIRU(fitNAIRU$par)) > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.

Re: [R] Problem with Princurve

2011-05-23 Thread Gavin Simpson
> > > This is even worse. It seems like the algorithm is quite sensitive to > starting values. Is this behavior expected or is there some flaw in the > algorithm? > > Ravi. > > > From: r-help-boun...@r-project.org [r-he

Re: [R] Problem with Princurve

2011-05-23 Thread Gavin Simpson
w.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Simpson [t] +44 (0)20 7679 0522 ECRC, UCL Geography, [f] +44 (0)20 7679 0565 Pearson Building, [e] gavin.simpsonATNOSPAMucl.ac.uk

Re: [R] metaMDS and envfit: Help reading output

2011-05-23 Thread Gavin Simpson
dimensions are these vectors > > being correlated? Is there code to produce the x, y, and z coordinates of > > each of the sites in Fish.mds? > > > > Thanks very much. > > Katie > > > > > >[[alternative HTML version deleted]] > > > >

Re: [R] R forum for only Statistics

2011-05-07 Thread Gavin Simpson
.org/posting-guide.html > > and provide commented, minimal, self-contained, reproducible code. > > > > > __ R-help@r-project.org mailing > list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting > guide

Re: [R] Using $ accessor in GAM formula

2011-05-06 Thread Gavin Simpson
> > > > Estimated degrees of freedom: > > 3.9393 total = 4.939297 > > > > GCV score: 0.002666985 > > > > > Just out of pure academic interest. :-) > > > > Hope your academic curiosity is now satisfied. :) > > > > HTH. > > &g

Re: [R] principal components

2011-04-04 Thread Gavin Simpson
h("EOF") It is at times like this that an R equivalent of: http://lmgtfy.com/ would be handy ;-) Or it slightly cruder cousin... G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Simpson [t] +44 (0)20 7679 0522 ECRC, UCL Geography

Re: [R] Dirichlet surface

2011-03-30 Thread Gavin Simpson
> 2011-03-30 10:42 keltezéssel, Kehl Dániel írta: > > It helped a lot indeed, thank you very much! > > Now I understand why it was a problem for persp! > > > > Daniel > > > > 2011-03-30 10:31 keltezéssel, Gavin Simpson írta: > >> On Wed, 2011-03-30 at 09:55

Re: [R] Dirichlet surface

2011-03-30 Thread Gavin Simpson
org mailing list > >> https://stat.ethz.ch/mailman/listinfo/r-help > >> PLEASE do read the posting guide > >> http://www.R-project.org/posting-guide.html > >> and provide commented, minimal, self-contained, reproducible code. > > > > David Winsemius, MD > > Heri

Re: [R] Quick recode of -999 to NA in R

2011-03-30 Thread Gavin Simpson
0,"a1gpar"] <- NA > > dat0[dat0$a1pias==-999.,"a1pias"] <- NA > > dat0[dat0$a1devt==-999.,"a1devt"] <- NA > > [[alternative HTML version deleted]] > > ______ > R-he

Re: [R] ordination in vegan

2011-03-28 Thread Gavin Simpson
> > > > [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] take me off the R help mailing list.

2011-03-21 Thread Gavin Simpson
A "please" wouldn't go amiss there! You can remove yourself; the first link in the footer of every email from the list is to the mailing list website. From there you can unsubscribe from the list. If you encounter probelms, contact the list admin; the several thousand members of this list do no

Re: [R] exploring dist()

2011-03-20 Thread Gavin Simpson
elp@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~

Re: [R] persuade tabulate function to count NAs in a data frame

2011-03-19 Thread Gavin Simpson
je tudomásul veszi és > hozzájárul, hogy az üzenetekhez más banki alkalmazott is hozzáférhet az EBH > folytonos munkamenetének biztosítása érdekében. > > > This e-mail and any attached files are confidential and/...{{dropped:19}} > >

Re: [R] Problem with vegan package instalation on linux

2011-03-09 Thread Gavin Simpson
[[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented,

Re: [R] Rioja package, creating transfer function, WA, "Error in FUN"

2011-03-02 Thread Gavin Simpson
e=sqlFetch(STEST, "Sheet8") > > odbcClose(STEST) > > Spe > > (The species data contains the abundance of 32 species over 90 sites, set > out like this) >F1AmpFlavAmpWri ArcCat ArcDis > 1 N1_1 22.2929936 0.000 0.000

Re: [R] different results in MASS's mca and SAS's corresp

2011-02-07 Thread Gavin Simpson
reproduce SAS exactly. And I'm a little frustrated I > can't remember the matrix command to get that multiplication done > without cbinding the 2 columns together that way. You might have been thinking of sweep(): sweep(scalermca[,1:2], 2, c(0.827094,-0.7644828), "*&qu

Re: [R] Please stop all e-mail

2011-02-04 Thread Gavin Simpson
_ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. -- %~%~%~%~%~%~%~%~

Re: [R] GAM quasipoisson in MuMIn

2011-02-04 Thread Gavin Simpson
.html > and provide commented, minimal, self-contained, reproducible code. -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Simpson [t] +44 (0)20 7679 0522 ECRC, UCL Geography, [f] +44 (0)20 7679 0565 Pearson Building, [e] gavin.simpsonATNOSPAMucl

Re: [R] Help needed

2011-01-28 Thread Gavin Simpson
gt; > and provide commented, minimal, self-contained, reproducible code. > > > > [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do

Re: [R] Quasi-poisson glm and calculating a qAIC and qAICc...trying to modilfy Bolker et al. 2009 function to work for a glm model

2011-01-28 Thread Gavin Simpson
Instead of resorting to these ad-hoc approaches to correct for overdispersion, you would be better off fitting models that model the overdispersion. Try a negative binomial (glm.nb() in MASS) or the zeroinfl() and hurdle() functions in the pscl package. Those all have proper log likelihoods and you c

Re: [R] agnes clustering and NAs

2011-01-28 Thread Gavin Simpson
I was referring to is : > > "In case of a matrix or data frame, each row corresponds to an observation, > and each column corresponds to a variable. All variables must be numeric. > Missing values (NAs) are allowed." > > So, I'm under the impression it handles

Re: [R] agnes clustering and NAs

2011-01-27 Thread Gavin Simpson
iling list > > https://stat.ethz.ch/mailman/listinfo/r-help > > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > > and provide commented, minimal, self-contained, reproducible code. > > _

Re: [R] Extrapolating values from a glm fit

2011-01-27 Thread Gavin Simpson
etc), so you compute them on the link scale and apply the inverse of the link function to get them on to the response scale. HTH G > > > > On Wed, Jan 26, 2011 at 11:41 PM, Gavin Simpson > wrote: > On Wed, 2011-01-26 at 19:25 -1000, Ahnate Lim wrote: >

  1   2   3   4   5   6   7   8   >