With an additional parameter (i.e., the variational parameter for Vi), the estimating equations proposed in [24] are(4)(5)where we employ the usual statistical convention of lower Roman variables, xij, to denote the observed version of the random data, Xij.