99999久久久久久亚洲,欧美人与禽猛交狂配,高清日韩av在线影院,一个人在线高清免费观看,啦啦啦在线视频免费观看www

熱線電話:13121318867

登錄
首頁精彩閱讀SAS、spss進(jìn)行Durbin-Watson檢驗
SAS、spss進(jìn)行Durbin-Watson檢驗
2016-12-22
收藏

SAS、spss進(jìn)行Durbin-Watson檢驗

1.首先要知道你用的哪一個SAS子程序。
2.如果是PROC NLIN,那么非常遺憾,它沒有現(xiàn)成的DW統(tǒng)計量,但是你可以在OUTPUT選項中在輸出數(shù)據(jù)集里面輸出殘差。
3.DW實際上是對殘差做一階自相關(guān)判斷,因此你完全可以根據(jù)公式用SAS中的DATA步來完成代碼開發(fā),參考代碼如下:
data nkwilling;
do i=1 to 100;
e=normal(0);
output;
end;
run;

data DW;
set nkwilling end=last;
e_lag=lag(e);
e_dif=sum(e,-e_lag);
t1=e_dif*e_dif;
t2=e*e;
if _n_=1 then do;dw1=t1;dw2=t2;end;
else do;dw1+t1;dw2+t2;end;
if last then dw=dw1/dw2;
run;
再參考DW有關(guān)自相關(guān)的范圍,我記得好像是0-4,作出判斷。
匆忙寫的,你再參考有關(guān)書籍做一下修改。

sas 和spss都能做
檢驗圖就是殘差圖
以預(yù)測值Y為橫軸,以y與預(yù)測值Y之間的誤差et為縱軸(或?qū)W生化殘差與擬和值或一個自變量),繪制殘差的散點圖。如果散點呈現(xiàn)出明顯的規(guī)律性,則認(rèn)為存在自相關(guān)性或者非線性或者非常數(shù)方差的問題。

DW是0<D<4,統(tǒng)計學(xué)意義如下:
①當(dāng)殘差與自變量互為獨立時,D=2 或 DW 越接近2,判斷無自相關(guān)性把握越大。
②當(dāng)相鄰兩點的殘差為正相關(guān)時,D<2,DW 越接近于0,正自相關(guān)性越強(qiáng)。
③當(dāng)相鄰兩點的殘差為負(fù)相關(guān)時,D>2,DW 越接近于4,負(fù)自相關(guān)性越強(qiáng)。
判斷。根據(jù)樣本容量n 和解釋變量的數(shù)目p 查DW 分布表,得下臨界值L D 和上臨界值U D ,
并依下列準(zhǔn)則判斷擾動項的自相關(guān)情形。
(1)如果0<DW< L D ,則拒絕零假設(shè),擾動項存在一階正自相關(guān)。DW 越接近于0,正自相關(guān)
性越強(qiáng)。
(2)如果L D <DW< U D ,則無法判斷是否有自相關(guān)。
(3)如果U D <DW<4- U D ,則接受零假設(shè),擾動項不存在一階正自相關(guān)。DW 越接近2,判斷
無自相關(guān)性把握越大。
(4)如果4- U D <DW<4- L D ,則無法判斷是否有自相關(guān)。
(5) 如果4- L D <DW<4,則拒絕零假設(shè),擾動項存在一階負(fù)自相關(guān)。DW 越接近于4,負(fù)自
相關(guān)性越強(qiáng)。


檢驗不難,據(jù)不完全統(tǒng)計,PROC REG/AUTOREG/MODEL都有選項輸出統(tǒng)計量和p-值。你要是非線性的,可以用PROC MODEL。其實Durbin-Watson檢驗的統(tǒng)計量也可以利用殘差根據(jù)公式手工算。

但是檢驗圖是個什么概念不才就一點也不懂了,一個模型不就只有一個Durbin-Watson值嗎?


option nocenter;
dm ‘log;clear;output;clear’;

