[程式] sas跑非線性迴歸

看板Statistics作者 ( )時間15年前 (2009/10/04 03:18), 編輯推噓1(108)
留言9則, 2人參與, 最新討論串1/1
------------------------------------------------------------------------ [軟體程式類別]:SAS [程式問題]: 非線性迴歸參數估計 [軟體熟悉度]: 中(3個月到1年) [問題敘述]: 大家好 我現在想要跑一個非線性迴歸 然後估計他的參數 只做一次試驗的話 參數的估計值會變動很大 (每次估計出來的值離預設值忽近忽遠) 所以我希望能重複試驗 然後算參數的平均 看看這樣會不會比較穩定一點 但是我不知道要怎麼樣做才能把我的"程式"整個重複 [就是從抽樣到計算出參數的部分重複做100次...] 希望有高手能幫忙一下 謝謝 [程式範例]: data a; do j = 1 to 100; keep x u y; do i = 1 to 100; x =rand('uniform'); u = rand('normal' ,0,1); y = 0.5 +0.8/(x-2) +u; end; output; ^^^^^^^ 這裡少打 end; run; proc nlin best =30 ; parms a = 0 b = 0 c =0 ; model y = a + b/(x-c); output out =p p=predict; run; ----------------------------------------------------------------------------- -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 61.217.235.223

10/04 12:04, , 1F
妳的a最後有多少觀察值,妳有檢驗嗎?
10/04 12:04, 1F

10/05 00:01, , 2F
有一萬 可是我是希望每一百筆資料可以跑一次回歸
10/05 00:01, 2F

10/05 00:01, , 3F
但是我不知道如何下手...Orz
10/05 00:01, 3F

10/05 03:50, , 4F
妳貼出來的程式 最後只會有1筆 麻煩確認一下
10/05 03:50, 4F
※ 編輯: lovescent 來自: 114.40.95.71 (10/06 22:52)

10/06 22:52, , 5F
少打的部份已更正 謝謝你...:P
10/06 22:52, 5F

10/06 23:34, , 6F
保留j變項 依照j變項進行非線性迴歸 此外 妳要的是參數估計值
10/06 23:34, 6F

10/06 23:34, , 7F
不是用output out= p=pre (這邊是輸出預測值)
10/06 23:34, 7F

10/06 23:35, , 8F
輸出估計值 通常是在proc nlin data= 這邊加outest=filename
10/06 23:35, 8F

10/06 23:35, , 9F
不確定 但你可以試試看
10/06 23:35, 9F
文章代碼(AID): #1AnwDyl5 (Statistics)