相关系数t检验 java_java实现t检验

描述

用于两个总体均值之间的比较。

Demo

public class TtestDemo {

/**

* two-sided test

* @param args

*/

public static void main(String[] args) {

double[] x= {300,280,344,385,372,360,288,321,376,290,301,283};

double[] y= {274,220,308,336,198,300,315,258,318,310,332,263};

Ttest t=new Ttest(x, y);

System.out.println(t.getPValue());

}

}

实现代码

package com.math.statistics;

import org.apache.commons.math3.stat.descriptive.moment.StandardDeviation;

import JSci.maths.statistics.TDistribution;

/***

* @author miaoyibo

*

*/

public class Ttest {

private double[] x;

private double[] y;

StandardDeviation standardDeviation =new StandardDeviation();

public Ttest(double[] x, double[] y) {

super();

this.x = x;

this.y = y;

}

public int getXSize() {

return x==null?0:x.length;

}

public int getYSize() {

return y==null?0:y.length;

}

public double getXMean() {

int n=x.length;

double sum=0;

for (double d : x) {

sum=sum+d;

}

return sum/n;

}

public double getYMean() {

int n=y.length;

double sum=0;

for (double d : y) {

sum=sum+d;

}

return sum/n;

}

public double getStandard(double[] x) {

return standardDeviation.evaluate(x);

}

public double calculateTvalue() {

double a=getXMean()-getYMean();

double q1=getStandard(x);

double q2=getStandard(y);

double s1=q1*q1;

double s2=q2*q2;

double b=Math.sqrt(s1/x.length+s2/y.length);

return a/b;

}

public int getDegreesOfFreedom() {

double q1=getStandard(x);

double q2=getStandard(y);

double s1=q1*q1;

double s2=q2*q2;

double a=(s1/x.length+s2/y.length)*(s1/x.length+s2/y.length);

double b=((s1/x.length)*(s1/x.length))/(x.length-1)+((s2/y.length)*(s2/y.length))/(y.length-1);

Double free=a/b;

return free.intValue();

}

public double getPValue() {

int free=getDegreesOfFreedom();

double t=calculateTvalue();

TDistribution td=new TDistribution(free);

double cumulative = td.cumulative(t);

double p;

if(t>0) {

p=(1-cumulative)*2;

}else {

p=cumulative*2;

}

return p;

}

}