[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 Jari Oksanen
From: r-sig-ecology-boun...@r-project.org [r-sig-ecology-boun...@r-project.org] on behalf of Paolo Piras [paolo.pi...@uniroma3.it] Sent: 03 October 2013 14:59 To: r-sig-ecology@r-project.org Subject: [R-sig-eco] vegan RsquareAdj() for lm models Dear list, I would like to easily compute the adjusted

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 Jari Oksanen
Subject: RE: [R-sig-eco] vegan RsquareAdj() for lm models Thankyou very much Jari! I think that it is nearly ok what I would like to have is the same as in RsquareAdj(vegan::rda(yy,xx)) that is a GLOBAL measure of the association BUT...I want it for a multiple-multivariate lm model that does

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

2013-10-03 Thread Paolo Piras
Piras; r-sig-ecology@r-project.org Oggetto: RE: [R-sig-eco] vegan RsquareAdj() for lm models Paolo, See ?RsquareAdj for the call interface. The default method can be called as RsquareAdj(x, n, m), and in the default method x is the unadjusted correlation, n is the number of observations and m