[R] nlme

[Santiago Bueno]

Dear people:


I am getting the following error when trying to run the piece of code
below. Any insights??



Error in nlme.formula(Btronc ~ a + b * dbh^2 * haut, data = cbind(dat,  :

  object 'a' not found



library(nlme)

start <- coef(lm(Btronc~I(dbh**2*haut),data=dat))

names(start) <- c("a","b")

summary(nlme(Btronc~a+b*dbh**2*haut, data=cbind(dat,g="a"), fixed=a+b_1,
start=start,

groups=~g, weights=varPower(form=~dbh)))


Best,


Santiago

	[[alternative HTML version deleted]]