Skip to content

CMIM output differences between v6 and v7

Hi,

the following data produces different outputs between v6 and v7. Which one operates correctly?

v6

library(praznik)

X = structure(list(a = c(1, 2, 5.3, 6, -2, 4, 8.3, 9.2, 10.1, 1, 
                         2, 5.3, 6, -2, 4, 8.3, 9.2, 10.1), b = structure(c(1L, 3L, 2L, 
                                                                            1L, 3L, 2L, 1L, 3L, 2L, 1L, 3L, 2L, 1L, 3L, 2L, 1L, 3L, 2L), class = "factor", .Label = c("one", 
                                                                                                                                                                      "three", "two")), c = structure(c(1L, 1L, 1L, 2L, 1L, 2L, 1L, 
                                                                                                                                                                                                        1L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 2L), .Label = c("TRUE", 
                                                                                                                                                                                                                                                                "FALSE"), class = "factor"), d = c(1L, 3L, 5L, 7L, 9L, 17L, 1L, 
                                                                                                                                                                                                                                                                                                   3L, 5L, 7L, 9L, 17L, 1L, 3L, 5L, 7L, 9L, 17L), const1 = structure(c(1L, 
                                                                                                                                                                                                                                                                                                                                                                       3L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 
                                                                                                                                                                                                                                                                                                                                                                       2L), .Label = c("c1", "c2"), class = "factor"), const2 = c(1, 
                                                                                                                                                                                                                                                                                                                                                                                                                                  2, 5.3, 6, -2, 4, 8.3, 9.2, 10.1, 1, 2, 5.3, 6, -2, 4, 8.3, 9.2, 
                                                                                                                                                                                                                                                                                                                                                                                                                                  10.1)), row.names = c(NA, -18L), class = "data.frame")
Y = structure(c(1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 
                2L, 1L, 2L, 1L, 2L), .Label = c("c1", "c2"), class = "factor")

CMIM(X, Y)
#> $selection
#> d c 
#> 4 3 
#> 
#> $score
#>            d            c 
#> 6.931472e-01 1.110223e-16

Created on 2019-10-07 by the reprex package (v0.3.0)

