low dimensional problem call

Author: mcomas

DESCRIPTION

@@ -1,7 +1,7 @@
 Package: coda.count
 Type: Package
 Title: Tools for Compositional Count Data Distributions
-Version: 0.4.0.5
+Version: 0.4.0.6
 Date: 2023-04-19
 Authors@R: c(person("Marc", "Comas-Cufí", role = c("aut", "cre"), 
                     email = "mcomas@imae.udg.edu"))

R/conditional_lrnm_low_dim.R

@@ -29,7 +29,7 @@ fit_one_dimensional_conditional_lrnm = function(X, C = X > 0, V = NULL, V.direct
     stop("Direction V should be specified or V.direction set to 'center' or 'equal'")
   }
   if(is.null(V)){
-    if('equal'){
+    if(V.direction == 'equal'){
       sel = rowSums(X==0)>0
       PART = 2*t(X==0)-1
       Vclr = matrix(0, nrow = ncol(X), ncol = nrow(X))

testing/comparing_methods_02.R

@@ -18,7 +18,7 @@ compare_methods = function(mu,sigma,n,m){
   fit3 = fit_lrnm(X, method = 'laplace')
   fit4 = fit_conditional_lrnm(X, method = 'mc')
   fit5 = fit_conditional_lrnm(X, method = 'hermite')
-  fit6 = fit_conditional_lrnm(X, low.dim = 1)
+  fit6 = fit_one_dimensional_conditional_lrnm(X)
   fit7 = fit_lrnm(X, method = 'vem')
 
   # fit6 = fit_vem_lrnm(X, B)
@@ -42,8 +42,8 @@ lmu = lapply(seq_along(lres), function(i){
   row.names(dres) = NULL
   dres
 })
-dmu = bind_rows(lmu, .id = 'id')
 library(tidyverse)
+dmu = bind_rows(lmu, .id = 'id')
 ggplot(data=dmu) +
   geom_label(aes(x = ilr1, y = ilr2, label = id)) +
   facet_wrap(~method)