Edson,

Combine os resultados dos objetos "an" e "m1" abaixo em um novo data.frame. Ou você pode extrair a matriz do modelo e usa-la novamente dentro de outro lm(), pois a lm() pode receber fórmula (em que cada termo possui seus graus de liberdade reconhecidos) ou pode receber uma matriz (em que o número de colunas-1 da matriz é o grau de liberdade). Veja o CMR.

n <- 100
da <- data.frame(x=runif(n), z=3*runif(n))
da$y <- da$x+da$z+rnorm(n)
m0 <- lm(y~1, da)
m1 <- lm(y~x+z, da)
an <- anova(m0, m1)
str(an)
str(m1)

reg <- model.matrix(m1)
m2 <- lm(y~reg, da)
anova(m2)

À disposição.
Walmes.

==========================================================================
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573
VoIP: (3361 3600) 1053 1173
e-mail: walmes@ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================