em_invgausscdf.m

EM分群,matlab程式碼,用來分群用的

function y = EM_invgausscdf(x,par)
%INVGAUSSPDF Inverse Gaussian cumulative distribution function.
%   Y = INVGAUSSPDF(X,MU,LAM) returns the cumulative of the inverse gaussian function 
%   with parameters MU and LAM, at the values in X.
%
%   The size of Y is the size of the X argument. (GENERALIZE)
%

%   References (CHECK OUT):
%      [1]  .....
% 
%   Created by Enrico Rossoni - 17/06/2004
mu = par(1);
lam = par(2);


% Initialize Y to zero.
y = zeros(size(x));

x1 = sqrt(lam./x);
x2 = x./mu;

y = normcdf(x1.*(x2-1), 0,1) + exp(2*lam./mu).* normcdf(-x1.*(x2+1), 0,1);