R/helpers.R

In OHDSI/SkeletonPredictionValidationStudy: Add the title of the package here

# Copyright 2020 Observational Health Data Sciences and Informatics
#
# This file is part of SkeletonPredictionValidationStudy
#
# Licensed under the Apache License, Version 2.0 (the "License");
# you may not use this file except in compliance with the License.
# You may obtain a copy of the License at
#
#     http://www.apache.org/licenses/LICENSE-2.0
#
# Unless required by applicable law or agreed to in writing, software
# distributed under the License is distributed on an "AS IS" BASIS,
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
# See the License for the specific language governing permissions and
# limitations under the License.


addCohortSettings <- function(
  covariateSettings,
  cohortDatabaseSchema,
  cohortTable
){

  if(class(covariateSettings) == 'covariateSettings'){
    covariateSettings <- list(covariateSettings)
  }

  # set the cohort table and database to the settings where the cohorts were generated
  for(i in 1:length(covariateSettings)){

    if('cohortTable' %in% names(covariateSettings[[i]])){
      covariateSettings[[i]]$cohortTable <- cohortTable
    }

    if('cohortDatabaseSchema' %in% names(covariateSettings[[i]])){
      covariateSettings[[i]]$cohortDatabaseSchema <- cohortDatabaseSchema
    }

  }

  return(covariateSettings)
}

createCohorts <- function(
  databaseDetails,
  outputFolder
) {
  if (!file.exists(outputFolder))
    dir.create(outputFolder)

  connection <- DatabaseConnector::connect(databaseDetails$connectionDetails)
  on.exit(DatabaseConnector::disconnect(connection))

  CohortGenerator::createCohortTables(connection = connection,
                                      cohortDatabaseSchema = databaseDetails$cohortDatabaseSchema,
                                      cohortTableNames = CohortGenerator::getCohortTableNames(
                                        cohortTable = databaseDetails$cohortTable
                                      )
  )
  cohortDefinitionSet <- CohortGenerator::getCohortDefinitionSet(packageName = "SkeletonPredictionValidationStudy",
                                                                 settingsFileName = "Cohorts.csv",
                                                                 cohortFileNameValue = "cohortId")
  CohortGenerator::generateCohortSet(connection = connection,
                                     cohortDatabaseSchema = databaseDetails$cohortDatabaseSchema,
                                     cohortTableNames = CohortGenerator::getCohortTableNames(
                                       cohortTable = databaseDetails$cohortTable
                                     ),
                                     cdmDatabaseSchema = databaseDetails$cdmDatabaseSchema,
                                     tempEmulationSchema = databaseDetails$tempEmulationSchema,
                                     cohortDefinitionSet = cohortDefinitionSet)


  # Check number of subjects per cohort:
  message("Counting cohorts")
  counts <- CohortGenerator::getCohortCounts(connection = connection,
                                             cohortDatabaseSchema = databaseDetails$cohortDatabaseSchema,
                                             cohortTable = databaseDetails$cohortTable)

  counts <- addCohortNames(counts)
  utils::write.csv(counts, file.path(outputFolder, "CohortCounts.csv"), row.names = FALSE)
}

addCohortNames <- function(data, IdColumnName = "cohortId", nameColumnName = "cohortName") {
  pathToCsv <- system.file("Cohorts.csv", package = "SkeletonPredictionValidationStudy")

  idToName <- utils::read.csv(pathToCsv)
  idToName <- idToName[order(idToName$cohortId), ]
  idToName <- idToName[!duplicated(idToName$cohortId), ]
  names(idToName)[1] <- IdColumnName
  names(idToName)[2] <- nameColumnName
  data <- merge(data, idToName, all.x = TRUE)
  # Change order of columns:
  idCol <- which(colnames(data) == IdColumnName)
  if (idCol < ncol(data) - 1) {
    data <- data[, c(1:idCol, ncol(data) , (idCol + 1):(ncol(data) - 1))]
  }
  return(data)
}