ROC - Receiver Operator Characteristic

ROC is a program for computing the Receiver Operator Characteristic for a sequence analsysis program. The program can also produce plots of the percentages of the negative and positive results as the score threshhold is varied.

Source files you may copy include:

• roc.f90, the source code;

ROC needs a search result file to analyze. Four "generic" files are included, which are tiny examples of "fake data" that make it easier to understand. Four other files of real data will show you how complicated it can be to extract the information we need. Search result files you may copy include:

• GENERIC_01.TXT, generic file #1;
• GENERIC_02.TXT, generic file #2;
• GENERIC_03.TXT, generic file #3;
• GENERIC_04.TXT, generic file #4;
• BLAST_PA2.TXT, results of a BLAST search on PA2;
• FASTA_PA2.TXT, results of a FASTA search on PA2;
• MAXSEGS_PA2.TXT, results of a MAXSEGS search on PA2;
• PEARSON_PA2.TXT, results of a PEARSON search on PA2;

Along with a search result file, ROC needs a file that lists the "good" items, that is, the items that do belong to the class of interest. (This file may include many extra names that don't show up in the search list.) "Good" files you may copy include:

• GOOD_01.TXT, for generic search #1;
• GOOD_02.TXT, for generic search #2;
• GOOD_03.TXT, for generic search #3;
• GOOD_04.TXT, for generic search #4;
• GOOD_PA2.TXT, for searches on PA2;

The list of routines includes:

• ROC controls the computation of the Receiver Operator Characteristic.
• ANALYSIS carries out the ROC analysis.
• ANALYSIS_PARAMS_GET gets the ROC analysis parameters.
• ANALYSIS_PARAMS_PRINT prints the analysis parameters.
• ANALYSIS_POST does some post-analysis tasks.
• ANALYSIS_PRE carries out some pre-analysis tasks.
• ANALYSIS_PRINT prints the analysis.
• CH_CAP capitalizes a single character.
• CH_EQI is a case insensitive comparison of two characters for equality.
• CH_TO_DIGIT returns the integer value of a base 10 digit.
• CHVEC_PERMUTE permutes a character vector in place.
• CHVEC_REVERSE reverses the elements of a character vector.
• DATA_COUNT counts the positive and negative data.
• DATA_COUNT_PRINT prints the data counts.
• DATA_GET retrieves the identifiers, scores, and good identifiers.
• DATA_PARAMS_GET gets the data parameters.
• DATA_PARAMS_PRINT prints the data parameters.
• DATA_PRINT prints the data.
• DATA_SELECT selects the data to be analyzed.
• DIGIT_TO_CH returns the character representation of a decimal digit.
• FILE_ADVANCE_TO_STRING searches ahead in a text file for a string.
• GET_UNIT returns a free FORTRAN unit number.
• GOOD_FILE_READ reads a "good" ID file.
• GOOD_GET retrieves the good identifiers.
• GOOD_PRINT prints the good identifiers.
• I_TO_S_LEFT converts an integer to a left-justified string.
• IVEC_CUM computes the cumulutive sum of the entries of a vector.
• IVEC_IDENTITY sets an integer vector to the identity vector A(I)=I.
• IVEC_PERMUTE permutes an integer vector in place.
• IVEC_REVERSE reverses the elements of an integer vector.
• IVEC_NEGATIVE_INDEX indexes negative integer vector entries.
• IVEC_POSITIVE_INDEX indexes positive integer vector entries.
• LGTI = STRNG1 is lexically greater than STRNG2.
• LLEI = STRNG1 is lexically less than or equal to STRNG2.
• LLTI = STRNG1 is lexically less than STRNG2.
• PARAMS_DEFAULT sets the parameters to default values.
• PERM_CHECK checks that a vector represents a permutation.
• PN_GRAPH_FILE_WRITE creates a positive/negative graphics file of the results.
• PN_ROC_INT computes the ROC integral.
• PNS_DATA_SET sets the data for a positive/negative/score graph.
• PNS_GRAPH_FILE_WRITE creates a positive/negative/score graphics file of the r
• R_SWAP switches two real values.
• R_TO_S_LEFT writes a real into a left justified character string.
• RVEC_BIN computes bins based on a given real vector.
• RVEC_MINMAX returns the minimum and maximum values in a real vector.
• RVEC_ORDER_TYPE determines if a real array is (non)strictly ascending/descending.
• RVEC_PERMUTE permutes a real vector in place.
• RVEC_REVERSE reverses the elements of a real vector.
• RVEC_SORT_HEAP_INDEX_A does an indexed heap ascending sort of a real vector.
• S_BLANKS_DELETE replaces consecutive blanks by one blank.
• S_CAP replaces any lowercase letters by uppercase ones in a string.
• S_EQI is a case insensitive comparison of two strings for equality.
• S_GTI = S1 is lexically greater than S2.
• S_LEFT flushes a string left.
• S_LTI = S1 is lexically less than S2.
• S_SWAP switches two strings.
• S_TO_I reads an integer value from a string.
• S_TO_R reads a real number from a string.
• SCORE_BLAST_FILE_READ extracts information from a BLAST search file.
• SCORE_FASTA_FILE_READ extracts information from a FASTA search file.
• SCORE_GENERIC_FILE_READ reads a "name, score" file.
• SCORE_MAXSEGS_FILE_READ extracts information from a MAXSEGS search file.
• SCORE_PEARSON_FILE_READ extracts information from a PEARSON search file.
• SORT_HEAP_EXTERNAL externally sorts a list of items into linear order.
• SVEC_PERMUTE permutes a string vector in place.
• SVEC_REVERSE reverses the elements of a string vector.
• SVECI_SEARCH_BINARY_A searches an ascending sorted vector of implicitly capitalized strings.
• SVECI_SORT_HEAP sorts a vector of implicitly capitalized strings using heap s
• TIMESTAMP prints the current YMDHMS date as a time stamp.
• WORD_FIND finds the word of a given index in a string.

Return to the biomedical software page.