Marinka Zitnik avatar Marinka Zitnik committed d4877c3

include standardization in testing significance of intensity profiles

Comments (0)

Files changed (2)

ypta/analysis/patches.py

     strain2inten_profiles = {}
     for strn, ip in inten_profiles:
         ip = np.array(ip)
+        # standardize intensity profiles of all tracks in strain
+        ip = utils.standardize(ip)
         ip[np.isnan(ip)] = 0
         strain2inten_profiles[strn] = ip
+    # standardize intensity profiles at the level of subplate (batch effects)    
     # divide strains into two groups (wt, non-wt)
     wt_strains, nwt_strains = utils.partition_subplate(strain2orf, wt_name)
     # perform stats test

ypta/run_analysis.py

 
 if 0:
     data = utils.load(data_dir)
-    utils.save(data, res_dir)
-    filtered_data = utils.filtering(data)
+    filtered_data = utils.filtering(data)
+    utils.save(filtered_data, res_dir)
Tip: Filter by directory path e.g. /media app.js to search for public/media/app.js.
Tip: Use camelCasing e.g. ProjME to search for ProjectModifiedEvent.java.
Tip: Filter by extension type e.g. /repo .js to search for all .js files in the /repo directory.
Tip: Separate your search with spaces e.g. /ssh pom.xml to search for src/ssh/pom.xml.
Tip: Use ↑ and ↓ arrow keys to navigate and return to view the file.
Tip: You can also navigate files with Ctrl+j (next) and Ctrl+k (previous) and view the file with Ctrl+o.
Tip: You can also navigate files with Alt+j (next) and Alt+k (previous) and view the file with Alt+o.