The present invention provides a set of in vitro and in silico methodologies for predicting in vivo pharmacokinetics and pharmacodynamics of multiple components; the methodologies comprise mathematical models for solving multiple unknowns which are linearly independent and/or interacting with each other. The present invention can be applied to develop phytomedicines which contain multiple active ingredients without prior identification, isolation and purification of these components.