Fix a bug in StatStMoE class

- Fix a bug in `StatStMoE$EStep`: when trying to integrate (case `calcE3 = TRUE`), a vector of characters was returned in some cases instead of numerics.
- Update NEWS.md.

NEWS.md

@@ -7,3 +7,8 @@
 * Added a `cran-comments.md` file.
 
 * Added a `NEWS.md` file to track changes to the package.
+
+## Bug fixes
+
+* Fix a bug in `StatStMoE$EStep`: when trying to integrate (case `calcE3 = TRUE`),
+  a vector of characters was returned in some case instead of numerics.

R/StatStMoE.R

@@ -228,7 +228,7 @@ StatStMoE <- setRefClass(
 
         if (calcE3) {
 
-          Integgtx[, k] <- sapply(mik[, k], function(x) try(integrate(f = fx, lower = -Inf, upper = x)$value, silent = TRUE))
+          Integgtx[, k] <- as.numeric(sapply(mik[, k], function(x) try(integrate(f = fx, lower = -Inf, upper = x)$value, silent = TRUE)))
           E3ik[, k] <<- wik[, k] - log((paramStMoE$nu[k] + dik[, k] ^ 2) / 2) - (paramStMoE$nu[k] + 1) / (paramStMoE$nu[k] + dik[, k] ^ 2) + psigamma((paramStMoE$nu[k] + 1) / 2) + ((paramStMoE$lambda[k] * dik[, k] * (dik[, k] ^ 2 - 1)) / sqrt((paramStMoE$nu[k] + 1) * ((paramStMoE$nu[k] + dik[, k] ^ 2) ^ 3))) * dt(mik[, k], paramStMoE$nu[k] + 1) / pt(mik[, k], paramStMoE$nu[k] + 1) + (1 / pt(mik[, k], paramStMoE$nu[k] + 1)) * Integgtx[, k]
         }