Skip to content

Commit c6eb7ff

Browse files
beckdacbeckdac
committed
new analyses
1 parent 6ff5e38 commit c6eb7ff

12 files changed

Lines changed: 158 additions & 0 deletions

‎16.top_pop_figure‎

Lines changed: 21 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,21 @@
1+
#!/bin/bash
2+
3+
source globals
4+
5+
DAT=$0.tab
6+
IN=13.figures.otus_percent.tab
7+
8+
OTUS="OTU_115 OTU_4 OTU_3 OTU_1 OTU_239 OTU_5 OTU_2"
9+
10+
head -1 $IN | awk '{ printf("OTUid\t%s\n", $0); }' > $DAT
11+
12+
for otu in $OTUS
13+
do
14+
echo $otu
15+
awk -F'\t' -v otu="$otu" '{ if ($1 == otu) { print $0; exit } }' $IN >> $DAT
16+
done
17+
18+
\rm -rf $0.15uM.tab $0.45uM.tab $0.75uM.tab $0.150uM.tab $0.225uM.tab
19+
awk -F'\t' -f $0.awk $DAT
20+
21+
/gscratch/esci/dacb/gnuplot/bin/gnuplot $0.gplt

‎16.top_pop_figure.150uM.tab‎

Lines changed: 7 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,7 @@
1+
OTUid OTU_115 OTU_4 OTU_3 OTU_1 OTU_239 OTU_5 OTU_2
2+
O2_50A_10 0.000246404041026273 0.00221763636923646 0.0522992577078264 0.190439523208181 0.0874426340591986 0.578649089845073 0.0210675455077463
3+
O2_50B_10 0.0198138357959423 0.00449595589896525 0.00718449234106005 0.501830012200081 0.0672585965387917 0.244453481541729 0.0155438073290859
4+
O2_50C_10 0.00175216866777733 0.00206813350950767 0.00301602803469868 0.161659102659849 0.00318837249382432 0.165307060378009 0.0178663755960246
5+
O2_50A_15 0.00159716060337178 0.00301685891748004 0.143655723158829 0.37435669920142 0.0375332741792369 0.225022182786158 0.033717834960071
6+
O2_50B_15 0.020653557787523 0.00841703211203848 0.0280096194652709 0.517682840571509 0.013085302022917 0.135238364690904 0.0326778893761494
7+
O2_50C_15 0.000475737392959087 0.00445463195225327 0.0235706253784275 0.710319176541822 0.00086497707810743 0.00320041518899749 0.0524608597872156

‎16.top_pop_figure.15uM.tab‎

Lines changed: 7 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,7 @@
1+
OTUid OTU_115 OTU_4 OTU_3 OTU_1 OTU_239 OTU_5 OTU_2
2+
O2_5A_10 0.0131554015645725 0.469507078472241 0.0212489510868094 0.00646942587228974 0.00108274909996481 0.00476409603984517 0.349809165471131
3+
O2_5B_10 0.000460809268849293 0.525278679891161 0.0595979987711753 0.0073071184060388 0.00111910822434828 0.00441060300184324 0.303256385499868
4+
O2_5C_10 0.000962592582689377 0.0130217385491591 0.0201877055536244 0.153132436696168 0.00264712960239579 0.00981309660686115 0.735420733174684
5+
O2_5A_15 0.000780284803953443 0.108654658950517 0.381884387801548 0.0121594381949412 0.000520189869302295 0.00435659015540672 0.353208921256259
6+
O2_5B_15 0.00128609844499006 0.241786507658132 0.233309949725243 0.0115748860049106 0.00239681982929966 0.00485209867882614 0.36630422074126
7+
O2_5C_15 0.0250296208530806 0.199200236966825 0.0313240521327014 0.0301392180094787 0.00192535545023697 0.00651658767772512 0.663803317535545

‎16.top_pop_figure.225uM.tab‎

