作者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)