Error in q6[i] <- RHfit$ax[26] + RHfit$bx[26] * q5[i] + q55[i] :
  更换参数长度为零