用C语言实现两个样本t-test

#include <stdio.h>
#include <math.h>

double calculateMean(double data[], int size) {
    double sum = 0;
    for (int i = 0; i < size; i++) {
        sum += data[i];
    }
    return sum / size;
}

double calculateVariance(double data[], int size, double mean) {
    double sum = 0;
    for (int i = 0; i < size; i++) {
        sum += pow(data[i] - mean, 2);
    }
    return sum / (size - 1);
}

double calculateTValue(double data1[], int size1, double data2[], int size2) {
    double mean1 = calculateMean(data1, size1);
    double mean2 = calculateMean(data2, size2);
    double variance1 = calculateVariance(data1, size1, mean1);
    double variance2 = calculateVariance(data2, size2, mean2);
    double pooledVariance = ((size1 - 1) * variance1 + (size2 - 1) * variance2) / (size1 + size2 - 2);
    double tValue = (mean1 - mean2) / sqrt(pooledVariance * (1.0 / size1 + 1.0 / size2));
    return tValue;
}

int main() {
    double data1[] = {1.2, 2.5, 3.1, 4.2, 5.7};
    double data2[] = {0.8, 2.1, 3.5, 4.9, 6.3};
    int size1 = sizeof(data1) / sizeof(data1[0]);
    int size2 = sizeof(data2) / sizeof(data2[0]);

    double tValue = calculateTValue(data1, size1, data2, size2);
    printf("t-value: %f\n", tValue);

    return 0;
}