1. 程式人生 > >基於正太分佈的一元離群點檢測

基於正太分佈的一元離群點檢測

以下為Matlab實現程式:

clc;
clear;
%利用最大似然估計,檢測一元離群點
%氣溫測試資料,單位為攝氏度,這裡假設資料服從正太分佈
TData=[24,24.1,28.9,28.9,29.0,29.1,29.1,29.2,29.2,29.3,29.4,29,28,29.2,29.3,29,29.5,29.3,29,29,29,29,29,29.4,29.1];
len=length(TData);
%利用最大似然估計計算均值和方差
phat=mle(TData);
mu=phat(1);
sigma=phat(2);
%定義離群點集
Outliers=[];
%落在 mu+/-3sigma區間範圍之外的資料即為離群點
for i=1:len
    a=abs(TData(i)-mu)/sigma;
    if a>3
        Outliers(length(Outliers)+1)=TData(i);
    end
end