# HG changeset patch # User bornea # Date 1461701816 14400 # Node ID e6e456d3ac1486b2750cb493715312aee8d23401 # Parent bd71998aec8d1b7a783f854e39359dcba74df980 Uploaded diff -r bd71998aec8d -r e6e456d3ac14 pre_process_protein_name_set.R --- a/pre_process_protein_name_set.R Tue Apr 26 15:40:35 2016 -0400 +++ b/pre_process_protein_name_set.R Tue Apr 26 16:16:56 2016 -0400 @@ -88,7 +88,9 @@ } protein_intensities_tukeys = get_protein_values(peptides_txt_mapped_log2,intensity_columns) protein_intensities_tukeys[protein_intensities_tukeys == 1] <- 0 - write.table(protein_intensities_tukeys, "./tukeys_output.txt", row.names = FALSE, col.names = TRUE, quote = FALSE, sep = "\t") + # Eliminates all pure NA rows. + protein_intensities_tukeys = protein_intensities_tukeys[rowSums(is.na(protein_intensities_tukeys)) != (ncol(protein_intensities_tukeys)-1),] + write.table(protein_intensities_tukeys, "./tukeys_output.txt", row.names = FALSE, col.names = TRUE, quote = FALSE, sep = "\t") }