逐步多元非线性回归

Question

我正在逐步开发回归关系。以下是我的代码和错误消息：

step(nls(q50 ~ aXarea^a1Xslope^a2Xpcp^a3Xforest^a4, 
  start=list(a=6, a1=0.2, a2=1, a3=10, a4=-0.5), direction='backward'))

Error in model.frame.default(formula = ~q50 + area + slope + pcp + forest,  : 
  variable lengths differ (found for '(direction)')

无法找出错误所在。