def setUp(self):
ImportedDataTest.setUp(self)
self.ex = flow.ChannelStatisticOp(name="Y",
channel="Y2-A",
by=['Well', 'Dox'],
function=flow.geom_mean).apply(
self.ex)
self.ex = flow.ChannelStatisticOp(name="Y",
channel="Y2-A",
by=['Well', 'Dox'],
function=flow.geom_sd_range).apply(
self.ex)
self.ex = flow.ChannelStatisticOp(name="V",
channel="V2-A",
by=['Well', 'Dox'],
function=flow.geom_mean).apply(
self.ex)
self.ex = flow.ChannelStatisticOp(name="V",
channel="V2-A",
by=['Well', 'Dox'],
function=flow.geom_sd_range).apply(
self.ex)
self.view = flow.Stats2DView(xstatistic=("Y", "geom_mean"),
x_error_statistic=("Y", "geom_sd_range"),
ystatistic=("V", "geom_mean"),
y_error_statistic=("V", "geom_sd_range"),
variable="Dox",
huefacet="Well")
def testStats2D(self):
import numpy as np
flow.Stats2DView(name="Stats2D",
by="Dox",
xchannel="V2-A",
xfunction=np.mean,
ychannel="Y2-A",
yfunction=np.mean).plot(self.ex)
def testBadYErrorStatistic(self):
self.ex = flow.ChannelStatisticOp(name = "V_bad",
channel = "V2-A",
by = ['Well'],
function = flow.geom_sd_range).apply(self.ex)
self.view = flow.Stats2DView(xstatistic = ("Y", "geom_mean"),
x_error_statistic = ("Y", "geom_sd_range"),
ystatistic = ("V", "geom_mean"),
y_error_statistic = ("V_bad", "geom_sd_range"),
variable = "Dox",
huefacet = "Well")
self.assertRaises(util.CytoflowViewError, self.view.plot, self.ex)
def testStats2D(self):
import numpy as np
ex2 = flow.ChannelStatisticOp(name="StatsV",
by=["Dox"],
channel="V2-A",
function=np.mean).apply(self.ex)
ex3 = flow.ChannelStatisticOp(name="StatsY",
by=["Dox"],
channel="Y2-A",
function=np.mean).apply(ex2)
flow.Stats2DView(xstatistic=("StatsV", "mean"),
ystatistic=("StatsY", "mean"),
variable="Dox").plot(ex3)
plt.vlines(x, y_lo, y_hi, color=color, **kwargs)
if __name__ == '__main__':
import cytoflow as flow
tube1 = flow.Tube(file='../../cytoflow/tests/data/Plate01/RFP_Well_A3.fcs',
conditions={"Dox": 10.0})
tube2 = flow.Tube(file='../../cytoflow/tests/data/Plate01/CFP_Well_A4.fcs',
conditions={"Dox": 1.0})
ex = flow.ImportOp(conditions={"Dox": "float"}, tubes=[tube1, tube2])
thresh = flow.ThresholdOp()
thresh.name = "Y2-A+"
thresh.channel = 'Y2-A'
thresh.threshold = 2005.0
ex2 = thresh.apply(ex)
s = flow.Stats2DView()
s.variable = "Dox"
s.xchannel = "V2-A"
s.xfunction = np.mean
s.ychannel = "Y2-A"
s.yfunction = np.mean
s.huefacet = "Y2-A+"
plt.ioff()
s.plot(ex2)
plt.show()