Session info
devtools::session_info()
#> ─ Session info ──────────────────────────────────────────────────────────
#>  setting  value                       
#>  version  R version 3.6.1 (2019-07-05)
#>  os       Arch Linux                  
#>  system   x86_64, linux-gnu           
#>  ui       X11                         
#>  language (EN)                        
#>  collate  en_US.UTF-8                 
#>  ctype    en_US.UTF-8                 
#>  tz       Europe/Berlin               
#>  date     2019-10-07                  
#> 
#> ─ Packages ──────────────────────────────────────────────────────────────
#>  package     * version    date       lib source                        
#>  assertthat    0.2.1      2019-03-21 [1] CRAN (R 3.6.0)                
#>  backports     1.1.5      2019-10-02 [1] CRAN (R 3.6.1)                
#>  callr         3.3.2      2019-09-22 [1] CRAN (R 3.6.1)                
#>  cli           1.1.0      2019-03-19 [1] CRAN (R 3.6.0)                
#>  crayon        1.3.4      2017-09-16 [1] CRAN (R 3.6.0)                
#>  desc          1.2.0      2018-05-01 [1] CRAN (R 3.6.0)                
#>  devtools      2.2.1      2019-09-24 [1] CRAN (R 3.6.1)                
#>  digest        0.6.21     2019-09-20 [1] CRAN (R 3.6.1)                
#>  ellipsis      0.3.0      2019-09-20 [1] CRAN (R 3.6.1)                
#>  evaluate      0.14       2019-05-28 [1] CRAN (R 3.6.0)                
#>  fs            1.3.1      2019-05-06 [1] CRAN (R 3.6.0)                
#>  glue          1.3.1      2019-03-12 [1] CRAN (R 3.6.0)                
#>  highr         0.8        2019-03-20 [1] CRAN (R 3.6.0)                
#>  htmltools     0.4.0      2019-10-04 [1] CRAN (R 3.6.1)                
#>  knitr         1.25       2019-09-18 [1] CRAN (R 3.6.1)                
#>  magrittr      1.5        2014-11-22 [1] CRAN (R 3.6.0)                
#>  memoise       1.1.0      2017-04-21 [1] CRAN (R 3.6.0)                
#>  pkgbuild      1.0.5      2019-08-26 [1] CRAN (R 3.6.1)                
#>  pkgload       1.0.2      2018-10-29 [1] CRAN (R 3.6.0)                
#>  praznik     * 6.0.0      2019-06-21 [1] CRAN (R 3.6.1)                
#>  prettyunits   1.0.2      2015-07-13 [1] CRAN (R 3.6.0)                
#>  processx      3.4.1      2019-07-18 [1] CRAN (R 3.6.1)                
#>  ps            1.3.0      2018-12-21 [1] CRAN (R 3.6.0)                
#>  R6            2.4.0      2019-02-14 [1] CRAN (R 3.6.0)                
#>  Rcpp          1.0.2      2019-07-25 [1] CRAN (R 3.6.1)                
#>  remotes       2.1.0      2019-06-24 [1] CRAN (R 3.6.0)                
#>  rlang         0.4.0      2019-06-25 [1] CRAN (R 3.6.1)                
#>  rmarkdown     1.16       2019-10-01 [1] CRAN (R 3.6.1)                
#>  rprojroot     1.3-2      2018-01-03 [1] CRAN (R 3.6.0)                
#>  sessioninfo   1.1.1      2018-11-05 [1] CRAN (R 3.6.0)                
#>  stringi       1.4.3      2019-03-12 [1] CRAN (R 3.6.0)                
#>  stringr       1.4.0      2019-02-10 [1] CRAN (R 3.6.0)                
#>  testthat      2.2.1      2019-07-25 [1] CRAN (R 3.6.1)                
#>  usethis       1.5.1.9000 2019-10-07 [1] Github (r-lib/usethis@3015465)
#>  withr         2.1.2      2018-03-15 [1] CRAN (R 3.6.0)                
#>  xfun          0.10       2019-10-01 [1] CRAN (R 3.6.1)                
#>  yaml          2.2.0      2018-07-25 [1] CRAN (R 3.6.0)                
#> 
#> [1] /home/pjs/R/x86_64-pc-linux-gnu-library/3.6
#> [2] /usr/lib/R/library

v7

library(praznik)

X = structure(list(a = c(1, 2, 5.3, 6, -2, 4, 8.3, 9.2, 10.1, 1, 
                         2, 5.3, 6, -2, 4, 8.3, 9.2, 10.1), b = structure(c(1L, 3L, 2L, 
                                                                            1L, 3L, 2L, 1L, 3L, 2L, 1L, 3L, 2L, 1L, 3L, 2L, 1L, 3L, 2L), class = "factor", .Label = c("one", 
                                                                                                                                                                      "three", "two")), c = structure(c(1L, 1L, 1L, 2L, 1L, 2L, 1L, 
                                                                                                                                                                                                        1L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 2L), .Label = c("TRUE", 
                                                                                                                                                                                                                                                                "FALSE"), class = "factor"), d = c(1L, 3L, 5L, 7L, 9L, 17L, 1L, 
                                                                                                                                                                                                                                                                                                   3L, 5L, 7L, 9L, 17L, 1L, 3L, 5L, 7L, 9L, 17L), const1 = structure(c(1L, 
                                                                                                                                                                                                                                                                                                                                                                       2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 
                                                                                                                                                                                                                                                                                                                                                                       2L), .Label = c("c1", "c2"), class = "factor"), const2 = c(1, 
                                                                                                                                                                                                                                                                                                                                                                                                                                  2, 5.3, 6, -2, 4, 8.3, 9.2, 10.1, 1, 2, 5.3, 6, -2, 4, 8.3, 9.2, 
                                                                                                                                                                                                                                                                                                                                                                                                                                  10.1)), row.names = c(NA, -18L), class = "data.frame")
Y = structure(c(1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 
                2L, 1L, 2L, 1L, 2L), .Label = c("c1", "c2"), class = "factor")

CMIM(X, Y)
#> $selection
#> d 
#> 4 
#> 
#> $score
#>         d 
#> 0.6931472

Created on 2019-10-07 by the reprex package (v0.3.0)

