A newly synthesized bis thiosemicarbazone ligand, (2Z,2′Z)-2,2′-((4S,5R)-4,5,6-trihydroxyhexane-1,2-diylidene)bis(N-phenylhydrazinecarbothioamide), was used to make a complex with Cu2+, Ni2+, Co2+ and Fe3+ for their simultaneous spectrophotometric determination using chemometric methods. By Job’s method, the ratio of metal to ligand in Ni2+ was found to be 1:2, whereas it was 1:4 for the others. The effect of pH on the sensitivity and selectivity of the formed complexes was studied according to the net analyte signal (NAS). Under optimum conditions, the calibration graphs were linear in the ranges of 0.10–3.83, 0.20–3.83, 0.23–5.23 and 0.32–8.12mgL−1 with the detection limits of 2, 3, 4 and 10μgL−1 for Cu2+, Co2+, Ni2+ and Fe3+ respectively. The OSC-PLS1 for Cu2+ and Ni2+, the PLS1 for Co2+ and the PC-FFANN for Fe3+ were selected as the best models. The selected models were successfully applied for the simultaneous determination of elements in some foodstuffs and vegetables.