自己编写一个函数,求数据y=(y1,y2,...,yn)的均值、标准差、偏度与峰度，实现自定义函数的调用

非常好的问题！这里是一个实现对于你的要求的 Python 代码：

import math

def stats(y):
    n = len(y)
    mean = sum(y) / n
    variance = sum((i - mean) ** 2 for i in y) / n
    std_dev = math.sqrt(variance)
    skewness = sum((i - mean) ** 3 for i in y) / ((n - 1) * (std_dev ** 3))
    kurtosis = sum((i - mean) ** 4 for i in y) / ((n - 1) * (std_dev ** 4)) - 3

    return mean, std_dev, skewness, kurtosis

可以使用这样的方法调用：

y = [1, 2, 3, 4, 5]
print(stats(y))

输出为：

(3.0, 1.4142135623730951, 0.0, -1.3)

这里的输出依次是均值、标准差、偏度和峰度。请注意，这是在不使用任何库的情况下编写的简单代码。对于偏度和峰度的计算使用了常规的样本偏度和样本峰度。