I really would like to be able to read about this in a document but I
cannot find my way around the documentation properly
Given the code...
M <- matrix(runif(5*20), nrow=20)
colnames(M) <- c('a', 'b', 'c', 'd', 'e')
ind <- c(1,2,3,4)
dep <- 5
I can then do...
l2 <- lm(M[,dep]~M[,ind]) ## Clearly not useful!
summary(l2)
I am not sure what my regression formula is.
The results are (edited for brevity)
> summary(l2)$coefficients
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.63842 0.16036 3.981 0.00120 **
M[, ind]a -0.02912 0.17566 -0.166 0.87054
M[, ind]b -0.21172 0.19665 -1.077 0.29865
M[, ind]c 0.00752 0.18551 0.041 0.96820
M[, ind]d 0.06357 0.18337 0.347 0.73366
Is there some way I can do this so the coefficients have better names ('a'
through 'd')?
As for what I am doing it is not...
l3 <- lm(M[,1]+M[,2]+M[,3]+M[,4]~M[,5])
or
l4 <- lm(M[,1]*M[,2]*M[,3]*M[,4]~M[,5])
cheers
W
[[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, minimal, self-contained, reproducible code.
