Browse code

better tests

deril2605 authored on 18/02/2024 17:17:58
Showing 3 changed files

... ...
@@ -187,7 +187,6 @@ qcServer <- function(input, output, session,parent_session, loadpage_input,get_d
187 187
                             address = file
188 188
         )
189 189
 
190
-        # return(plot)
191 190
         
192 191
       } else if (loadpage_input()$BIO == "PTM"){
193 192
         
... ...
@@ -196,10 +195,6 @@ qcServer <- function(input, output, session,parent_session, loadpage_input,get_d
196 195
                             which.PTM = protein,
197 196
                             summaryPlot = input$summ,
198 197
                             address = file)
199
-        # print("outsidedataProcessPlotsPTM ")
200
-        # print(length(plotly_output))
201
-        # # return which plot here, first?
202
-        # return(plotly_output[[1]])
203 198
         
204 199
       } else{
205 200
         plot <- dataProcessPlots(data = preprocess_data(),
... ...
@@ -190,7 +190,6 @@ getAnnot1 <- function(input) {
190 190
 }
191 191
 
192 192
 getData <- function(input) {
193
-  print("+++++++++ In getDataDeril +++++++++")
194 193
   show_modal_spinner()
195 194
   ev_maxq = getEvidence(input)
196 195
   pg_maxq = getProteinGroups(input)
... ...
@@ -208,14 +207,11 @@ getData <- function(input) {
208 207
   cat(file=stderr(), "Reached in get_data\n")
209 208
 
210 209
   cat(file=stderr(), paste("File type is",input$filetype,"\n"))
211
-  cat(file=stderr(), "herrr\n")
212 210
   if(is.null(input$filetype)) {
213 211
     return(NULL)
214 212
   }
215
-  cat(file=stderr(), "herrr\n")
216 213
   if(input$filetype == 'sample') {
217 214
     if(input$BIO != "PTM" && input$DDA_DIA =='LType' && input$LabelFreeType == "SRM_PRM") {
218
-      cat(file=stderr(), "herrr\n")
219 215
       mydata = MSstats::DDARawData
220 216
     }
221 217
     else if(input$BIO != "PTM" &&  input$DDA_DIA == 'LType' && input$LabelFreeType == "DDA") {
... ...
@@ -344,8 +340,6 @@ getData <- function(input) {
344 340
                                              mod_id=input$mod_id_spec)
345 341
 
346 342
     } else if (input$filetype=='sky') {
347
-      print(":checll")
348
-      cat(file=stderr(), "checll\n")
349 343
       mydata = read_excel(input$ptm_input$datapath)
350 344
       mydata_protein = try(read_excel(input$ptm_protein_input$datapath),silent=TRUE)
351 345
       annotation = try(read.csv(input$ptm_annot$datapath),silent=TRUE)
... ...
@@ -855,10 +855,10 @@ test_that("get summary 1 PTM PTMTMT:Yes", {
855 855
     mock_input$filetype = "sample"
856 856
     mock_input$BIO <- "PTM"
857 857
     mock_input$DDA_DIA <- "TMT"
858
-    # stub(getSummary1,"getData",mockGetData(mock_input))
859
-    
860
-    output <- getSummary1(mock_input,MSstatsPTM::raw.input.tmt)
858
+    stub(getSummary1,"getData",mockGetData(mock_input))
861 859
     
860
+    output <- getSummary1(mock_input,getData(mock_input))
861
+
862 862
     expected_names <- c("Number of Conditions","Number of PTM Mixtures","Number of PTM Biological Replicates","Number of PTM MS runs","Number of PTM Technical Replicates","Number of Unmod Mixtures","Number of Unmod Biological Replicates","Number of Unmod MS runs","Number of Unmod Technical Replicates")
863 863
     expect_type(output,"list")
864 864
     expect_identical(rownames(output), expected_names)
... ...
@@ -870,9 +870,9 @@ test_that("get summary 1 PTM PTMTMT:No", {
870 870
     mock_input$filetype = "sample"
871 871
     mock_input$BIO <- "PTM"
872 872
     mock_input$DDA_DIA <- "LType"
873
-    # stub(getSummary1,"getData",mockGetData(mock_input))
873
+    stub(getSummary1,"getData",mockGetData(mock_input))
874 874
     
875
-    output <- getSummary1(mock_input,MSstatsPTM::raw.input)
875
+    output <- getSummary1(mock_input,getData(mock_input))
876 876
     
877 877
     expected_names <- c("Number of Conditions","Number of PTM Biological Replicates","Number of PTM MS runs","Number of Unmod Conditions","Number of Unmod Biological Replicates","Number of Unmod MS runs")
878 878
     expect_type(output,"list")
... ...
@@ -886,9 +886,9 @@ test_that("get summary 1 Other:DDA", {
886 886
     mock_input$filetype = "sample"
887 887
     mock_input$DDA_DIA <- "LType"
888 888
     mock_input$LabelFreeType <- "DDA"
889
-    # stub(getSummary1,"getData",mockGetData(mock_input))
889
+    stub(getSummary1,"getData",mockGetData(mock_input))
890 890
     
891
-    output <- getSummary1(mock_input,MSstats::DDARawData)
891
+    output <- getSummary1(mock_input,getData(mock_input))
892 892
     
893 893
     expected_names <- c("Number of Conditions","Number of Biological Replicates","Number of Technical Replicates","Number of Fractions","Number of MS runs")
894 894
     expect_type(output,"list")
... ...
@@ -919,9 +919,9 @@ test_that("get summary 2 PTM PTMTMT:Yes", {
919 919
     mock_input$filetype = "sample"
920 920
     mock_input$BIO <- "PTM"
921 921
     mock_input$DDA_DIA <- "TMT"
922
-    # stub(getSummary2,"getData",mockGetData(mock_input))
922
+    stub(getSummary2,"getData",mockGetData(mock_input))
923 923
     
924
-    output <- getSummary2(mock_input,MSstatsPTM::raw.input.tmt)
924
+    output <- getSummary2(mock_input,getData(mock_input))
925 925
     print(rownames(output))
926 926
     
927 927
     expected_names <- c("Number of PTMs","Number of PTM Features",
... ...
@@ -940,9 +940,9 @@ test_that("get summary 2 PTM PTMTMT:No", {
940 940
     mock_input$filetype = "sample"
941 941
     mock_input$BIO <- "PTM"
942 942
     mock_input$DDA_DIA <- "LType"
943
-    # stub(getSummary2,"getData",mockGetData(mock_input))
943
+    stub(getSummary2,"getData",mockGetData(mock_input))
944 944
     
945
-    output <- getSummary2(mock_input,MSstatsPTM::raw.input)
945
+    output <- getSummary2(mock_input,getData(mock_input))
946 946
     print(rownames(output))
947 947
     
948 948
     expected_names <- c("Number of PTMs","Number of PTM Features","Number of Features/PTM","PTM Intensity Range","Number of Unmod Proteins","Number of Protein Peptides","Number of Protein Features","Number of Features/Peptide","Number of Peptides/Protein","Protein Intensity Range") 
... ...
@@ -957,9 +957,9 @@ test_that("get summary 2 Other:DDA", {
957 957
     mock_input$BIO <- "Protein"
958 958
     mock_input$DDA_DIA <- "LType"
959 959
     mock_input$LabelFreeType <- "DDA"
960
-    # stub(getSummary2,"getData",mockGetData(mock_input))
960
+    stub(getSummary2,"getData",mockGetData(mock_input))
961 961
     
962
-    output <- getSummary2(mock_input,MSstats::DDARawData)
962
+    output <- getSummary2(mock_input,getData(mock_input))
963 963
     expected_names <- c("Number of Proteins","Number of Peptides","Number of Features","Number of Peptides/Protein","Number of Features/Peptide","Intensity Range") 
964 964
     expect_type(output,"list")
965 965
     expect_identical(rownames(output), expected_names)
... ...
@@ -992,7 +992,7 @@ test_that("preprocessData QC, PTM and PTMTMT: No", {
992 992
     mock_input$norm = "equalizeMedians"
993 993
     mock_input$log = "2"
994 994
 
995
-    # stub(preprocessData,"getData",mockGetData(mock_input))
995
+    stub(preprocessData,"getData",mockGetData(mock_input))
996 996
     stub(preprocessData,"loadpage_input",mock_input)
997 997
     stub(preprocessData,"qc_input",mock_input,2)
998 998
 
... ...
@@ -1006,14 +1006,13 @@ test_that("preprocessData QC, PTM and PTMTMT: No", {
1006 1006
         mockery::stub(MSstatsShiny::lf_summarization_loop, "remove_modal_progress", NULL);
1007 1007
         MSstatsShiny::lf_summarization_loop(...)})
1008 1008
 
1009
-    output <- preprocessData(mock_input,mock_input,MSstatsPTM::raw.input)
1009
+    output <- preprocessData(mock_input,mock_input,getData(mock_input))
1010 1010
     expected_names <- c("PTM","PROTEIN")
1011 1011
     expect_type(output,"list")
1012 1012
     expect_identical(names(output), expected_names)
1013 1013
   })
1014 1014
 })
1015 1015
  
1016
-# err
1017 1016
 test_that("preprocessData QC, PTM and PTMTMT: Yes", {
1018 1017
   suppressWarnings({
1019 1018
     mock_input$BIO <- "PTM"
... ...
@@ -1023,7 +1022,7 @@ test_that("preprocessData QC, PTM and PTMTMT: Yes", {
1023 1022
     mock_input$log = "2"
1024 1023
     mock_input$summarization = "msstats"
1025 1024
 
1026
-    # stub(preprocessData,"getData",mockGetData(mock_input))
1025
+    stub(preprocessData,"getData",mockGetData(mock_input))
1027 1026
     stub(preprocessData,"loadpage_input",mock_input)
1028 1027
     stub(preprocessData,"qc_input",mock_input,2)
1029 1028
 
... ...
@@ -1037,7 +1036,7 @@ test_that("preprocessData QC, PTM and PTMTMT: Yes", {
1037 1036
         mockery::stub(MSstatsShiny::tmt_summarization_loop, "remove_modal_progress", NULL);
1038 1037
         MSstatsShiny::tmt_summarization_loop(...)})
1039 1038
 
1040
-    output <- preprocessData(mock_input,mock_input,MSstatsPTM::raw.input.tmt)
1039
+    output <- preprocessData(mock_input,mock_input,getData(mock_input))
1041 1040
     expected_names <- c("PTM","PROTEIN")
1042 1041
     expect_type(output,"list")
1043 1042
     expect_identical(names(output), expected_names)
... ...
@@ -1083,7 +1082,7 @@ test_that("preprocessData QC Other", {
1083 1082
     mock_input$norm = "equalizeMedians"
1084 1083
     mock_input$log = "2"
1085 1084
 
1086
-    # stub(preprocessData,"getData",mockGetData(mock_input))
1085
+    stub(preprocessData,"getData",mockGetData(mock_input))
1087 1086
     stub(preprocessData,"loadpage_input",mock_input)
1088 1087
     stub(preprocessData,"qc_input",mock_input,2)
1089 1088
 
... ...
@@ -1097,7 +1096,7 @@ test_that("preprocessData QC Other", {
1097 1096
         mockery::stub(MSstatsShiny::lf_summarization_loop, "remove_modal_progress", NULL);
1098 1097
         MSstatsShiny::lf_summarization_loop(...)})
1099 1098
 
1100
-    output <- preprocessData(mock_input,mock_input,MSstats::DDARawData)
1099
+    output <- preprocessData(mock_input,mock_input,getData(mock_input))
1101 1100
     expected_names <- c("FeatureLevelData","ProteinLevelData","SummaryMethod")
1102 1101
     expect_type(output,"list")
1103 1102
     expect_identical(names(output), expected_names)
... ...
@@ -1190,7 +1189,7 @@ test_that("dataComparison statmodel PTM PTMTMT: Yes", {
1190 1189
         mockery::stub(MSstatsShiny::tmt_summarization_loop, "remove_modal_progress", NULL);
1191 1190
         MSstatsShiny::tmt_summarization_loop(...)})
1192 1191
 
1193
-    output <- dataComparison(mock_input,mock_input,mock_input,dummy_matrix,preprocessData(mock_input,mock_input,MSstatsPTM::raw.input.tmt))
1192
+    output <- dataComparison(mock_input,mock_input,mock_input,dummy_matrix,preprocessData(mock_input,mock_input,getData(mock_input)))
1194 1193
     expected_names <- c("PTM.Model","PROTEIN.Model","ADJUSTED.Model")
1195 1194
     expect_type(output,"list")
1196 1195
     expect_identical(names(output), expected_names)
... ...
@@ -1238,7 +1237,7 @@ test_that("dataComparison statmodel PTM PTMTMT: No", {
1238 1237
         mockery::stub(MSstatsShiny::lf_model, "remove_modal_progress", NULL);
1239 1238
         MSstatsShiny::lf_model(...)})
1240 1239
 
1241
-    output <- dataComparison(mock_input,mock_input,mock_input,dummy_matrix,preprocessData(mock_input,mock_input,MSstatsPTM::raw.input))
1240
+    output <- dataComparison(mock_input,mock_input,mock_input,dummy_matrix,preprocessData(mock_input,mock_input,getData(mock_input)))
1242 1241
     expected_names <- c("PTM.Model","PROTEIN.Model","ADJUSTED.Model")
1243 1242
     expect_type(output,"list")
1244 1243
     expect_identical(names(output), expected_names)
... ...
@@ -1341,7 +1340,7 @@ test_that("dataComparison statmodel Other", {
1341 1340
         mockery::stub(MSstatsShiny::lf_summarization_loop, "remove_modal_progress", NULL);
1342 1341
         MSstatsShiny::lf_summarization_loop(...)})
1343 1342
 
1344
-    output <- dataComparison(mock_input,mock_input,mock_input,dummy_matrix,preprocessData(mock_input,mock_input,MSstats::DDARawData))
1343
+    output <- dataComparison(mock_input,mock_input,mock_input,dummy_matrix,preprocessData(mock_input,mock_input,getData(mock_input)))
1345 1344
     expected_names <- c("ComparisonResult","ModelQC","FittedModel")
1346 1345
     expect_type(output,"list")
1347 1346
     expect_identical(names(output), expected_names)