test/testthat/testrthpearson.R

In matloff/Rth: Rth: Parallel R through Thrust

describe("rthpearson", {
  it("implements Pearson product moment correlation", {
    x <- c(44.4, 45.9, 41.9, 53.3, 44.7, 44.1, 50.7, 45.2, 60.1)
    y <- c( 2.6,  3.1,  2.5,  5.0,  3.6,  4.0,  5.2,  2.8,  3.8)

    rthval <- rthpearson(x, y)
    rval <- cor(x, y)

    expect_equal(rthval, rval, tolerance = 0.1)
  })
})