OK, what is the trick to extracting the overall p value from an lm object?
It shows up in the summary(lm(model)) output but I can't seem to extract it:

> test2 = apply(aa, 1, function(x) summary(lm(x[,1] ~ 0 + x[,3] + x[,6])))
> test2[[1]]

Call:
lm(formula = x[, 1] ~ 0 + x[, 3] + x[, 6])

[omitted summary output]
F-statistic: 40.94 on 2 and 7 DF,  p-value: 0.0001371

It does not seem to be obtainable from anova(lm(model)) either, only the p
values for the individual predictors.
Stumped.

Jim Bouldin
Research Ecologist

	[[alternative HTML version deleted]]

