library( systemfit )
hg.formula<-hg~x*x*x-y-1
dg.formula<-dg~y-7
start.values<-c(x=2,y=7)
model <- list( hg.formula, dg.formula )
ppine<-list(hg=0,dg=0)
model.ols <- nlsystemfit( "OLS", model, start.values,data=ppine )  执行后显示：

Error in diag(diag(qr.solve(S, tol = solvtol))) :
  在为'diag'函数选择方法时评估'x'参数出了错: Error in diag(qr.solve(S, tol = solvtol)) :
  在为'diag'函数选择方法时评估'x'参数出了错: Error in qr.default(a, tol = tol) : 外接函数调用时不能有NA/Na

请指点