compute_median() change

R/compute_functions.R

@@ -244,14 +244,15 @@ compute_median <- function(file, file_name, variable, weight = NULL, na.rm = FAL
   }
 
 
-  # the original function from Hsmisc is incredibly buggy. It would return a result with NAs and na.rm = F.
+  # the original function from Hsmisc is incredibly buggy. It would return a result with NAs and na.rm = F.  ## OLD COMMENT 
   if(!na.rm & (any(is.na(file[[variable]])) | any(is.na(weight_)) )){
     return(NA)
   }
 
-  unname(matrixStats::weightedMedian(x = file[[variable]],
-                      w = weight_,
-                      na.rm = na.rm))
+  unname(wtd.quantile(x = var, weights = weight_, probs = 0.5,
+                      na.rm = na.rm, normwt = TRUE))  
+  # previous functions matrixStats::weightedMedian wasn't outputing the same values as STATA, when weights were being applied. 
+  # The argument "ties" was also working is a very wierd way. We will resort to the same method as in compute_percentiles() but with probs = 0.5. 
 
 }