com.github.jessemull.microflexbigdecimal.stat.CentralMoment Maven / Gradle / Ivy

Go to download
Show more of this group Show more artifacts with this name
Show all versions of microflex-bigdecimal Show documentation
Microplate library for parsing wet lab data.
The newest version!
/**
 * Licensed to the Apache Software Foundation (ASF) under one
 * or more contributor license agreements.  See the NOTICE file
 * distributed with this work for additional information
 * regarding copyright ownership.  The ASF licenses this file
 * to you under the Apache License, Version 2.0 (the
 * "License"); you may not use this file except in compliance
 * with the License.  You may obtain a copy of the License at
 *
 *   http://www.apache.org/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing,
 * software distributed under the License is distributed on an
 * "AS IS" BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY
 * KIND, either express or implied.  See the License for the
 * specific language governing permissions and limitations
 * under the License.
 */

/* -------------------------------- Package -------------------------------- */

package com.github.jessemull.microflexbigdecimal.stat;

/* ----------------------------- Dependencies ------------------------------ */

import java.math.BigDecimal;
import java.math.MathContext;
import java.util.List;

/**
 * This class pth central moment for plate stacks, plates, wells and well sets.
 *
 * 


 * 
 * Statistical operations can be performed on stacks, plates, sets and wells using
 * standard or aggregated functions. Standard functions calculate the desired
 * statistic for each well in the stack, plate or set. Aggregated functions aggregate
 * the values from all the wells in the stack, plate or set and perform the statistical
 * operation on the aggregated values. Both standard and aggregated functions can
 * be performed on a subset of data within the stack, plate, set or well.
 * 
 * 


 * 
 * The methods within the MicroFlex library are meant to be flexible and the
 * descriptive statistic object supports operations using a single stack, plate,
 * set or well as well as collections and arrays of stacks, plates, sets or wells. 
 *  
 * 
 *    
 *    
 *    
 *    
 *    
 *       
 *       
 *       
 *       
 *    
 *    
 *       
 *       
 *       
 *       
 *    
 * Operation Beginning
Index Length of
Subset Input/Output

 *          
 *             
 *                
 *             
 *          Standard  
 *       
 *          
 *             
 *                
 *             
 *          +/-  
 *       
 *          
 *             
 *                
 *             
 *          +/-  
 *       
 *          
 *             
 *                
 *             
 *             
 *                
 *             
 *          Accepts a single well, set, plate or stack as input
Calculates the statistic for each well in a well, set, plate or stack  
 *       

 *          
 *             
 *                
 *             
 *          Aggregated  
 *       
 *          
 *             
 *                
 *             
 *          +/-  
 *       
 *          
 *             
 *                
 *             
 *          +/-  
 *       
 *          
 *             
 *                
 *             
 *              
 *                
 *             
 *          Accepts a single well/set/plate/stack or a collection/array of wells/sets/plates/stacks as input
Aggregates the data from all the wells in a well/set/plate/stack and calculates the statistic using the aggregated data  
 *       
 * 
 * @author Jesse L. Mull
 * @update Updated Oct 18, 2016
 * @address http://www.jessemull.com
 * @email [email protected]
 */
public class CentralMoment extends QuantileStatisticContext {
    
    /**
     * Calculates the median.
     * @param    List    the list
     * @param    int                 pth central moment
     * @param    MathContext         the math context
     * @return                       the result
     */
    public BigDecimal calculate(List list, int p, MathContext mc) {

    	BigDecimal length = new BigDecimal(list.size());
    	BigDecimal mean = new BigDecimal(0);
    	
    	for(BigDecimal bd : list) {
    		mean = mean.add(bd, mc);
    	}
    	
    	mean = mean.divide(length, mc);
    	
    	BigDecimal result = new BigDecimal(0);
    	
    	for(BigDecimal bd : list) {
    		BigDecimal difference = bd.subtract(mean);
    		result = result.add(difference.pow(p, mc));
    	}
    	
    	return result.divide(length, mc);
        
    }
    
    /**
     * Calculates the median of the values between the beginning and 
     * ending indices.
     * @param    List    the list
     * @param    int             beginning index of subset
     * @param    int             length of subset
     * @param    int             pth central moment
     * @param    MathContext     the math context
     * @return                   the result
     */
    public BigDecimal  calculate(List list, int begin, int length, int p, MathContext mc) {
        return calculate(list.subList(begin, begin + length), p, mc);
    }

}