stat_txt2.py

import numpy as np
import pandas as pd
import sys

def calculate_statistics(file_path):
    # 读取txt文件并提取第三列数据
    data = pd.read_csv(file_path, delimiter='\t', header=None)
    column_3 = data.iloc[:, 2].values

    # 计算均值、标准差、众数和中位数
    mean = np.mean(column_3)
    std = np.std(column_3)

    median = np.median(column_3)

    # 计算大于0.5的百分比
    percent_0_354 = (column_3 > 0.354).mean() * 100

    # 计算0.125-0.5的百分比
    percent_0_177_to_0_354 = ((column_3 >= 0.177) & (column_3 <= 0.354)).mean() * 100
    percent_0_0884_to_0_177=((column_3 >= 0.0884) & (column_3 <= 0.177)).mean()*100
    percent_0_0442_to_0_0884=((column_3 >= 0.0442) & (column_3 <= 0.0884)).mean()*100
    # 计算小于0.5的百分比
    percent_lt_0_5 = (column_3 < 0.0442).mean() * 100

    results = {
        "均值": mean,
        "标准差": std,

        "中位数": median,
        "大于0.354的百分比": percent_0_354,
        "0.177-0.354的百分比": percent_0_177_to_0_354,
        "0.0884-0.177":percent_0_0884_to_0_177,
        "0.0442-0.0884":percent_0_0442_to_0_0884,
        "小于0.0442的百分比": percent_lt_0_5
    }

    return results

if __name__ == '__main__':
    file_path = sys.argv[1]  # 从命令行参数获取文件路径
    results = calculate_statistics(file_path)

    for key, value in results.items():
        print(key, ":", value)