Session info
devtools::session_info()
#> ─ Session info ──────────────────────────────────────────────────────────
#>  setting  value                       
#>  version  R version 3.6.1 (2019-07-05)
#>  os       Arch Linux                  
#>  system   x86_64, linux-gnu           
#>  ui       X11                         
#>  language (EN)                        
#>  collate  en_US.UTF-8                 
#>  ctype    en_US.UTF-8                 
#>  tz       Europe/Berlin               
#>  date     2019-10-07                  
#> 
#> ─ Packages ──────────────────────────────────────────────────────────────
#>  package     * version    date       lib source                        
#>  assertthat    0.2.1      2019-03-21 [1] CRAN (R 3.6.0)                
#>  backports     1.1.5      2019-10-02 [1] CRAN (R 3.6.1)                
#>  callr         3.3.2      2019-09-22 [1] CRAN (R 3.6.1)                
#>  cli           1.1.0      2019-03-19 [1] CRAN (R 3.6.0)                
#>  crayon        1.3.4      2017-09-16 [1] CRAN (R 3.6.0)                
#>  desc          1.2.0      2018-05-01 [1] CRAN (R 3.6.0)                
#>  devtools      2.2.1      2019-09-24 [1] CRAN (R 3.6.1)                
#>  digest        0.6.21     2019-09-20 [1] CRAN (R 3.6.1)                
#>  ellipsis      0.3.0      2019-09-20 [1] CRAN (R 3.6.1)                
#>  evaluate      0.14       2019-05-28 [1] CRAN (R 3.6.0)                
#>  fs            1.3.1      2019-05-06 [1] CRAN (R 3.6.0)                
#>  glue          1.3.1      2019-03-12 [1] CRAN (R 3.6.0)                
#>  highr         0.8        2019-03-20 [1] CRAN (R 3.6.0)                
#>  htmltools     0.4.0      2019-10-04 [1] CRAN (R 3.6.1)                
#>  knitr         1.25       2019-09-18 [1] CRAN (R 3.6.1)                
#>  magrittr      1.5        2014-11-22 [1] CRAN (R 3.6.0)                
#>  memoise       1.1.0      2017-04-21 [1] CRAN (R 3.6.0)                
#>  pkgbuild      1.0.5      2019-08-26 [1] CRAN (R 3.6.1)                
#>  pkgload       1.0.2      2018-10-29 [1] CRAN (R 3.6.0)                
#>  praznik     * 7.0.0      2019-10-01 [1] CRAN (R 3.6.1)                
#>  prettyunits   1.0.2      2015-07-13 [1] CRAN (R 3.6.0)                
#>  processx      3.4.1      2019-07-18 [1] CRAN (R 3.6.1)                
#>  ps            1.3.0      2018-12-21 [1] CRAN (R 3.6.0)                
#>  R6            2.4.0      2019-02-14 [1] CRAN (R 3.6.0)                
#>  Rcpp          1.0.2      2019-07-25 [1] CRAN (R 3.6.1)                
#>  remotes       2.1.0      2019-06-24 [1] CRAN (R 3.6.0)                
#>  rlang         0.4.0      2019-06-25 [1] CRAN (R 3.6.1)                
#>  rmarkdown     1.16       2019-10-01 [1] CRAN (R 3.6.1)                
#>  rprojroot     1.3-2      2018-01-03 [1] CRAN (R 3.6.0)                
#>  sessioninfo   1.1.1      2018-11-05 [1] CRAN (R 3.6.0)                
#>  stringi       1.4.3      2019-03-12 [1] CRAN (R 3.6.0)                
#>  stringr       1.4.0      2019-02-10 [1] CRAN (R 3.6.0)                
#>  testthat      2.2.1      2019-07-25 [1] CRAN (R 3.6.1)                
#>  usethis       1.5.1.9000 2019-10-07 [1] Github (r-lib/usethis@3015465)
#>  withr         2.1.2      2018-03-15 [1] CRAN (R 3.6.0)                
#>  xfun          0.10       2019-10-01 [1] CRAN (R 3.6.1)                
#>  yaml          2.2.0      2018-07-25 [1] CRAN (R 3.6.0)                
#> 
#> [1] /home/pjs/R/x86_64-pc-linux-gnu-library/3.6
#> [2] /usr/lib/R/library