Самозагрузка для модели nlme

Question

Я пытаюсь использовать начальную загрузку, чтобы получить ошибки вокруг моей оценки параметров для фиксированных эффектов в следующей модели. Это просто оценка количества раз, которое животное может пересечь дорогу, основываясь на расстоянии дороги от ручья.

Когда я запускаю модель nlme, она сходится и все хорошо. Я пробовал несколько разных способов сделать загрузку, но безуспешно. Я пробовал использовать как загрузочный пакет, так и просто разработать немного кода, который пересчитывает мои данные и переносит оценки параметров в новые векторы.

Ниже приведена моя попытка позднее и полученные сообщения об ошибках. Любая помощь будет принята с благодарностью.

bv <- numeric(100)
cv <- numeric(100)
for(i in 1:100){
  ss <- sample(1:130, replace=T)
  y <- nwfcross[ss]
  x <- nwfdist[ss]
  modelb <- nlme(y~a*exp(-b*x), fixed=a+b~1, 
         random=a+b~1|nwfid, start=c(a=300,b=0.016))
  bv[i] <- coef(modelb)[1]
  cv[i] <- coef(modelb)[2]
}

Error in nlme.formula(y ~ a * exp(-b * x), fixed = a + b ~ 1, random = a +  : 
  Maximum number of iterations reached without convergence
In addition: There were 50 or more warnings (use warnings() to see the first 50)

Предупреждающие сообщения:

1: Singular precision matrix in level -1, block 1
2: Singular precision matrix in level -1, block 1
3: Singular precision matrix in level -1, block 1….