作者TheHell (勒黑唷)
看板Statistics
标题[问题] 请问有关SAS与R中的GAM模式建构
时间Wed May 2 14:02:09 2007
不好意思
我最近在尝试利用R与SAS进行GAM的运算
但同时都遇到了一些瓶颈
并且在网路上没有搜寻到相关的范例可供参考
想请问不知版上有没有有进行过类似运算的朋友可以帮我看一下
在R的部分
我使用的是gam package
(因为mgcv不能用loess)
我总共有四个变数
分别为Recruitment <---y变数
SST SOI 及SSB <---x变数
我尝试进行了28种排列组合
其中包括尝试在不同x变数上使用loess或是cubic spline
或在model中加入interaction的考量
最後用ANOVA table去比较这28个model
然後发现其中model4的AIC是最小的
我想知道每个变数的影响力有多大
所以便在R中打了summary(model4)
不过我却发现interaction在R中
他不会帮我去做显着的测定
瓶颈一:在R中我要如何得知interaction是否有显着影响力?
在SAS的部分
主要是因为R里得不到我想知道的interaction影响力
所以想说用SAS做也许output会比较详细
但是我却在第一步model输入的时候就於无法成功
我仿造我在R中得到的AIC最小的model
尝试在SAS建构
於是我输入以下指令
proc gam data=test;
model recu=loess(sst,ssb) spline(soi) param(sst*soi) param(ssb*soi)
param(sst*soi*ssb);
output out=PredGAM;
run;
不过也是在model输入部分error一堆
瓶颈二:请问我model的建构问题在哪里?
SAS似乎是不能在loess中放入两个变数
不知这是他软体本身的限制
还是可以经由其他的设定?
不好意思 落落长一大篇
我目前在网上找到的GAM说明或范例指令相当少
所以想说来这边请教
谢谢大家耐心看完
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.112.69.148
※ 编辑: TheHell 来自: 140.112.69.148 (05/02 14:18)
※ 编辑: TheHell 来自: 140.112.69.148 (05/02 14:20)