Via nls que não altera as estimativas dos parâmetros você pode fazer:
mg=nls(r ~ a[t] + b[t]*x,data=d,start=list(a=c(1,1),b=c(1,1)))
ma=nls(r ~ a + b[t]*x,data=d,start=list(a=1,b=c(1,1)))
anova(mg,ma)
mb=nls(r ~ a[t] + b*x,data=d,start=list(a=c(1,1),b=1))
anova(mg,mb)
(S,f,P)
Allaman
\begin{signature}
<<>>=
Prof. Dr. Ivan Bezerra Allaman
Universidade Estadual de Santa Cruz
Departamento de Ciências Exatas e Tecnológicas
Ilhéus/BA - Brasil
Fone: +55 73 3680-5076
E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com
@
\end{signature}