>>>>> "Manuel" == Manuel Morales <Manuel.A.Morales at williams.edu> writes: > nls(y ~ a[fac]*x^b, start=list(a=c(1,1), b=0.25)) Did you mean a[f]? nls(y ~ a[f]*x^b, start=list(a=c(1,1), b=0.25)) Mike