Skip to content
Prev 4097 / 7420 Next
Original Raw

vegan RsquareAdj() for lm models

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 multiple predictors, while it works for univariate response and multiple predictors. 

For example:
library(vegan)

yy<-matrix(rnorm(200,0,1),ncol=4)
xx<-matrix(rnorm(200,0,1),ncol=4)
RsquareAdj(lm(yy~xx-1))
RsquareAdj(lm(yy[,1]~xx-1))



There some specific reason for this behavior?
Thanks in advance for any advice
best regards
Poalo

Thread (5 messages)

Paolo Piras vegan RsquareAdj() for lm models Oct 3 Jari Oksanen vegan RsquareAdj() for lm models Oct 3 Paolo Piras vegan RsquareAdj() for lm models Oct 3 Jari Oksanen vegan RsquareAdj() for lm models Oct 3 Paolo Piras vegan RsquareAdj() for lm models Oct 3