The best fit model found was log(P(Tuberculosis)/P(sarcoidosis))  = 1.214 + 8.884 * IL1 + 0.075 * IFN + 0.32 * IL12p70 + 0.039 * TNF − 0.097 * IL5 − 0.015 * IL8 − 0.73 * IL13 − 0.606 * IL10 − 2.404 * IL2 − 4.231 * IL4. The gene discussed is IL13; the disease is sarcoidosis.