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}