Lines changed: 7 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,7 @@
1+
OTUid OTU_115 OTU_4 OTU_3 OTU_1 OTU_239 OTU_5 OTU_2
2+
O2_75A_10 0.000157996603073034 0.00231070031994312 0.00315993206146068 0.656140142986926 0.00152071730457795 0.247758423193901 0.0142591934273413
3+
O2_75B_10 0.000194225042082092 0.00183434761966421 0.0637273943631577 0.419806638180327 0.306422374724848 0.149941732487375 0.0183003150761794
4+
O2_75C_10 0.00236502193204099 0.00238523579470801 0.00430555274807463 0.714741970043056 0.00119261789735401 0.196903236239413 0.0459056821167957
5+
O2_75A_15 0.000433254705627497 0.00226255235161026 0.0331680546863717 0.490588745005536 0.208491792230299 0.125355028161556 0.0613296105521591
6+
O2_75B_15 0 0 0.0155936094900985 0.510044963168468 0.169185879651775 0.220989189706304 0.0315698842437578
7+
O2_75C_15 0.00014808233377758 0.00133274100399822 0.0103657633644306 0.450614541685177 0.00251739967421887 0.424255886272768 0.031393454760847

‎16.top_pop_figure.45uM.tab‎

Lines changed: 7 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,7 @@
1+
OTUid OTU_115 OTU_4 OTU_3 OTU_1 OTU_239 OTU_5 OTU_2
2+
O2_15A_10 0.000203533338760889 0.00199462671985671 0.0491940079785069 0.488764959700399 0.00122120003256533 0.00588211349018969 0.410893104290483
3+
O2_15B_10 8.77943855490441e-05 0.00215096244595158 0.115427668400606 0.272821053093655 0.00173393911459362 0.00796734048857576 0.543908167072716
4+
O2_15C_10 0.0017664136556018 0.102032250166148 0.517769072020707 0.00834236944279268 0.000787015985169121 0.00363776277589283 0.234863059218581
5+
O2_15A_15 0.00193254497768371 0.00257672663691161 0.0346017576956702 0.391846500713201 0.00156444117241062 0.00630377766530162 0.46643353425666
6+
O2_15B_15 0.000241266164833044 0.00241266164833044 0.585118702953098 0.0123045744064852 0.00101331789229878 0.00414977803512835 0.342549700829956
7+
O2_15C_15 0.000368091707419677 0.00347057895567124 0.328916232844297 0.0226113477414945 0.00152495135931009 0.00625755902613451 0.592049219119735

‎16.top_pop_figure.75uM.tab‎

Lines changed: 7 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,7 @@
1+
OTUid OTU_115 OTU_4 OTU_3 OTU_1 OTU_239 OTU_5 OTU_2
2+
O2_25A_10 0.00441885715728606 0.0216145241522107 0.0199479837386057 0.00583289144761761 0.0012877812286948 0.00653990859278337 0.566371234502437
3+
O2_25B_10 0.00575535043165128 0.00567672542575441 0.00929347569701068 0.564684792351359 0.00110075008255626 0.00575535043165128 0.349928451244634
4+
O2_25C_10 0.000253292806484296 0.00405268490374873 0.183320668693009 0.206876899696049 0.00151975683890578 0.00698665991219183 0.515873015873016
5+
O2_25A_15 0.0525107917600131 0.0274301950713076 0.617452598218677 0.102180208731763 0.000764985519916944 0.00289601661111415 0.159881973662641
6+
O2_25B_15 0.00600555886440341 0.0130037720865595 0.0250645225332539 0.604973198332341 0.00198530871550526 0.00565812983918999 0.318145721659718
7+
O2_25C_15 0.013517359134889 0.0120470498956555 0.4314171883893 0.255169797002466 0.0011857332574464 0.00426863972680706 0.248624549421362

‎16.top_pop_figure.awk‎

