[R-sig-eco] simple quetion about posthoc in adonis

2013-03-01 Thread Paolo Piras
I saw in some forum some considerations about this topic. Maybe I propose a brutal (perhaps incorrect?) solution; given more than two groups in a model subjected to permutation in adonis.and if the global model is significant if one contrasts pairwise all possible pair comparisons

[R-sig-eco] multivariate continuous response and ordinal predictor(s)

2013-05-08 Thread Paolo Piras
Dear all, I need to run a model with multivariate continuous responses and one (or more) ORDINAL (i.e. 1,2,3, etc.) predictor variables; these are not factors because are ordinal; the more intuitive solution could be to apply a standard lm() but I ask you if some more appropriate strategies

Re: [R-sig-eco] simple question about nmds

2013-06-10 Thread Paolo Piras
Hi, actually I think you CAN use them (as predictors or response either) because, albeit in its unique way, NMDS **IS** a method to summarize a highly multidimensional phenomenon. best paolo Da: r-sig-ecology-boun...@r-project.org

[R-sig-eco] vegan RsquareAdj() for lm models

2013-10-03 Thread Paolo Piras
Dear list, I would like to easily compute the adjusted R-square in a lm model without intercept (excluding the intercept is essential for my analysis) I found that RsquareAdj() in vegan returns NA if the argument is a multiple-multivariate lm model thus including multivariate responses and

Re: [R-sig-eco] vegan RsquareAdj() for lm models

2013-10-03 Thread Paolo Piras
...@oulu.fi Inviato: giovedì 3 ottobre 2013 14.25 A: Paolo Piras; r-sig-ecology@r-project.org Oggetto: RE: [R-sig-eco] vegan RsquareAdj() for lm models Specific reason is that nobody has implemented this. These things don't come by automatic writing, but somebody must do them. What would you expect

Re: [R-sig-eco] vegan RsquareAdj() for lm models

2013-10-03 Thread Paolo Piras
to have non-centred RDA: I care too much about people and i don't want to do evil. If you really *know* that you need non-centred RDA. then you know how to change those lines of code in rda.default. Cheers, Jari Oksanen From: Paolo Piras [paolo.pi

[R-sig-eco] envfit() in vegan

2013-10-16 Thread Paolo Piras
Dear list, I write you because I do not understand the behavior of envfit() in vegan. Basically, it takes a matrix coming from an ordination procedure and it fits on it another matrix (often an environmental matrix). The projections of points onto vectors have maximum correlation with

Re: [R-sig-eco] envfit() in vegan

2013-10-17 Thread Paolo Piras
environmental variables and the ordination scores best paolo Da: Jari Oksanen jari.oksa...@oulu.fi Inviato: giovedì 17 ottobre 2013 9.32 A: Paolo Piras Cc: r-sig-ecology@r-project.org Oggetto: Re: [R-sig-eco] envfit() in vegan On 17/10/2013, at 02:49 AM, Paolo

Re: [R-sig-eco] Comparing results of two CCAs

2014-03-14 Thread Paolo Piras
Hi, maybe partial least squares: you can run two separate partial least squares analyses and then comparing vectors. best paolo Da: r-sig-ecology-boun...@r-project.org r-sig-ecology-boun...@r-project.org per conto di Eliot Miller eliotmil...@umsl.edu

Re: [R-sig-eco] Comparing results of two CCAs

2014-03-14 Thread Paolo Piras
...@gmail.com eliot.is...@gmail.com per conto di Eliot Miller eliotmil...@umsl.edu Inviato: venerd? 14 marzo 2014 15.19 A: Paolo Piras; highs...@highstat.com Cc: r-sig-ecology@r-project.org Oggetto: Re: [R-sig-eco] Comparing results of two CCAs The partial least squares sounds really promising

Re: [R-sig-eco] correlation between dissimilarity matrices

2014-05-21 Thread Paolo Piras
Hi Carlos Very briefly... you could look at rda() in vegan (not that in calibrate) partial and constrianed rda seems to be appropriate for your case. Maybe you should code a bit for controlling your nested stratification variable even if it is not very clear to me the nature of the data (sites vs

[R-sig-eco] 3d box plot

2015-11-08 Thread Paolo Piras
Hi folks, anyone could address me towards a R function/package able to do a 3d boxplot similar to http://158.132.155.107/oess/POSH/StatSoft/popups/popup125.gif Thanks in advance best paolo ___ R-sig-ecology mailing list R-sig-ecology@r-project.org

Re: [R-sig-eco] 3d box plot

2015-11-10 Thread Paolo Piras
times (the third dimension) best paolo Da: Jeremy Chacon <chaco...@umn.edu> Inviato: luned� 9 novembre 2015 16.21 A: Paolo Piras Cc: R-SIG list; r-sig-ecology@r-project.org Oggetto: Re: [R-sig-eco] 3d box plot This doesn't directly answer your question,

[R-sig-eco] 3d fitting plane

2016-01-19 Thread Paolo Piras
Hi folks, I look for a fast way to estimate a 3d fitting plane to my 3d data. I do not want z~y+x as this is a regression model. I just want the equation of the best plane that fits the data in 3d. Maybe using princomp() and Total least squares? Looking around I found some solutions but nothing

Re: [R-sig-eco] 3d fitting plane

2016-01-19 Thread Paolo Piras
pr$rotation[3,3]),alpha=0.5,col=3,bbox=F) decorate3d() It seems fine; any advice is welcome best paolo Da: R-sig-ecology <r-sig-ecology-boun...@r-project.org> per conto di Paolo Piras <paolo.pi...@uniroma3.it> Inviato: martedì 19 gennaio 20