commit/galaxy-central: kanwei: Update test-data files for meme 4.6.0
1 new changeset in galaxy-central: http://bitbucket.org/galaxy/galaxy-central/changeset/e448ba0e8843/ changeset: r5129:e448ba0e8843 user: kanwei date: 2011-02-25 20:58:35 summary: Update test-data files for meme 4.6.0 affected #: 4 files (2.8 KB) --- a/test-data/meme/meme/meme_output_html_1.html Fri Feb 25 13:35:27 2011 -0500 +++ b/test-data/meme/meme/meme_output_html_1.html Fri Feb 25 14:58:35 2011 -0500 @@ -4,6 +4,8 @@ <meta http-equiv="Content-Type" content="text/html; charset=UTF-8"><title>MEME</title><style type="text/css"> + + /* START INCLUDED FILE "meme.css" */ /* The following is the content of meme.css */ body { background-color:white; font-size: 12px; font-family: Verdana, Arial, Helvetica, sans-serif;} @@ -13,6 +15,7 @@ padding: 0px; width: 12px; height: 13px; + cursor: pointer; background-image: url("help.gif"); background-image: url("data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAAwAAAANAQMAAACn5x0BAAAAAXNSR0IArs4c6QAAAAZQTFRFAAAAnp6eqp814gAAAAF0Uk5TAEDm2GYAAAABYktHRACIBR1IAAAACXBIWXMAAAsTAAALEwEAmpwYAAAAB3RJTUUH2gMJBQgGYqhNZQAAACZJREFUCNdj+P+BoUGAoV+AYeYEEGoWYGgTYGgRAAm2gRGQ8f8DAOnhC2lYnqs6AAAAAElFTkSuQmCC"); } @@ -123,8 +126,12 @@ padding:0px 10px; } + /* END INCLUDED FILE "meme.css" */ + </style><script type="text/javascript"> + + /* START INCLUDED FILE "motif_logo.js" */ //====================================================================== // start Alphabet object //====================================================================== @@ -337,7 +344,9 @@ //====================================================================== // start Pspm object //====================================================================== - function Pspm(pssm, name) { + function Pspm(pssm, name, ltrim, rtrim) { + if (ltrim === undefined) ltrim = 0; + if (rtrim === undefined) rtrim = 0; //variable prototype this.alph_length = 0; this.motif_length = 0; @@ -345,6 +354,8 @@ this.name = (typeof name == "string" ? name : ""); this.nsites = 0; this.evalue = 0; + this.ltrim = ltrim; + this.rtrim = rtrim; //function prototype this.copy = Pspm_copy; this.reverse_complement = Pspm_reverse_complement; @@ -352,10 +363,11 @@ this.get_stack_ic = Pspm_get_stack_ic; this.get_motif_length = Pspm_get_motif_length; this.get_alph_length = Pspm_get_alph_length; + this.get_left_trim = Pspm_get_left_trim; + this.get_right_trim = Pspm_get_right_trim; this.toString = Pspm_to_string; //construct var pspm_header = /letter-probability matrix:\s+alength=\s+(\d+)\s+w=\s+(\d+)(\s+nsites=\s+(\S+))?(\s+E=\s+(\S+))?\s*/; - var is_prob = /^((1(\.0+)?)|(0(\.\d+)?))$/; var is_empty = /^\s*$/; var lines = pssm.split(/\s*\n\s*/); var read_pssm = false; @@ -392,11 +404,14 @@ if (!parts.hasOwnProperty(part_index) || parts[part_index] === undefined) continue; var prob = parts[part_index]; - if (!is_prob.test(prob)) continue; if (col_num >= this.alph_length) { throw "TOO_MANY_COLS"; } this.pspm[line_num][col_num] = (+prob); + //check the probability is within bounds + if (this.pspm[line_num][col_num] > 1 || this.pspm[line_num][col_num] < 0) { + throw "NUM_NOT_PROB"; + } col_num++; } if (col_num != this.alph_length) { @@ -417,6 +432,11 @@ //so far only a shallow copy, need to copy everything clone.alph_length = (0+this.alph_length); clone.motif_length = (0+this.motif_length); + clone.name = (""+this.name); + clone.nsites = (0+this.nsites); + clone.evalue = (0+this.evalue); + clone.ltrim = (0+this.ltrim); + clone.rtrim = (0+this.rtrim); clone.pspm = new Array(this.motif_length); for (row = 0; row < this.motif_length; row++) { clone.pspm[row] = new Array(this.alph_length); @@ -460,6 +480,10 @@ row[c_index] = row[g_index]; row[g_index] = temp; } + //swap triming + var temp_trim = this.ltrim; + this.ltrim = this.rtrim; + this.rtrim = temp_trim; //note that ambigs are ignored because they don't effect motifs return this; //allow function chaining... } @@ -515,6 +539,14 @@ return this.alph_length; } + function Pspm_get_left_trim() { + return this.ltrim; + } + + function Pspm_get_right_trim() { + return this.rtrim; + } + function Pspm_to_string() { var str = ""; for (row_index in this.pspm) { @@ -1036,6 +1068,62 @@ pos = pos - sym_height; } } + + function draw_dashed_line(ctx, pattern, start, x1, y1, x2, y2) { + var x, y, len, i; + var dx = x2 - x1; + var dy = y2 - y1; + var tlen = Math.pow(dx*dx + dy*dy, 0.5); + var theta = Math.atan2(dy,dx); + var mulx = Math.cos(theta); + var muly = Math.sin(theta); + var lx = []; + var ly = []; + for (i = 0; i < pattern; ++i) { + lx.push(pattern[i] * mulx); + ly.push(pattern[i] * muly); + } + i = start; + x = x1; + y = y1; + len = 0; + ctx.beginPath(); + while (len + pattern[i] < tlen) { + ctx.moveTo(x, y); + x += lx[i]; + y += ly[i]; + ctx.lineTo(x, y); + len += pattern[i]; + i = (i + 1) % pattern.length; + x += lx[i]; + y += ly[i]; + len += pattern[i]; + i = (i + 1) % pattern.length; + } + if (len < tlen) { + ctx.moveTo(x, y); + x += mulx * (tlen - len); + y += muly * (tlen - len); + ctx.lineTo(x, y); + } + ctx.stroke(); + } + + function draw_trim_background(ctx, metrics, pspm, offset) { + var lwidth = metrics.stack_width * pspm.get_left_trim(); + var rwidth = metrics.stack_width * pspm.get_right_trim(); + var mwidth = metrics.stack_width * pspm.get_motif_length(); + var rstart = mwidth - rwidth; + ctx.save();//s8 + ctx.translate(offset * metrics.stack_width, 0); + ctx.fillStyle = "rgb(240, 240, 240)"; + if (pspm.get_left_trim() > 0) ctx.fillRect(0, 0, lwidth, metrics.stack_height); + if (pspm.get_right_trim() > 0) ctx.fillRect(rstart, 0, rwidth, metrics.stack_height); + ctx.fillStyle = "rgb(51, 51, 51)"; + if (pspm.get_left_trim() > 0) draw_dashed_line(ctx, [3], 0, lwidth-0.5, 0, lwidth-0.5, metrics.stack_height); + if (pspm.get_right_trim() > 0) draw_dashed_line(ctx, [3], 0, rstart+0.5, 0, rstart+0.5, metrics.stack_height); + ctx.restore();//s8 + } function draw_logo_on_canvas(logo, canvas, show_names, scale) { var draw_name = (typeof show_names == "boolean" ? show_names : (logo.get_rows() > 1)); @@ -1089,6 +1177,10 @@ //translate across past the scale ctx.translate(metrics.y_label_height + metrics.y_label_spacer + metrics.y_num_width + metrics.y_tic_width, 0); + //draw the trimming background + if (pspm.get_left_trim() > 0 || pspm.get_right_trim() > 0) { + draw_trim_background(ctx, metrics, pspm, offset); + } //draw letters ctx.translate(0, metrics.y_num_height / 2); for (var col_index = 0; col_index < logo.get_columns(); col_index++) { @@ -1211,15 +1303,8 @@ element.parentNode.replaceChild(canvas, element); } - //example code for creating a logo - //function do_load() { - // var alphabet = new Alphabet(document.getElementById("alphabet").value, document.getElementById("bgfreq").value); - // var pspm = new Pspm(document.getElementById("pspm1").value); - // var pspm_rc = pspm.copy().reverse_complement(alphabet); - // draw_logo(pspm, alphabet, "Motif 1", "logo1"); - // draw_logo(pspm_rc, alphabet, "Motif 1 Reverse Complement", "logo_rc1"); - //} - + /* END INCLUDED FILE "motif_logo.js" */ + function setup() { var motif_count = 1; @@ -1377,15 +1462,15 @@ } </script></head> -<body onload="javascript:setup()"><form enctype="application/x-www-form-urlencoded" method="post" target="_new" action="http://dan/meme_request.cgi"> +<body onload="javascript:setup()"><form enctype="application/x-www-form-urlencoded" method="post" target="_new" action="http://straub/meme/cgi-bin/meme_request.cgi"><!--+++++++++++++++++++++++START DATA+++++++++++++++++++++++++++++++++++++--> -<input type="hidden" name="version" value="MEME version 4.5.0"> +<input type="hidden" name="version" value="MEME version 4.6.0"><input type="hidden" name="alphabet" id="alphabet" value="ACDEFGHIKLMNPQRSTVWY"><input type="hidden" name="strands" value="none"><input type="hidden" name="bgfreq" id="bgfreq" value="A 0.291 C 0.229 D 0.001 E 0.001 F 0.001 G 0.255 H 0.001 I 0.001 K 0.001 L 0.001 M 0.001 N 0.001 P 0.001 Q 0.001 R 0.001 S 0.001 T 0.215 V 0.001 W 0.001 Y 0.001"> -<input type="hidden" name="name" value="dan.dat"> +<input type="hidden" name="name" value="dataset_3.dat"><input type="hidden" name="combinedblock" value=" chr21_19617074_19617124_+ 1.22e-03 1 50 1 39 3.06e-05 chr21_26934381_26934431_+ 2.21e-03 1 50 1 27 5.52e-05 @@ -2179,12 +2264,12 @@ </tbody></table></div> -<p>Time 1.6 secs.</p> +<p>Time 2.1 secs.</p></div><div class="bar"><div style="text-align:right;"><a href="#top_buttons">Top</a></div><div class="subsection"> -<a name="version"></a><h5>MEME version</h5>4.5.0 (Release date: Thu Oct 14 08:44:28 EST 2010) +<a name="version"></a><h5>MEME version</h5>4.6.0 (Release date: Thu Jan 20 14:06:48 PST 2011) </div><div class="subsection"><a name="reference"></a><h5>Reference</h5> @@ -2196,7 +2281,7 @@ </div><a name="sequences"></a><a href="javascript:showHidden('trainingseq')" id="trainingseq_activator">show training set...</a><div class="subsection" id="trainingseq_data" style="display:none;"><h5>Training set:</h5> - Datafile = dan.dat<br> + Datafile = dataset_3.dat<br> Alphabet = ACDEFGHIKLMNPQRSTVWY<br><table><col><col> @@ -2387,19 +2472,19 @@ </div><a href="javascript:hideShown('trainingseq')" style="display:none;" id="trainingseq_deactivator">hide training set...</a><div class="subsection"><a name="command"></a><h5>Command line summary</h5> -<textarea rows="1" style="width:100%;" readonly>meme dan.dat -o /dataset_dan_files -nostatus </textarea><br> +<textarea rows="1" style="width:100%;" readonly>meme dataset_3.dat -o dataset_4_files -nostatus </textarea><br> Letter frequencies in dataset<br><div style="margin-left:25px;">A: 0.294 C: 0.231 D: 0.000 E: 0.000 F: 0.000 G: 0.257 H: 0.000 I: 0.000 K: 0.000 L: 0.000 M: 0.000 N: 0.000 P: 0.000 Q: 0.000 R: 0.000 S: 0.000 T: 0.217 V: 0.000 W: 0.000 Y: 0.000</div><br> Background letter frequencies (from dataset with add-one prior applied):<br><div style="margin-left:25px;">A: 0.291 C: 0.229 D: 0.001 E: 0.001 F: 0.001 G: 0.255 H: 0.001 I: 0.001 K: 0.001 L: 0.001 M: 0.001 N: 0.001 P: 0.001 Q: 0.001 R: 0.001 S: 0.001 T: 0.215 V: 0.001 W: 0.001 Y: 0.001</div><br></div><div class="subsection"> -<a name="stopped"></a><h5>Stopping Reason</h5>Stopped because nmotifs = 1 reached. Program ran on <i>dan</i>. +<a name="stopped"></a><h5>Stopping Reason</h5>Stopped because nmotifs = 1 reached. Program ran on <i>scofield</i>. </div><a href="javascript:showHidden('model')" id="model_activator">show model parameters...</a><div class="subsection" id="model_data" style="display:none;"><h5>Model parameters</h5><textarea style="width:100%;" rows="28" readonly> -host = dan +host = scofield type = zoops nmotifs = 1 evalue_threshold = inf --- a/test-data/meme/meme/meme_output_txt_1.txt Fri Feb 25 13:35:27 2011 -0500 +++ b/test-data/meme/meme/meme_output_txt_1.txt Fri Feb 25 14:58:35 2011 -0500 @@ -1,7 +1,7 @@ ******************************************************************************** MEME - Motif discovery tool ******************************************************************************** -MEME version 4.5.0 (Release date: Thu Oct 14 08:44:28 EST 2010) +MEME version 4.6.0 (Release date: Thu Jan 20 14:06:48 PST 2011) For further information on how to interpret these results or to get a copy of the MEME software please access http://meme.nbcr.net. @@ -28,7 +28,7 @@ ******************************************************************************** TRAINING SET ******************************************************************************** -DATAFILE= dan.dat +DATAFILE= dataset_3.dat ALPHABET= ACDEFGHIKLMNPQRSTVWY Sequence name Weight Length Sequence name Weight Length ------------- ------ ------ ------------- ------ ------ @@ -55,7 +55,7 @@ This information can also be useful in the event you wish to report a problem with the MEME software. -command: meme dan.dat -o dan/dataset_dan_files -nostatus +command: meme dataset_3.dat -o dataset_4_files -nostatus model: mod= zoops nmotifs= 1 evt= inf object function= E-value of product of p-values @@ -267,7 +267,7 @@ -Time 1.56 secs. +Time 2.13 secs. ******************************************************************************** @@ -320,6 +320,6 @@ Stopped because nmotifs = 1 reached. ******************************************************************************** -CPU: dan +CPU: scofield ******************************************************************************** --- a/test-data/meme/meme/meme_output_xml_1.xml Fri Feb 25 13:35:27 2011 -0500 +++ b/test-data/meme/meme/meme_output_xml_1.xml Fri Feb 25 14:58:35 2011 -0500 @@ -157,8 +157,8 @@ ]><!-- Begin document body --> -<MEME version="4.5.0" release="Thu Oct 14 08:44:28 EST 2010"> -<training_set datafile="/dan_110.dat" length="30"> +<MEME version="4.6.0" release="Thu Jan 20 14:06:48 PST 2011"> +<training_set datafile="dataset_3.dat" length="30"><alphabet id="amino-acid" length="20"><letter id="letter_A" symbol="A"/><letter id="letter_C" symbol="C"/> @@ -183,6 +183,8 @@ </alphabet><ambigs><letter id="letter_B" symbol="B"/> +<letter id="letter_J" symbol="J"/> +<letter id="letter_O" symbol="O"/><letter id="letter_U" symbol="U"/><letter id="letter_X" symbol="X"/><letter id="letter_Z" symbol="Z"/> @@ -245,8 +247,8 @@ </letter_frequencies></training_set><model> -<command_line>meme /dan.dat -o /dataset_dan_files -nostatus </command_line> -<host>dan</host> +<command_line>meme dataset_3.dat -o dataset_4_files -nostatus </command_line> +<host>scofield</host><type>zoops</type><nmotifs>1</nmotifs><evalue_threshold>inf</evalue_threshold> @@ -300,7 +302,7 @@ </background_frequencies></model><motifs> -<motif id="motif_1" name="1" width="11" sites="25" ic="40.0" re="13.8" llr="239" e_value="2.4e-011" bayes_threshold="5.33554" elapsed_time="1.563760"> +<motif id="motif_1" name="1" width="11" sites="25" ic="40.0" re="13.8" llr="239" e_value="2.4e-011" bayes_threshold="5.33554" elapsed_time="2.128133"><scores><alphabet_matrix><alphabet_array> --- a/tools/meme/meme.xml Fri Feb 25 13:35:27 2011 -0500 +++ b/tools/meme/meme.xml Fri Feb 25 14:58:35 2011 -0500 @@ -319,7 +319,7 @@ <param name="input1" value="meme_input_1.fasta" ftype="fasta" dbkey="hg19"/><param name="options_type_selector" value="basic"/><param name="non_commercial_use" value="True"/> - <output name="html_outfile" file="meme/meme/meme_output_html_1.html" lines_diff="14"/> + <output name="html_outfile" file="meme/meme/meme_output_html_1.html" lines_diff="12"/><output name="txt_outfile" file="meme/meme/meme_output_txt_1.txt" lines_diff="12"/><output name="xml_outfile" file="meme/meme/meme_output_xml_1.xml" lines_diff="8"/></test> Repository URL: https://bitbucket.org/galaxy/galaxy-central/ -- This is a commit notification from bitbucket.org. You are receiving this because you have the service enabled, addressing the recipient of this email.
participants (1)
-
Bitbucket