作者lovescent ( )
站內Statistics
標題[程式] sas跑非線性迴歸
時間Sun Oct 4 03:18:17 2009
------------------------------------------------------------------------
[軟體程式類別]: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
1F：→ tew:妳的a最後有多少觀察值，妳有檢驗嗎? 10/04 12:04
2F：→ lovescent:有一萬 可是我是希望每一百筆資料可以跑一次回歸 10/05 00:01
3F：→ lovescent:但是我不知道如何下手...Orz 10/05 00:01
4F：推 tew:妳貼出來的程式 最後只會有1筆 麻煩確認一下 10/05 03:50
※ 編輯: lovescent 來自: 114.40.95.71 (10/06 22:52)
5F：→ lovescent:少打的部份已更正 謝謝你...:P 10/06 22:52
6F：→ tew:保留j變項 依照j變項進行非線性迴歸 此外 妳要的是參數估計值 10/06 23:34
7F：→ tew:不是用output out= p=pre (這邊是輸出預測值) 10/06 23:34
8F：→ tew:輸出估計值 通常是在proc nlin data= 這邊加outest=filename 10/06 23:35
9F：→ tew:不確定 但你可以試試看 10/06 23:35