Browse code

propogate changes from add_edge_data to KL_compare and KEGG_lincs

Shana White authored on 30/10/2017 18:13:30
Showing 3 changed files

... ...
@@ -117,7 +117,8 @@ function(pathwayid, cell_line = NA,
117 117
                                 data_type = data_type, pert_time = pert_time)
118 118
         if (!is.na(user_data)[1,1]) {
119 119
             edges_plus_data <- add_edge_data(expanded_edges, KEGG_mappings, 
120
-                                            user_data, c(11,13), 
120
+                                            user_data, c(11,13),
121
+                                            map_type = "SYMBOL", 
121 122
                                             only_mapped = only_mapped)
122 123
             edge_map <- edge_mapping_info(edges_plus_data, data_added = TRUE, 
123 124
                                     significance_markup = significance_markup,
... ...
@@ -7,14 +7,10 @@
7 7
 #' SHSY5Y,SKL,SW480,VCAP)
8 8
 #' @param cell_line2 A cell line such that cell_line1 != cell_line2
9 9
 #' @param refine_by_cell_line A logical indicator
10
-#' @param add_L1000_edge_data A logical indicator 
11 10
 #' @param data_type Choose from data types: (100_full, 100_bing, 50_lm)
12 11
 #' @param pert_time Choose from (6,24,48,96,120,144,168)
13 12
 #' @param only_mapped A logical indicator; if set to FALSE will return 'de-novo'
14 13
 #' edges that 'exist' in data but are not documented in KEGG
15
-#' @param significance_markup A logical indicator; if set to TRUE will color
16
-#'  edges based on direction and significance of correlation (as determined by 
17
-#'  user-data-analysis)
18 14
 #' @param layered_nodes A logical indicator; if set to TRUE will create a graph 
19 15
 #' with 'stacked' nodes that the user can manipulate when multiple nodes are 
20 16
 #' mapped to one location
... ...
@@ -36,8 +32,6 @@
36 32
 KL_compare <-
37 33
     function(pathwayid, cell_line1 = NA, cell_line2 = NA,
38 34
              refine_by_cell_line = TRUE,
39
-             add_L1000_edge_data = TRUE,  
40
-             significance_markup = TRUE,
41 35
              data_type = "100_full",
42 36
              pert_time = 96,
43 37
              only_mapped = TRUE,
... ...
@@ -144,18 +138,29 @@ KL_compare <-
144 138
         if (!is.na(cell_data1)[1,1] & !is.na(cell_data2)[1,1]) {
145 139
             edges_plus_data1 <- add_edge_data(expanded_edges, KEGG_mappings, 
146 140
                                               cell_data1, c(15,16), 
141
+                                              map_type = "SYMBOL",
147 142
                                               only_mapped = only_mapped)
148 143
             edges_plus_data2 <- add_edge_data(expanded_edges, KEGG_mappings, 
149 144
                                               cell_data2, c(15,16), 
145
+                                              map_type = "SYMBOL",
150 146
                                               only_mapped = only_mapped)
151
-            edges_plus_data1$unique_ID <- paste0(edges_plus_data1$entry1symbol, ":", edges_plus_data1$entry2symbol, ":", edges_plus_data1$edgeID)
152
-            edges_plus_data2$unique_ID <- paste0(edges_plus_data2$entry1symbol, ":", edges_plus_data2$entry2symbol, ":", edges_plus_data2$edgeID)
147
+            names(edges_plus_data1)[18] <- "SE_1"
148
+            names(edges_plus_data2)[18] <- "SE_2"
149
+            edges_plus_data1$unique_ID <- paste0(edges_plus_data1$entry1symbol, 
150
+                                                 ":", edges_plus_data1$entry2symbol, 
151
+                                                 ":", edges_plus_data1$edgeID)
152
+            edges_plus_data2$unique_ID <- paste0(edges_plus_data2$entry1symbol, 
153
+                                                 ":", edges_plus_data2$entry2symbol, 
154
+                                                 ":", edges_plus_data2$edgeID)
153 155
             
154
-            edges_compare <- merge(edges_plus_data1, edges_plus_data2[,c(18,19,21)], by = "unique_ID")
156
+            edges_compare <- merge(edges_plus_data1, edges_plus_data2[,c(18,19,21)], 
157
+                                    by = "unique_ID")
155 158
             
156 159
             for(i in 1:nrow(edges_compare)){
157
-                if (!is.na(edges_compare$log_OR_1[i]) & !is.na(edges_compare$log_OR_2[i])){
158
-                    edges_compare$test[i] <- (edges_compare$log_OR_1[i] - edges_compare$log_OR_2[i])/
160
+                if (!is.na(edges_compare$log_OR_1[i]) & 
161
+                    !is.na(edges_compare$log_OR_2[i])){
162
+                    edges_compare$test[i] <- (edges_compare$log_OR_1[i] - 
163
+                                                  edges_compare$log_OR_2[i])/
159 164
                         sqrt(edges_compare$SE_1[i]^2 + edges_compare$SE_2[i]^2)
160 165
                     edges_compare$summary_score[i] <- exp(abs(edges_compare$test[i]))
161 166
                     
... ...
@@ -5,8 +5,7 @@
5 5
 \title{Combines all other package functions for one-step cell line comparison}
6 6
 \usage{
7 7
 KL_compare(pathwayid, cell_line1 = NA, cell_line2 = NA,
8
-  refine_by_cell_line = TRUE, add_L1000_edge_data = TRUE,
9
-  significance_markup = TRUE, data_type = "100_full", pert_time = 96,
8
+  refine_by_cell_line = TRUE, data_type = "100_full", pert_time = 96,
10 9
   only_mapped = TRUE, get_data = FALSE, convert_KEGG_IDs = TRUE,
11 10
   graph_title = "default", tidy_edge = TRUE, layered_nodes = FALSE)
12 11
 }
... ...
@@ -22,12 +21,6 @@ SHSY5Y,SKL,SW480,VCAP)}
22 21
 
23 22
 \item{refine_by_cell_line}{A logical indicator}
24 23
 
25
-\item{add_L1000_edge_data}{A logical indicator}
26
-
27
-\item{significance_markup}{A logical indicator; if set to TRUE will color
28
-edges based on direction and significance of correlation (as determined by 
29
-user-data-analysis)}
30
-
31 24
 \item{data_type}{Choose from data types: (100_full, 100_bing, 50_lm)}
32 25
 
33 26
 \item{pert_time}{Choose from (6,24,48,96,120,144,168)}