如题，举个简单的例子，在可加模型y=x1*b1+x2*b2+s(z1)+s(z2)+v  中，想得到 v=My中的M，即残差生成矩阵，在R中mgcv包或者gam可以实现么？或者可以通过编程实现么？ 求各位指导。