Solubility prediction of organic ionic compounds in both aqueous and organic solvents is important for understanding and optimizing lithographic performances. In this study, we proposed computational methods to predict solubility of organic ionic compounds. To compare the predicted solubility with the experimental one, we applied a multiple linear regression model by changing a set of explanatory variables. We conclude that the variables of solvation free energies of cation-anion pair, cation and anion, which are Delta G(AB)degrees, Delta G(A)degrees and Delta G(B)degrees respectively, will be sufficient to describe the relationship between the predicted and experimental solubility values. We expect that the more accurate empirical model for quantitative prediction of solubility of organic ionic compounds by expanding these regression models and further optimizing the parameters based on larger set of experimental values will be reserved