em_invert_pars.m

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

function [p] = EM_Invert_pars(mu, sig, fun_type)


switch fun_type
    
    case 'InvGauss'  
              
        p = [mu mu^3/sig^2];   
    
    case 'Gam'  
                
        p = [mu^2/sig^2 sig^2/mu];
        
    case 'LogNorm' 
        
        p = [log(mu)-0.5*log(1+(sig/mu)^2) sqrt(log(1+(sig/mu)^2))];        
        
end