R Under development (unstable) (2023-07-19 r84711 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2023 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> 
> source("helpers.R")
> 
> library(reclin2)
Loading required package: data.table

Attaching package: 'reclin2'

The following object is masked from 'package:base':

    identical

> 
> x <- data.table(a = c(1,1,2,2), b = c(1,2,1,2))
> y <- data.table(a = c(3,3,2,2), b = c(1,2,1,2))
> 
> pairs <- pair_minsim(x, y, on = c("b", "a"), minsim = 2)
> expect_equal(names(pairs), c(".x", ".y", "simsum"))
> expect_equal(pairs$.x, c(3,4))
> expect_equal(pairs$.y, c(3,4))
> expect_equal(attr(pairs, "x"), x)
> expect_equal(attr(pairs, "y"), y)
> 
> pairs <- pair_minsim(x, y, on = c("a", "b"), minsim = 2, 
+   keep_simsum = FALSE, add_xy = FALSE)
> expect_equal(names(pairs), c(".x", ".y"))
> expect_equal(attr(pairs, "x"), NULL)
> expect_equal(attr(pairs, "y"), NULL)
> 
> pairs <- pair_minsim(x, y, on = c("b", "a"), minsim = 3)
> expect_equal(names(pairs), c(".x", ".y", "simsum"))
> expect_equal(nrow(pairs), 0)
> expect_equal(attr(pairs, "x"), x)
> expect_equal(attr(pairs, "y"), y)
> 
> pairs <- pair_minsim(x, on = c("b", "a"), minsim = 1,
+   deduplication = TRUE)
> expect_equal(names(pairs), c(".x", ".y", "simsum"))
> expect_equal(pairs$.x, c(1,1,2,3))
> expect_equal(pairs$.y, c(2,3,4,4))
> 
> # Regression test for issue #9: pairs are not generated when there are missing
> # values
> x <- data.table(a = c(1,1,2,2), b = c(1,2,1,2))
> y <- data.table(a = c(3,3,2,2), b = c(1,2,1,2))
> x$a[1] <- NA
> pairs <- pair_minsim(x, y, on = c("b", "a"), minsim = 1)
> expect_equal(pairs$.x, c(1L, 1L, 2L, 2L, 3L, 3L, 3L, 4L, 4L, 4L))
> expect_equal(pairs$.y, c(1L, 3L, 2L, 4L, 1L, 3L, 4L, 2L, 3L, 4L))
> expect_equal(pairs$simsum, c(1, 1, 1, 1, 1, 2, 1, 1, 1, 2))
> 
> 
> # ===== ON_BLOCKING
> 
> x <- data.table(a = c(1,1,2,2), b = c(1,2,1,2))
> y <- data.table(a = c(3,3,2,2), b = c(1,2,1,2))
> 
> pairs <- pair_minsim(x, y, on = "a", on_blocking = "b")
> expect_equal(names(pairs), c(".x", ".y", "simsum"))
> expect_equal(pairs$.x, c(1,1,3,3,2,2,4,4))
> expect_equal(pairs$.y, c(1,3,1,3,2,4,2,4))
> expect_equal(pairs$simsum, c(0,0,0,1,0,0,0,1))
> expect_equal(attr(pairs, "x"), x)
> expect_equal(attr(pairs, "y"), y)
> 
> # Unsing non-existent columns
> expect_error(pair_minsim(x, y, on = "a", on_blocking = "foo"))
> expect_error(pair_minsim(x, y, on = c("a", "foo"), on_blocking = "b"))
> 
> pairs <- pair_minsim(x, y, on = "a", on_blocking = "b", minsim = 1)
> expect_equal(names(pairs), c(".x", ".y", "simsum"))
> expect_equal(pairs$.x, c(3,4))
> expect_equal(pairs$.y, c(3,4))
> expect_equal(pairs$simsum, c(1,1))
> expect_equal(attr(pairs, "x"), x)
> expect_equal(attr(pairs, "y"), y)
> 
> pairs <- pair_minsim(x, y, on = "a", on_blocking = "b", minsim = 2)
> expect_equal(names(pairs), c(".x", ".y", "simsum"))
> expect_equal(nrow(pairs), 0)
> expect_equal(attr(pairs, "x"), x)
> expect_equal(attr(pairs, "y"), y)
> 
> # Missing values
> x <- data.table(a = c(1,NA,2,2), b = c(1,2,1,2))
> y <- data.table(a = c(3,3,2,2), b = c(1,2,NA,2))
> pairs <- pair_minsim(x, y, on = "a", on_blocking = "b", minsim = 1)
> expect_equal(names(pairs), c(".x", ".y", "simsum"))
> expect_equal(pairs$.x, c(4))
> expect_equal(pairs$.y, c(4))
> expect_equal(pairs$simsum, c(1))
> expect_equal(attr(pairs, "x"), x)
> expect_equal(attr(pairs, "y"), y)
> 
> 
> proc.time()
   user  system elapsed 
   0.37    0.06    0.42