diff --git a/R/obj.gauss.R b/R/obj.gauss.R index 47e1d60..43c393a 100644 --- a/R/obj.gauss.R +++ b/R/obj.gauss.R @@ -18,8 +18,8 @@ #' @author Francois Pelletier obj.gauss <- function(param,Y,meanf,variancef,dmean,dsd,Q=diag(4)) { - eqn.gauss(param,Y,meanf,variancef,dmean,dsd) %*% Q %*% - eqn.gauss(param,Y,meanf,variancef,dmean,dsd) + as.vector(eqn.gauss(param,Y,meanf,variancef,dmean,dsd) %*% Q %*% + eqn.gauss(param,Y,meanf,variancef,dmean,dsd)) }