Package fork.lib.math.algebra.elementary.function.v1.distr

Examples of fork.lib.math.algebra.elementary.function.v1.distr.ParameterEstimation

315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
* @param sigma   standard deviation
* @return  
*/
public FunctionV1 kernelDensityEstimation(double sigma){
    try{
        FunctionV1 fn= new ConstantFunction(0);
        Iterator<Double> it= cs.keySet().iterator();
        while(it.hasNext()){
            Double v= it.next();
            Integer f= cs.get(v);
            double coeff= (double) f.intValue()/size();
            fn= fn.add( (new NormalDistribution( v.doubleValue(), sigma)).
                    multiply(new ConstantFunction(coeff)) );
        }
        return fn;
    }catch(FunctionException e){
        e.printStackTrace();
View Full Code Here
119
120
121
122
123
124
125
126
127
128
129
* @param coeffs
* @return
* @throws FunctionException
*/
public static double logLikelihood(Distribution data, ArrayList<DistributionFunction> distrs, ArrayList<Double> coeffs) throws FunctionException{
    FunctionV1 fn= new ConstantFunction(0);
    for( int i=0; i<distrs.size(); i++ ){
        fn = fn.add( distrs.get(i).multiply(new ConstantFunction(coeffs.get(i))) );
    }
    return data.logLikelihood(fn);
}
View Full Code Here
93
94
95
96
97
98
99
100
101
102
103
104
    double mm= Double.POSITIVE_INFINITY;
    double ms= Double.POSITIVE_INFINITY;
    if(res!=null){
        ArrayList<DistributionFunction> distrs= res.distrs;
        for( int i=0; i<distrs.size(); i++ ){
            DistributionFunction df= distrs.get(i);
            NormalDistributionParam param= (NormalDistributionParam)df.getParameters();
            if(param.mu<mm){
                mm= param.mu;
                ms= param.sigma;
            }
        }
View Full Code Here
321
322
323
324
325
326
327
328
329
330
331
        Iterator<Double> it= cs.keySet().iterator();
        while(it.hasNext()){
            Double v= it.next();
            Integer f= cs.get(v);
            double coeff= (double) f.intValue()/size();
            fn= fn.add( (new NormalDistribution( v.doubleValue(), sigma)).
                    multiply(new ConstantFunction(coeff)) );
        }
        return fn;
    }catch(FunctionException e){
        e.printStackTrace();
View Full Code Here
450
451
452
453
454
455
456
457
458
459
460



public static void main(String[] args) throws Exception { //debug
    NormalDistribution n= new NormalDistribution(10,5);
    Distribution d= new Distribution(
            0,0,0,0,1,1,1,1,1.2,1.2,1.3,2,2,3,5,6,7,
            10,10,10,10,11,11,11,11,12,12,12,14,13,13
            );
    d.add(20, 5); d.add(21, 10); d.add(22, 8);
View Full Code Here
83
84
85
86
87
88
89
90
91
92
93
94
* @return
*/
protected double threshold() throws EMException, DistributionException, FunctionException{
    Distribution d= new Distribution(data.getValues());
    DistributionFunction[] fn= new DistributionFunction[]{
        new NormalDistribution(0,1),
        new NormalDistribution(0,1)
    };
    EM em= new EM(d, fn, null);
    EMResultEntry res= em.getOptimalResultEntry();
    double mm= Double.POSITIVE_INFINITY;
    double ms= Double.POSITIVE_INFINITY;
View Full Code Here
84
85
86
87
88
89
90
91
92
93
94
protected ArrayList<DistributionFunction> findStartDistrs() throws EMException, DistributionException{
    ArrayList<DistributionFunction> ret= new ArrayList<>();
    ArrayList<Distribution> subs= data.quantileSubsets(types.length);
    for(int i=0;i<types.length;i++){
        if(types[i] instanceof NormalDistribution){
            ParameterEstimation pe = new ParameterEstimation(subs.get(i), new NormalDistribution(0, 1));
            ret.add(pe.funcOut);
        } else if(types[i] instanceof NormalDistribution){
            throw new EMException();
        } else {
            throw new EMException();
View Full Code Here
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
    d.add(7, 5); d.add(8, 8); d.add(9, 10); d.add(10,5);
    d.add(20, 5); d.add(21, 10); d.add(22, 8);
    d.add(120, 5); d.add(121, 10); d.add(122, 8);
   
    DistributionFunction[] types= new DistributionFunction[]{
        new NormalDistribution(0,1),
        new NormalDistribution(0,1),
        new NormalDistribution(0,1),
        //new NormalDistribution(0,1),
        new NormalDistribution(0,1)
    };
   
   
    EM em= new EM(d, types, null);
    em.start();
View Full Code Here
94
95
96
97
98
99
100
101
102
103
104
    double ms= Double.POSITIVE_INFINITY;
    if(res!=null){
        ArrayList<DistributionFunction> distrs= res.distrs;
        for( int i=0; i<distrs.size(); i++ ){
            DistributionFunction df= distrs.get(i);
            NormalDistributionParam param= (NormalDistributionParam)df.getParameters();
            if(param.mu<mm){
                mm= param.mu;
                ms= param.sigma;
            }
        }
View Full Code Here
84
85
86
87
88
89
90
91
92
93
94
protected ArrayList<DistributionFunction> findStartDistrs() throws EMException, DistributionException{
    ArrayList<DistributionFunction> ret= new ArrayList<>();
    ArrayList<Distribution> subs= data.quantileSubsets(types.length);
    for(int i=0;i<types.length;i++){
        if(types[i] instanceof NormalDistribution){
            ParameterEstimation pe = new ParameterEstimation(subs.get(i), new NormalDistribution(0, 1));
            ret.add(pe.funcOut);
        } else if(types[i] instanceof NormalDistribution){
            throw new EMException();
        } else {
            throw new EMException();
View Full Code Here
TOP

Related Classes of fork.lib.math.algebra.elementary.function.v1.distr.ParameterEstimation

  • fork.lib.math.algebra.elementary.set.continuous.Region
  • fork.lib.math.applied.dm.peak.EMPeakDetector
  • fork.lib.math.applied.optim.EM
  • fork.lib.math.applied.optim.EMCycle$EMParameterEstimation
  • fork.lib.math.applied.stat.Distribution
  • fork.lib.math.algebra.elementary.function.FunctionException
    • Copyright © 2018 www.massapicom. All rights reserved.
    All source code are property of their respective owners. Java is a trademark of Sun Microsystems, Inc and owned by ORACLE Inc. Contact coftware#gmail.com.