Elucidation of the systems biology foundation underlying the effect of Fangji, which are multi-herbal traditional Chinese medicine (TCM) formulas, is one of the major aims in the field. The numerous bioactive ingredients of a Fangji deal with the multiple targets of a complex disease, which is influenced by a number of genes and their interactions with the environment. Genome-wide association study (GWAS) is an unbiased approach for dissecting the genetic variants underlying complex diseases and individual response to a given treatment. GWAS has great potential for the study of systems biology from the point of view of genomics, but the capacity using current analysis models is largely handicapped, as evidenced by missing heritability. Recent development of a full genetic model, in which gene-gene interactions (dominance and epistasis) and gene-environment interactions are all considered, has addressed these problems. This approach has been demonstrated to substantially increase model power, remarkably improving the detection of association of GWAS and the construction of the molecular architecture. This analysis does not require a very large sample size, which is often difficult to meet for a GWAS of treatment response. Furthermore, this analysis can integrate other omic information and allow for variations of Fangji, which is very promising for Fangjiomic study and detection of the sophisticated molecular architecture of the function of Fangji, as well as for the delineation of the systems biology of personalized medicine in TCM in an unbiased and comprehensive manner.