Lines changed: 41 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,41 @@
1+
// {
2+
++line;
3+
for (i = 1; i <= NF; ++i) {
4+
data[line, i] = $i;
5+
}
6+
nf = NF;
7+
}
8+
END {
9+
for (i = 1; i <= nf; ++i) {
10+
split(data[1, i], a, "_");
11+
if (a[2] == "5A" || a[2] == "5B" || a[2] == "5C")
12+
of = "16.top_pop_figure.15uM.tab";
13+
else if (a[2] == "15A" || a[2] == "15B" || a[2] == "15C")
14+
of = "16.top_pop_figure.45uM.tab";
15+
else if (a[2] == "25A" || a[2] == "25B" || a[2] == "25C")
16+
of = "16.top_pop_figure.75uM.tab";
17+
else if (a[2] == "50A" || a[2] == "50B" || a[2] == "50C")
18+
of = "16.top_pop_figure.150uM.tab";
19+
else if (a[2] == "75A" || a[2] == "75B" || a[2] == "75C")
20+
of = "16.top_pop_figure.225uM.tab";
21+
else {
22+
files[1] = "16.top_pop_figure.15uM.tab";
23+
files[2] = "16.top_pop_figure.45uM.tab";
24+
files[3] = "16.top_pop_figure.75uM.tab";
25+
files[4] = "16.top_pop_figure.150uM.tab";
26+
files[5] = "16.top_pop_figure.225uM.tab";
27+
for (f in files) {
28+
of = files[f]
29+
for (j = 1; j <= line; ++j) {
30+
printf("%s\t", data[j, i]) >> of;
31+
}
32+
printf("\n") >> of;
33+
}
34+
continue;
35+
}
36+
for (j = 1; j <= line; ++j) {
37+
printf("%s\t", data[j, i]) >> of;
38+
}
39+
printf("\n") >> of;
40+
}
41+
}

‎16.top_pop_figure.gplt‎

