tests/testthat/test-dssNlme.R

In IulianD/dsSwissKnifeClient: DataSHIELD Tools and Utilities - client side

test_that("dssNlme_groupedData works", {
  datashield.aggregate(opals["local1"], as.symbol('partialData("Orthodont", NULL, NULL, "nlme")'))
  dssNlme_groupedData(newobj = 'grouped', formula =  distance ~ age | Subject, data ='Orthodont', async = FALSE, datasources = opals["local1"] )
  expect_true(suppressWarnings('grouped' %in% datashield.symbols(opals["local1"])[[1]]))
  lme.model <- dssNlme_lme(fixed = distance ~ age, data = 'grouped', random = ~ 1, async = FALSE, datasources = opals["local1"])
  expect_lte(lme.model$local1$sigma - 1.431592, 1e-03)
})