>>>>> "Manuel" == Manuel Morales <manuel.a.mora...@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 ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.