关于R语言多水平线性回归分析

我再使用R语言分析嵌套模型的时候，遇到了一下问题，希望可以有大神指导一下。
数据是学校，班级，学生组成的嵌套模型。在模型构建的开始，我先确定模型中应该使用几水平。首先我先把三个水平均作为随机项放进模型中
model2.1<-lmer(disease~
age+gender+height+weight+
(1|school)+(1|class)+(1|stu),
data=mydata,REML = FALSE,
control = lmerControl(
optimizer = “optimx”, optCtrl=list(method=‘nlminb’)))
summary(model2.1)
但是不幸的是，我发现，模型显示了警告：
boundary (singular) fit: see ?isSingular
Warning messages:
1: In optwrap(optimizer, devfun, getStart(start, rho$pp),lower=rho$lower, :
convergence code 1 from optimx
2: Model failed to converge with 2 negative eigenvalues: -5.5e-05 -2.2e-03
于是我把模型改成两水平，仅仅包含学校和班级两个随机项，但是仍然警告：
Warning messages:
1: In optwrap(optimizer, devfun, getStart(start, rho$pp),lower=rho$lower, :
convergence code 1 from optimx
2: In checkConv(attr(opt, “derivs”), opt$par,ctrl=control$checkConv, :
unable to evaluate scaled gradient
3: In checkConv(attr(opt, “derivs”), opt$par,ctrl=control$checkConv, :
Model failed to converge: degenerate Hessian with 1 negative eigenvalues
4: Model failed to converge with 1 negative eigenvalue: -5.1e-04
最终我把模型改成仅包含学校这一个随机项的，结果没有报错，但是AIC还是蛮大的，5000左右。
请问有人知道这是什么原因吗？谢谢指点~