umcg.genetica.methylation.ConvertBetaAndMvalues Maven / Gradle / Ivy

Go to download

Show more of this group Show more artifacts with this name
Show all versions of genetica-libraries Show documentation

There is a newer version: 1.0.7

/*
 * To change this template, choose Tools | Templates
 * and open the template in the editor.
 */
package umcg.genetica.methylation;

import cern.colt.matrix.tdouble.DoubleMatrix1D;
import cern.colt.matrix.tdouble.DoubleMatrix2D;
import cern.jet.math.tdouble.DoubleFunctions;

/**
 *
 * @author MarcJan
 */
public class ConvertBetaAndMvalues {
    
    public static void transformToMvalue(double[][] rawData){
        double minValue = Double.MAX_VALUE;
        double maxValue = Double.MIN_VALUE;
        
        int probeCount = rawData.length;
        int sampleCount = rawData[probeCount-1].length;
        
        for (int p=0; pmaxValue) {
                    maxValue = rawData[p][s];
                }
            }
        }
        
        double multiplier = 1.0d / Math.log10(2.0d);
        minValue = (double) (Math.log10((minValue/(1-minValue))) * multiplier);
        maxValue = (double) (Math.log10((maxValue/(1-maxValue))) * multiplier);
        
        for (int p=0; p