fix the test for updated output

Author: qinyun-lin

R/test_cop.R

@@ -334,24 +334,25 @@ fig <- ggplot2::ggplot(figTable, ggplot2::aes(x = figTable$ModelLabel)) +
 
   if (to_return == "print") {
     cat(crayon::bold("Coefficient of Proportionality (COP):\n\n"))
-    cat("This function calculates (Oster's) delta* and the exact value of delta.")
+    cat("This function calculates a correlation-based coefficient of\nproportionality (delta) as well as Oster's delta*.")
     cat("\n")
     if (negest == 1) {
       cat("Using the absolute value of the estimated effect, result can be interpreted\nby symmetry.\n")
     }
-    cat(sprintf("Delta* is %.3f (assuming no covariates in the baseline model M1),\nthe exact delta is %.3f, with a bias of %.3f%%.\n", 
+    cat("\n")
+    cat(sprintf("Delta* is %.3f (assuming no covariates in the baseline model M1),\nthe correlation-based delta is %.3f, with a bias of %.3f%%.\n", 
               delta_star, delta_exact, delta_pctbias))
     cat("\n")
-    cat(sprintf("With delta*, the coefficient in the final model will be %.3f.\nWith the exact delta, the coefficient will be %.3f.\n",  
+    cat(sprintf("With delta*, the coefficient in the final model will be %.3f.\nWith the correlation-based delta, the coefficient will be %.3f.\n",  
               eff_x_M3_oster, eff_x_M3))
     cat("\n")
     cat("Use to_return = \"raw_output\" to see more specific results and graphic\npresentation of the result.\n")
     cat("\n")
-    cat("This function calculates a correlation-based based coefficient of\nproportionality (delta) as well as Oster's delta*.\n")
+    cat("This function also calculates conditional RIR that invalidates the statistical inference.\n")
     cat("\n")
-    cat(sprintf("If the replacement cases have a fixed value, then RIR = %.3f.\n", cond_RIR_fixedY))
-    cat(sprintf("If the replacement cases follow a null distribution, then RIR = %.3f.\n", cond_RIR_null))
-    cat(sprintf("If the replacement cases satisfy rxy|Z = 0, then RIR = %.3f.\n", cond_RIR_rxyz))
+    cat(sprintf("If the replacement data points have a fixed value, then RIR = %.3f.\n", cond_RIR_fixedY))
+    cat(sprintf("If the replacement data points follow a null distribution, then RIR = %.3f.\n", cond_RIR_null))
+    cat(sprintf("If the replacement data points satisfy rxy|Z = 0, then RIR = %.3f.\n", cond_RIR_rxyz))
   }
 
 }

R/test_pse.R

@@ -147,7 +147,7 @@ test_pse <- function(est_eff,
         cat(sprintf("The correlation between X and CV is %.3f, and the correlation between\nY and CV is %.3f.\n\n", rxcv, rycv))
         cat(sprintf("Conditional on the covariates, the correlation between X and CV is %.3f,\nand the correlation between Y and CV is %.3f.\n\n", rxcvGz, rycvGz))
         cat(sprintf("Including such CV, the coefficient changes to %.3f, and standard error\nis %.3f.\n\n", eff_x_M3, se_x_M3))
-        cat("Use to_return = \"raw_ouput\" to see more specific results.")
+        cat("Use to_return = \"raw_output\" to see more specific results.")
     }
 
 

tests/testthat/test-helper_output_print.R

@@ -20,15 +20,13 @@ test_that("output_print produces the correct language when index is 'RIR' with e
     # Define expected lines
     expected_lines <- c(
         "Robustness of Inference to Replacement (RIR):",
-        "TO INVALIDATE:",
-        "",
         "RIR = 5",
         "",
         "The estimated effect is 0.6, and specified threshold for inference is 0.5.",
         "To invalidate the inference based on your estimate, 30% of the (0.6)",
         "estimate would have to be due to bias. This implies that to invalidate",
         "the inference one would expect to have to replace 5 (30%) observation",
-        "with cases for which the effect is 0.1 (RIR = 5).",
+        "with data points for which the effect is 0.1 (RIR = 5).",
         "",
         "See Frank et al. (2013) for a description of the method.",
         "",
@@ -64,15 +62,13 @@ test_that("output_print produces the correct language for 'TO SUSTAIN' when inde
     # Define expected lines
     expected_lines <- c(
         "Robustness of Inference to Replacement (RIR):",
-        "TO SUSTAIN:",
-        "",
         "RIR = 10",
         "",
         "The estimated effect is 0.3, and specified threshold for inference is 0.4.",
         "To reach that threshold, 20% of the (0.3) estimate would have to be due",
         "to bias. This implies that to sustain an inference one would expect to have",
-        "to replace 10 (20%) observations with effect of 0.2 with cases with effect",
-        "of 0.5 (RIR = 10).",
+        "to replace 10 (20%) observations with effect of 0.2 with data points with",
+        "effect of 0.5 (RIR = 10).",
         "",
         "See Frank et al. (2013) for a description of the method.",
         "",

tests/testthat/test-test_cop.R

@@ -82,21 +82,25 @@ test_that("test_cop outputs the correct language and values when to_return is 'p
 
     # Define expected lines
     expected_lines <- c(
-        "This function calculates delta* and the exact value of delta.",       
-        "delta* is 0.387 (assuming no covariates in the baseline model M1), ", 
-        "                the exact delta is 1.172, with a bias of -66.980%.",  
-        "With delta*, the coefficient in the final model will be 0.447. ",     
-        "                With the exact delta, the coefficient will be 0.300.",
-        "Use to_return = raw_ouput to see more specific results ",             
-        "        and graphic presentation of the result.",                     
+        "Coefficient of Proportionality (COP):",
+        "",
+        "This function calculates a correlation-based coefficient of",
+        "proportionality (delta) as well as Oster's delta*.",
+        "",
+        "Delta* is 0.387 (assuming no covariates in the baseline model M1),", 
+        "the correlation-based delta is 1.172, with a bias of -66.980%.",  
+        "",
+        "With delta*, the coefficient in the final model will be 0.447.",     
+        "With the correlation-based delta, the coefficient will be 0.300.",
+        "",
+        "Use to_return = \"raw_output\" to see more specific results and graphic",             
+        "presentation of the result.",                     
         "",                                                                    
-        "This function also calculates conditional RIR that ",                 
-        "        invalidates the statistical inference.",                      
-        "If the replacement cases have a fixed value, then RIR = 80.03469 .",  
-        "If the replacement cases follow a null distribution, ",               
-        "        then RIR = 55.31744 .",                                       
-        "If the replacement cases satisfy rxy|Z = 0, then RIR = 56.31727 .",   
-        "" 
+        "This function also calculates conditional RIR that invalidates the statistical inference.",  
+        "",                                                                    
+        "If the replacement data points have a fixed value, then RIR = 80.035.",  
+        "If the replacement data points follow a null distribution, then RIR = 55.317.",                                       
+        "If the replacement data points satisfy rxy|Z = 0, then RIR = 56.317."
     )
 
     # Check each line against expected

tests/testthat/test-test_pse.R

@@ -77,9 +77,20 @@ test_that("test_pse print output is correct", {
     to_return <- "print"
 
     # Expected string
-    expected_output <- paste(
-        "This function calculates the conditions that set the\n            estimated effect approximately equal to the threshold\n            while preserving the standard error.\nThe correlation between X and CV is 0.214,\n                    and the correlation between Y and CV is 0.313.\nConditional on the covariates, the correlation between X\n                    and CV is 0.248, and the correlation between Y and CV\n                    is 0.372.\nIncluding such CV, the coefficient changes to 0.097,\n                    and standard error is 0.054.\nUse to_return = raw_ouput to see more specific results.",
-        sep = ""
+    expected_output <- c(
+        "This function calculates the conditions that set the estimated effect",
+        "approximately equal to the threshold while preserving the standard error.",
+        "",
+        "The correlation between X and CV is 0.214, and the correlation between",
+        "Y and CV is 0.313.",
+        "",
+        "Conditional on the covariates, the correlation between X and CV is 0.248,",
+        "and the correlation between Y and CV is 0.372.",
+        "",
+        "Including such CV, the coefficient changes to 0.097, and standard error",
+        "is 0.054.",
+        "",
+        "Use to_return = \"raw_output\" to see more specific results."
     )
 
     # Capture the output of the print statement
@@ -88,9 +99,9 @@ test_that("test_pse print output is correct", {
     )
 
     # Convert output array into a single string for easier grepl checking
-    output_string <- paste(output, collapse = "\n")
+    # output_string <- paste(output, collapse = "\n")
 
     # Check if the output is as expected
-    expect_equal(output_string, expected_output)
+    expect_equal(output, expected_output)
 })