如何用R语言做半参数回归，是用mgcv包还是其他的，请提供一个具体的实例（含数据和整个运行代码，并且有适当的文字叙述和对结果的解释）
三天有效，过期关闭主题。