proc import datafile=”c:\example.xls” replace
out=one;
getnames=yes;
data one;  set one;
proc nlin data=one;
parms b1=0.2 b2=-0.2 b3=-0.4 ;
AOld = a1;anew=a2;
hdold=hd1;temp = AOld / ANew;
do anew = (1+a1) to a2 by 1;
HdNew = exp(temp*log(hdold)+(1-temp)*(b1+b2/aold+b3*hdold));                                   AOld = ANew;
hdold=hdnew;
end;
model hd2 = HdNew;
output out=two predicted=hd2hat;

我應(yīng)該如何在上面代碼中添加Durbin-Watson檢驗?zāi)兀?/span>

用proc reg就好了,在option那里加一個“DW


數(shù)據(jù)分析咨詢請掃描二維碼

若不方便掃碼,搜微信號:CDAshujufenxi

數(shù)據(jù)分析師資訊
更多

OK
客服在線
立即咨詢
客服在線
立即咨詢
') } function initGt() { var handler = function (captchaObj) { captchaObj.appendTo('#captcha'); captchaObj.onReady(function () { $("#wait").hide(); }).onSuccess(function(){ $('.getcheckcode').removeClass('dis'); $('.getcheckcode').trigger('click'); }); window.captchaObj = captchaObj; }; $('#captcha').show(); $.ajax({ url: "/login/gtstart?t=" + (new Date()).getTime(), // 加隨機(jī)數(shù)防止緩存 type: "get", dataType: "json", success: function (data) { $('#text').hide(); $('#wait').show(); // 調(diào)用 initGeetest 進(jìn)行初始化 // 參數(shù)1:配置參數(shù) // 參數(shù)2:回調(diào),回調(diào)的第一個參數(shù)驗證碼對象,之后可以使用它調(diào)用相應(yīng)的接口 initGeetest({ // 以下 4 個配置參數(shù)為必須,不能缺少 gt: data.gt, challenge: data.challenge, offline: !data.success, // 表示用戶后臺檢測極驗服務(wù)器是否宕機(jī) new_captcha: data.new_captcha, // 用于宕機(jī)時表示是新驗證碼的宕機(jī) product: "float", // 產(chǎn)品形式,包括:float,popup width: "280px", https: true // 更多配置參數(shù)說明請參見:http://docs.geetest.com/install/client/web-front/ }, handler); } }); } function codeCutdown() { if(_wait == 0){ //倒計時完成 $(".getcheckcode").removeClass('dis').html("重新獲取"); }else{ $(".getcheckcode").addClass('dis').html("重新獲取("+_wait+"s)"); _wait--; setTimeout(function () { codeCutdown(); },1000); } } function inputValidate(ele,telInput) { var oInput = ele; var inputVal = oInput.val(); var oType = ele.attr('data-type'); var oEtag = $('#etag').val(); var oErr = oInput.closest('.form_box').next('.err_txt'); var empTxt = '請輸入'+oInput.attr('placeholder')+'!'; var errTxt = '請輸入正確的'+oInput.attr('placeholder')+'!'; var pattern; if(inputVal==""){ if(!telInput){ errFun(oErr,empTxt); } return false; }else { switch (oType){ case 'login_mobile': pattern = /^1[3456789]\d{9}$/; if(inputVal.length==11) { $.ajax({ url: '/login/checkmobile', type: "post", dataType: "json", data: { mobile: inputVal, etag: oEtag, page_ur: window.location.href, page_referer: document.referrer }, success: function (data) { } }); } break; case 'login_yzm': pattern = /^\d{6}$/; break; } if(oType=='login_mobile'){ } if(!!validateFun(pattern,inputVal)){ errFun(oErr,'') if(telInput){ $('.getcheckcode').removeClass('dis'); } }else { if(!telInput) { errFun(oErr, errTxt); }else { $('.getcheckcode').addClass('dis'); } return false; } } return true; } function errFun(obj,msg) { obj.html(msg); if(msg==''){ $('.login_submit').removeClass('dis'); }else { $('.login_submit').addClass('dis'); } } function validateFun(pat,val) { return pat.test(val); }