目前研究维生素D与血糖血脂等代谢指标的关系，想控制年龄和性别，使用偏相关分析无法控制年龄，而多重线性回归需要因果关系，应该怎么选择统计方法呢？