本文整理汇总了C#中Workspace.filterByMinimumCoveringSet方法的典型用法代码示例。如果您正苦于以下问题:C# Workspace.filterByMinimumCoveringSet方法的具体用法?C# Workspace.filterByMinimumCoveringSet怎么用?C# Workspace.filterByMinimumCoveringSet使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类Workspace
的用法示例。
在下文中一共展示了Workspace.filterByMinimumCoveringSet方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的C#代码示例。
示例1: Main
//.........这里部分代码省略.........
start = DateTime.Now.Ticks;
Console.WriteLine( "Assembling protein groups..." );
ws.assembleProteinGroups();
Console.WriteLine( "\nFinished assembling protein groups; " + new TimeSpan( DateTime.Now.Ticks - start ).TotalSeconds + " seconds elapsed." ); { ws.validate( rtConfig.MaxFDR, rtConfig.MinDistinctPeptides, rtConfig.MaxAmbiguousIds ); }
start = DateTime.Now.Ticks;
Console.WriteLine( "Assembling peptide groups..." );
ws.assemblePeptideGroups();
Console.WriteLine( "\nFinished assembling peptide groups; " + new TimeSpan( DateTime.Now.Ticks - start ).TotalSeconds + " seconds elapsed." ); { ws.validate( rtConfig.MaxFDR, rtConfig.MinDistinctPeptides, rtConfig.MaxAmbiguousIds ); }
start = DateTime.Now.Ticks;
Console.WriteLine( "Assembling clusters..." );
ws.assembleClusters();
Console.WriteLine( "\nFinished assembling clusters; " + new TimeSpan( DateTime.Now.Ticks - start ).TotalSeconds + " seconds elapsed." ); { ws.validate( rtConfig.MaxFDR, rtConfig.MinDistinctPeptides, rtConfig.MaxAmbiguousIds ); }
start = DateTime.Now.Ticks;
// Determine the minimum number of protein clusters needed
// to explain all the peptides.
Console.WriteLine( "Assembling minimum covering set for clusters..." );
int clusterCount = 0;
foreach( ClusterInfo c in ws.clusters )
{
++clusterCount;
Console.Write( "Assembling minimum covering set for cluster " + clusterCount + " of " + ws.clusters.Count + " (" + c.proteinGroups.Count + " protein groups, " + c.results.Count + " results) \r" );
ws.assembleMinimumCoveringSet( c );
}
Console.WriteLine( "\nFinished assembling minimum covering sets; " + new TimeSpan( DateTime.Now.Ticks - start ).TotalSeconds + " seconds elapsed." ); { ws.validate( rtConfig.MaxFDR, rtConfig.MinDistinctPeptides, rtConfig.MaxAmbiguousIds ); }
if( rtConfig.MinAdditionalPeptides > 0 )
{
start = DateTime.Now.Ticks;
Console.WriteLine( "Filtering workspace by minimum covering set..." );
ws.filterByMinimumCoveringSet( rtConfig.MinAdditionalPeptides );
Console.WriteLine( "\nFinished filtering by minimum covering set; " + new TimeSpan( DateTime.Now.Ticks - start ).TotalSeconds + " seconds elapsed." ); { ws.validate( rtConfig.MaxFDR, rtConfig.MinDistinctPeptides, rtConfig.MaxAmbiguousIds ); }
}
start = DateTime.Now.Ticks;
Console.Write( "Assembling source groups..." );
ws.assembleSourceGroups();
Console.WriteLine( " finished; " + new TimeSpan( DateTime.Now.Ticks - start ).TotalSeconds + " seconds elapsed." );
#endregion
#region Post-filtering validation
// Check to see if the protein, peptides, sources, and the clusters have
// all required fields and satisfy the user-defined criteria.
try
{
ws.validate( rtConfig.MaxFDR, rtConfig.MinDistinctPeptides, rtConfig.MaxAmbiguousIds );
} catch( Exception e )
{
Console.Error.WriteLine( "Error while validating workspace: " + e.Message );
return;
}
#endregion
#region SNP annotation stuff (TODO: move this into a Workspace function)
if( rtConfig.FlagUnknownMods.Length > 0 )
{
string[] toks = rtConfig.FlagUnknownMods.Split( new char[] { ' ' } );
if( toks.Length % 2 == 0 )
{
ws.knownModMasses = new float[toks.Length / 2];
ws.knownModResidues = new char[toks.Length / 2];
int index = 0;
for( int i = 0; i < toks.Length - 1; i = i + 2 )