广义可加模型（即GAM）中偏残差图究竟表示什么意思？

ywh19860616 发表于 2013-4-12 14:15
您好，想请教下，在估计glm模型时通常需要先
选择link function，如果我的数据只适合于normal，
而模型 ...

谢谢您的问题。我只提供我的一些理解，不一定正确。
这要看你的数据量的大小。如果您的数据不多（比如只有不到30个数据点），出现非正态是非常正常的，现在大量的研究，比如方差分析，即便不是正态性照样用。
如果您的数据量较大，比如超过300个点，而残差不满足正态性，那么你就需要注意了。可能模型的选用存在问题。
对于30-300之间的样本量，很难说使用的模型正确与否了。我的经验。
一般你使用glm，估计别人也不会对你的正态性问题进行提问，因为如果glm都不能使用的话，普通的线性回归或者非线性回归就更不能使用了。

peijianshi 发表于 2013-4-13 09:58
谢谢您的问题。我只提供我的一些理解，不一定正确。
这要看你的数据量的大小。如果您的数据不多（比如只 ...

谢谢您的回复，答案非常精彩。
我的样本量是60多个，的确是这样，我又利用了不同的正态性检验方法进行
检验，所得到的结果也存在差异，所以我想或许可以不非常关注这个问题。

在glm模型中常用的模型拟合优度指标有deviance，pearson等，在选择link function
为normal时，好像两者是相等吧？

ywh19860616 发表于 2013-4-13 12:07
谢谢您的回复，答案非常精彩。
我的样本量是60多个，的确是这样，我又利用了不同的正态性检验方法进行
...

呵呵，我对这些指标没有做过研究。以前都是根据预测值来计算决定系数R^2或者AIC，澳大利亚一位朋友说使用R^2一点意义都没有，我也不知道他说的在理不。一般我使用R^2，
cor(y.obervation, y.prediction)^2就大致等于R^2了，越接近一越好。

peijianshi 发表于 2013-4-14 13:53
呵呵，我对这些指标没有做过研究。以前都是根据预测值来计算决定系数R^2或者AIC，澳大利亚一位朋友说使用 ...

好的，非常感谢您的解答

peijianshi 发表于 2012-3-27 18:57
我已经清楚啦。偏残差是指因变量减去期望值后再减去其他变量的s(x)后的残差，对于简单广义可加模型而言（高 ...

楼主看来是无师自通啊，牛！
我也接触R用来做GAM，虽然很快能照着书自己做一遍，可是里面很多过程和参数的选择都还是一头雾水，但部分还是用软件默认的。
我有个问题想问下大家
c<-gam(y~s(x1)+s(x2)+s(x3,x4))
在这个方程里面，s(x3,x4)有实际的物理意义吗？是通常意义上的交互作用吗？
用plot(c)这个函数画出来的图感觉像个等值线图，不知道怎么读那个图，楼主不知道能不能赐教一二，不甚感激啦

zhouyuanshen 发表于 2013-5-14 10:41
楼主看来是无师自通啊，牛！
我也接触R用来做GAM，虽然很快能照着书自己做一遍，可是里面很多过程和参数 ...

关于交互，我个人认为，使用起来要慎重。最初的交互作用更多来源于医学统计上的药物之间的相互影响，但是以后出现的滥用，许多风马牛不相及的事情，有些人总喜欢搞些交互作用，许多是伪交互。交互只能使分析结果更为混乱。处于在逻辑上具有明显的可能互作，否则尽量避免交互分析。
广义可加模型缺点在于对交互作用的分析。你可以通过Hastie & Tibshirani (1986)的原始论文后的一些评议中看出。作者则主张可以使用x1*x2这类去试试。
te(x1, x2)似乎要好于s(x1*x2)。具体含义需要自己琢磨。广义可加模型的使用最大的问题就是移动参数的选取带有经验的主观性，同时对是否超出了拟合的合理范围，比如造成过度拟合的判断。非常难说。这极大制约了广义可加模型的使用。呵呵。

前辈，您好，我刚刚接触gam不久，想问一下，在mgcv拟合得到模型之后，如何得到（R代码）非参数项的lgRR的曲线图，直接使用plot得到的只有s(x)与x的关系图，怎么得到lgRR与x的关系图，找了好久没有找到资料，还望解答。

楼上大神，想清问您有没有非参数GARCH的程序？

前辈，我的毕业论文用到非参数GARCH这一块，可是对于程序我琢磨了好久没弄出来，请您赐教，多谢多谢！！！！非常感谢！！！！！

还想再请问前辈，lo（）是什么函数啊？在哪个程序包里呢？