add functionality to readdesign

R/readdesign.R

@@ -20,7 +20,25 @@
 
 
 
-readdesign <- function(design = designfile, designtype = destype) {
+readdesign <- function(design = designfile, designtype = NULL) {
+
+  # if (is.null(designtype)) {
+  #
+  #
+  #
+  #   # Check if the string ends with ".ngd"
+  #   if (grepl("\\.ngd$", designfile)) {
+  #     # Code to execute if condition is true
+  #     designtype = "ngene"
+  #     print("I guessed it is an ngene file")
+  #   } else {
+  #     # Code to execute if condition is false
+  #     designtype = "spdesign"
+  #     print("I assume it is a spdesign")
+  #   }
+  #
+  # }
+
   design <- switch(designtype,
                    "ngene" = suppressWarnings(readr::read_delim(design,
                                                                 delim = "\t",

tests/manual-tests/test_decisiongroups.R

@@ -256,9 +256,9 @@ apollo_probabilities=function(apollo_beta, apollo_inputs, functionality="estimat
 # ################################################################# #
 
 ### Optional starting values search
-# apollo_beta=apollo_searchStart(apollo_beta, apollo_fixed,apollo_probabilities, apollo_inputs)
+apollo_beta=apollo_searchStart(apollo_beta, apollo_fixed,apollo_probabilities, apollo_inputs)
 
-model = apollo_estimate(apollo_beta, apollo_fixed, apollo_probabilities, apollo_inputs ,estimate_settings =list(estimationRoutine="bhhh" , hessianRoutine= "maxLik"))
+model = apollo_estimate(apollo_beta, apollo_fixed, apollo_probabilities, apollo_inputs ) #,estimate_settings =list(estimationRoutine="bhhh" , hessianRoutine= "maxLik")
 
 # ################################################################# #
 #### MODEL OUTPUTS                                               ####