绘制残差图

def solve():
    import numpy as np
    from sklearn.datasets import load_iris
    iris = load_iris()
    X = iris.data
    y = iris.target
    # 重复测量：使用4个特征列作为4个时间点
    from scipy import stats
    stat, p = stats.friedmanchisquare(____)
    print(____)