def get_latex_documentbody(fs):
"""
This is obsolete because I am now using pure R output.
The latex documentbody should have a bunch of tikz pieces in it.
Each tikz piece should have been generated from R.
"""
Q_mut, Q_sels = get_qmut_qsels(fs)
# compute the statistics
ER_ratios, NSR_ratios, ER_NSR_ratios = get_statistic_ratios(Q_mut, Q_sels)
M = zip(*(ER_ratios, NSR_ratios, ER_NSR_ratios))
column_headers = ('ER.ratio', 'NSR.ratio', 'ER.times.NSR.ratio')
table_string = RUtil.get_table_string(M, column_headers)
nsels = len(Q_sels)
# define the R scripts
scripts = []
for name in column_headers:
scripts.append(get_r_tikz_script(nsels, name))
# get the tikz codes from R, for each histogram
retcode, r_out, r_err, tikz_code_list = RUtil.run_plotter_multiple_scripts(
table_string, scripts, 'tikz',
width=3, height=2)
if retcode:
raise RUtil.RError(r_err)
#
# show some timings
print 'R did not fail, but here is its stderr:'
print r_err
#
# write the latex code
out = StringIO()
#print >> out, '\\pagestyle{empty}'
for tikz_code in tikz_code_list:
print >> out, tikz_code
# return the latex code, consisting mainly of a bunch of tikz plots
return out.getvalue()
def get_latex_documentbody(fs):
"""
This is obsolete because I am now using pure R output.
The latex documentbody should have a bunch of tikz pieces in it.
Each tikz piece should have been generated from R.
"""
Q_mut, Q_sels = get_qmut_qsels(fs)
# compute the statistics
ER_ratios, NSR_ratios, ER_NSR_ratios = get_statistic_ratios(Q_mut, Q_sels)
M = zip(*(ER_ratios, NSR_ratios, ER_NSR_ratios))
column_headers = ('ER.ratio', 'NSR.ratio', 'ER.times.NSR.ratio')
table_string = RUtil.get_table_string(M, column_headers)
nsels = len(Q_sels)
# define the R scripts
scripts = []
for name in column_headers:
scripts.append(get_r_tikz_script(nsels, name))
# get the tikz codes from R, for each histogram
retcode, r_out, r_err, tikz_code_list = RUtil.run_plotter_multiple_scripts(
table_string, scripts, 'tikz', width=3, height=2)
if retcode:
raise RUtil.RError(r_err)
#
# show some timings
print 'R did not fail, but here is its stderr:'
print r_err
#
# write the latex code
out = StringIO()
#print >> out, '\\pagestyle{empty}'
for tikz_code in tikz_code_list:
print >> out, tikz_code
# return the latex code, consisting mainly of a bunch of tikz plots
return out.getvalue()
