Results: The diagnostic power of the binary logistic regression model incorporating skeleton-related symptoms, anemia, the SUVmax of bone lesions, bone changes, the location of bone lesions, ALP, LDH, CEA, and CA19-9 was significantly higher than that of the model using only clinical factors (p = 0.008). The gene discussed is CEACAM5; the disease is anemia (phenotype).