tests/testthat/test_R2.R

In stephenslab/flashr2: Empirical Bayes Matrix Factorization

context("test_R2")

test_that("test R2 computation matches R^2 when EL2=EL^2 and EF2=EF^2 ", {
  set.seed(1)
  Y = matrix(nrow=5,ncol=20,rnorm(100))
  data = flash_set_data(Y)

  f = flash_init_fn(data,"udv_svd",1)
  R = Y - flash_get_fitted_values(f)
  R2 = R^2
  expect_equal(R2,flash_get_R2(data,f))
})