Lines changed: 50 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,50 @@
1+
set terminal pdf color solid enhanced
2+
set output '16.top_pop_figure.pdf'
3+
#set datatfile separator "\t"
4+
set border 3 front linetype -1 linewidth 1.000
5+
set boxwidth 0.8 absolute
6+
set style fill solid 1.00 noborder
7+
set grid nopolar
8+
set grid noxtics nomxtics ytics nomytics noztics nomztics \
9+
nox2tics nomx2tics noy2tics nomy2tics nocbtics nomcbtics
10+
set grid layerdefault linetype 0 linewidth 1.000, linetype 0 linewidth 1.000
11+
set key bmargin center horizontal Left reverse noenhanced autotitles columnhead nobox
12+
#unset key
13+
set style histogram rowstacked title offset character 2, 0.25, 0
14+
set datafile missing '-'
15+
set style data histograms
16+
#set xtics border in scale 0,0 nomirror rotate by -45 offset character 0, 0, 0 autojustify
17+
set xtics border in scale 0,0 nomirror rotate by -90 offset character 0, 0, 0 autojustify
18+
set xtics norangelimit font ",8"
19+
set xtics ()
20+
set ytics border in scale 0,0 mirror norotate offset character 0, 0, 0 autojustify
21+
set ztics border in scale 0,0 nomirror norotate offset character 0, 0, 0 autojustify
22+
set cbtics border in scale 0,0 mirror norotate offset character 0, 0, 0 autojustify
23+
set rtics axis in scale 0,0 nomirror norotate offset character 0, 0, 0 autojustify
24+
set xlabel "[O_2] in enrichment"
25+
set xlabel offset character 0, -2, 0 font "" textcolor lt -1 norotate
26+
set ylabel "Fraction of population"
27+
set yrange [ 0 : 100 ] noreverse nowriteback
28+
29+
# setup color blind palette
30+
#set palette model RGB defined (0 "#88CCEE", 1 "#44AA99", 1 "#117733", 2 "#999933", 2 "#DDCC77", 3 "#CC6677", 4 "#882255", 5 "#AA4499", 6 "grey")
31+
set linetype 1 lc rgb "#88CCEE" lw 1
32+
set linetype 2 lc rgb "#44AA99" lw 1
33+
set linetype 3 lc rgb "#117733" lw 1
34+
set linetype 4 lc rgb "#999933" lw 1
35+
set linetype 5 lc rgb "#DDCC77" lw 1
36+
set linetype 6 lc rgb "#CC6677" lw 1
37+
set linetype 7 lc rgb "#882255" lw 1
38+
set linetype 8 lc rgb "#AA4499" lw 1
39+
set linetype 9 lc rgb "gray80" lw 1
40+
set linetype cycle 9
41+
42+
#unset colorbox
43+
#unset key
44+
45+
plot \
46+
newhistogram "15 {/Symbol m}M" lt 2, '16.top_pop_figure.15uM.tab' using 2:xtic(1) t col, '' u 3:xtic(1) t col, 4:xtic(1) t col, 5:xtic(1) t col, 6:xtic(1) t col, 7:xtic(1) t col, 8:xtic(1) t col, \
47+
newhistogram "45 {/Symbol m}M" lt 2, '16.top_pop_figure.45uM.tab' using 2:xtic(1) t col, '' u 3:xtic(1) t col, 4:xtic(1) t col, 5:xtic(1) t col, 6:xtic(1) t col, 7:xtic(1) t col, 8:xtic(1) t col, \
48+
newhistogram "75 {/Symbol m}M" lt 2, '16.top_pop_figure.75uM.tab' using 2:xtic(1) t col, '' u 3:xtic(1) t col, 4:xtic(1) t col, 5:xtic(1) t col, 6:xtic(1) t col, 7:xtic(1) t col, 8:xtic(1) t col, \
49+
newhistogram "150 {/Symbol m}M" lt 2, '16.top_pop_figure.150uM.tab' using 2:xtic(1) t col, '' u 3:xtic(1) t col, 4:xtic(1) t col, 5:xtic(1) t col, 6:xtic(1) t col, 7:xtic(1) t col, 8:xtic(1) t col, \
50+
newhistogram "225 {/Symbol m}M" lt 2, '16.top_pop_figure.225uM.tab' using 2:xtic(1) t col, '' u 3:xtic(1) t col, 4:xtic(1) t col, 5:xtic(1) t col, 6:xtic(1) t col, 7:xtic(1) t col, 8:xtic(1) t col

‎16.top_pop_figure.pdf‎

Whitespace-only changes.

‎16.top_pop_figure.tab‎

