internal void CreateNoSwitchablesHlaAssignment()
{
BestHlaAssignmentSoFar = BestSoFar <double, TrueCollection> .GetInstance(SpecialFunctions.DoubleGreaterThan);
TrueCollection trueCollection = TrueCollection.GetInstance(QmrrModelMissingAssignment.KnownHlaSet);
double scoreOnNoSwitchables = QmrrModelMissingAssignment.LogLikelihoodOfCompleteModelConditionedOnKnownHlas(trueCollection);
BestHlaAssignmentSoFar.Compare(scoreOnNoSwitchables, trueCollection);
}
private void ReportPerHla(OptimizationParameterList qmrrParams, Dictionary <string, BestSoFar <double, TrueCollection> > peptideToBestHlaAssignmentSoFar,
string directory, string name)
{
string fileName = string.Format(@"{0}\NoisyOr.PeptideHlaProbability.{1}.new.txt", directory, name);
using (StreamWriter output = File.CreateText(fileName))
{
output.WriteLine(SpecialFunctions.CreateTabString("Peptide", "HLA", "LogOdds", "Probability"));
foreach (QmrrPartialModel qmrrPartialModel in QmrrPartialModelCollection)
{
QmrrModelMissingAssignment aQmrrModelMissingAssignment = QmrrModelMissingAssignment.GetInstance(ModelLikelihoodFactories, qmrrPartialModel, qmrrParams);
BestSoFar <double, TrueCollection> bestHlaAssignment = peptideToBestHlaAssignmentSoFar[qmrrPartialModel.Peptide];
Set <Hla> trueCollectionFullAsSet = new Set <Hla>(bestHlaAssignment.Champ);
double loglikelihoodFull = bestHlaAssignment.ChampsScore;
foreach (Hla hla in trueCollectionFullAsSet)
{
bool known = qmrrPartialModel.KnownHlaSet.Contains(hla);
if (!known)
{
Set <Hla> allLessOne = trueCollectionFullAsSet.SubtractElement(hla);
TrueCollection trueCollectionWithout = TrueCollection.GetInstance(allLessOne);
double loglikelihoodWithout = aQmrrModelMissingAssignment.LogLikelihoodOfCompleteModelConditionedOnKnownHlas(trueCollectionWithout);
double logOdds = loglikelihoodFull - loglikelihoodWithout;
double probabilityFull = Math.Exp(loglikelihoodFull);
double probabilityWithout = Math.Exp(loglikelihoodWithout);
double probability = probabilityFull / (probabilityFull + probabilityWithout);
Debug.Assert(logOdds >= 0); // real assert
output.WriteLine(SpecialFunctions.CreateTabString(qmrrPartialModel.Peptide, hla, logOdds, probability));
}
}
}
}
}