rdiagnoise_d.m

来自「计量工具箱」· M 代码 · 共 29 行

% PURPOSE: demo of rdiagnose()
%          Regression Diagnostic routines
%          
%---------------------------------------------------
% USAGE: rdiagnose_d
%---------------------------------------------------


% generate data set
n = 100;
k = 3;
x = randn(n,k);
x(:,1) = ones(n,1);

beta = ones(k,1);

y = x*beta + randn(n,1);

% add outliers
y(50,1) = 10;
y(70,1) = -10;


result = rdiagnose(y,x);

prt = 1;
result = rdiagnose(y,x,prt);