Lines changed: 8 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,8 @@
1+
OTUid O2_15A_10 O2_15B_10 O2_15C_10 O2_25A_10 O2_25B_10 O2_25C_10 O2_50A_10 O2_50B_10 O2_50C_10 O2_5A_10 O2_5B_10 O2_5C_10 O2_75A_10 O2_75B_10 O2_75C_10 O2_15A_15 O2_15B_15 O2_15C_15 O2_25A_15 O2_25B_15 O2_25C_15 O2_50A_15 O2_50B_15 O2_50C_15 O2_5A_15 O2_5B_15 O2_5C_15 O2_75A_15 O2_75B_15 O2_75C_15
2+
OTU_115 0.000203533338760889 8.77943855490441e-05 0.0017664136556018 0.00441885715728606 0.00575535043165128 0.000253292806484296 0.000246404041026273 0.0198138357959423 0.00175216866777733 0.0131554015645725 0.000460809268849293 0.000962592582689377 0.000157996603073034 0.000194225042082092 0.00236502193204099 0.00193254497768371 0.000241266164833044 0.000368091707419677 0.0525107917600131 0.00600555886440341 0.013517359134889 0.00159716060337178 0.020653557787523 0.000475737392959087 0.000780284803953443 0.00128609844499006 0.0250296208530806 0.000433254705627497 0 0.00014808233377758
3+
OTU_4 0.00199462671985671 0.00215096244595158 0.102032250166148 0.0216145241522107 0.00567672542575441 0.00405268490374873 0.00221763636923646 0.00449595589896525 0.00206813350950767 0.469507078472241 0.525278679891161 0.0130217385491591 0.00231070031994312 0.00183434761966421 0.00238523579470801 0.00257672663691161 0.00241266164833044 0.00347057895567124 0.0274301950713076 0.0130037720865595 0.0120470498956555 0.00301685891748004 0.00841703211203848 0.00445463195225327 0.108654658950517 0.241786507658132 0.199200236966825 0.00226255235161026 0 0.00133274100399822
4+
OTU_3 0.0491940079785069 0.115427668400606 0.517769072020707 0.0199479837386057 0.00929347569701068 0.183320668693009 0.0522992577078264 0.00718449234106005 0.00301602803469868 0.0212489510868094 0.0595979987711753 0.0201877055536244 0.00315993206146068 0.0637273943631577 0.00430555274807463 0.0346017576956702 0.585118702953098 0.328916232844297 0.617452598218677 0.0250645225332539 0.4314171883893 0.143655723158829 0.0280096194652709 0.0235706253784275 0.381884387801548 0.233309949725243 0.0313240521327014 0.0331680546863717 0.0155936094900985 0.0103657633644306
5+
OTU_1 0.488764959700399 0.272821053093655 0.00834236944279268 0.00583289144761761 0.564684792351359 0.206876899696049 0.190439523208181 0.501830012200081 0.161659102659849 0.00646942587228974 0.0073071184060388 0.153132436696168 0.656140142986926 0.419806638180327 0.714741970043056 0.391846500713201 0.0123045744064852 0.0226113477414945 0.102180208731763 0.604973198332341 0.255169797002466 0.37435669920142 0.517682840571509 0.710319176541822 0.0121594381949412 0.0115748860049106 0.0301392180094787 0.490588745005536 0.510044963168468 0.450614541685177
6+
OTU_239 0.00122120003256533 0.00173393911459362 0.000787015985169121 0.0012877812286948 0.00110075008255626 0.00151975683890578 0.0874426340591986 0.0672585965387917 0.00318837249382432 0.00108274909996481 0.00111910822434828 0.00264712960239579 0.00152071730457795 0.306422374724848 0.00119261789735401 0.00156444117241062 0.00101331789229878 0.00152495135931009 0.000764985519916944 0.00198530871550526 0.0011857332574464 0.0375332741792369 0.013085302022917 0.00086497707810743 0.000520189869302295 0.00239681982929966 0.00192535545023697 0.208491792230299 0.169185879651775 0.00251739967421887
7+
OTU_5 0.00588211349018969 0.00796734048857576 0.00363776277589283 0.00653990859278337 0.00575535043165128 0.00698665991219183 0.578649089845073 0.244453481541729 0.165307060378009 0.00476409603984517 0.00441060300184324 0.00981309660686115 0.247758423193901 0.149941732487375 0.196903236239413 0.00630377766530162 0.00414977803512835 0.00625755902613451 0.00289601661111415 0.00565812983918999 0.00426863972680706 0.225022182786158 0.135238364690904 0.00320041518899749 0.00435659015540672 0.00485209867882614 0.00651658767772512 0.125355028161556 0.220989189706304 0.424255886272768
8+
OTU_2 0.410893104290483 0.543908167072716 0.234863059218581 0.566371234502437 0.349928451244634 0.515873015873016 0.0210675455077463 0.0155438073290859 0.0178663755960246 0.349809165471131 0.303256385499868 0.735420733174684 0.0142591934273413 0.0183003150761794 0.0459056821167957 0.46643353425666 0.342549700829956 0.592049219119735 0.159881973662641 0.318145721659718 0.248624549421362 0.033717834960071 0.0326778893761494 0.0524608597872156 0.353208921256259 0.36630422074126 0.663803317535545 0.0613296105521591 0.0315698842437578 0.031393454760847

0 commit comments

Comments
 (0)