public void GenerateVennDiagrams()
{
var methodName = MethodBase.GetCurrentMethod().Name;
TestUtils.ShowStarting(methodName);
// DIA
const string dir = @"H:\Research\EDRN\Ic\DIA_Replicate";
const string rep1 = dir + @"\EDRN_Serum_07_DIA_1_01_13Nov13_Samwise_13-07-28_IcTda.tsv";
//const string rep2 = dir + @"\EDRN_Serum_07_DIA_1_02_13Nov13_Samwise_13-07-28_IcTda.tsv";
//const string rep3 = dir + @"\EDRN_Serum_07_DIA_1_03_13Nov13_Samwise_13-07-28_IcTda.tsv";
const string rep4 = dir + @"\EDRN_Serum_07_DIA_1_04_13Nov13_Samwise_13-07-28_IcTda.tsv";
//const string rep5 = dir + @"\EDRN_Serum_07_DIA_1_05_18Nov13_Samwise_13-07-28_IcTda.tsv";
const string resultPath1 = rep1;
const string resultPath2 = rep4;
var result1 = new TsvFileParser(resultPath1);
var result2 = new TsvFileParser(resultPath2);
const double pepQValueThreshold = 0.01;
var vennDiagram = new VennDiagram <string>(result1.GetPeptidesAboveQValueThreshold(pepQValueThreshold),
result2.GetPeptidesAboveQValueThreshold(pepQValueThreshold));
Console.WriteLine("{0}\t{1}\t{2}",
vennDiagram.Set1Only.Count, // + vennDiagram.Intersection.Count,
vennDiagram.Intersection.Count,
vennDiagram.Set2Only.Count //+ vennDiagram.Intersection.Count
);
}
public void GenerateVennDiagrams()
{
var methodName = MethodBase.GetCurrentMethod().Name;
Utils.ShowStarting(methodName);
// DIA
const string dir = @"H:\Research\DDAPlus\NTT2";
//const string dda1 = dir + @"\20140701_yeast_DDA_01_IcTda.tsv";
//const string dda2 = dir + @"\20140701_yeast_DDA_02_2_IcTda.tsv";
const string ddaPlus1 = dir + @"\20140701_yeast_DDAp_binCharge_01_IcTda.tsv";
const string ddaPlus2 = dir + @"\20140701_yeast_DDAp_binCharge_02_IcTda.tsv";
const string resultPath1 = ddaPlus1;
const string resultPath2 = ddaPlus2;
if (!File.Exists(resultPath1))
{
Assert.Ignore(@"Skipping test {0} since file not found: {1}", methodName, resultPath1);
}
if (!File.Exists(resultPath2))
{
Assert.Ignore(@"Skipping test {0} since file not found: {1}", methodName, resultPath2);
}
var result1 = new TsvFileParser(resultPath1);
var result2 = new TsvFileParser(resultPath2);
const double pepQValueThreshold = 0.01;
var vennDiagram = new VennDiagram <string>(result1.GetPeptidesAboveQValueThreshold(pepQValueThreshold),
result2.GetPeptidesAboveQValueThreshold(pepQValueThreshold));
Console.WriteLine("{0}\t{1}\t{2}",
vennDiagram.Set1Only.Count, // + vennDiagram.Intersection.Count,
vennDiagram.Intersection.Count,
vennDiagram.Set2Only.Count //+ vennDiagram.Intersection.Count
);
}