异方差 Python实现

最新推荐文章于 2024-08-23 19:51:01 发布

python-ing

最新推荐文章于 2024-08-23 19:51:01 发布

阅读量2.8k

点赞数

分类专栏： python 轴承分析

本文链接：https://blog.youkuaiyun.com/weixin_43174381/article/details/103411555

版权

借鉴一位博主的文章，结合两篇，用于实现Python代码，本人亲测，可以跑通。
这位博主的代码字母有连一起的，以下是我亲测的代码：
http://www.sohu.com/a/203047458_654419
http://www.sohu.com/a/201764621_654419

import pandas as pd

import numpy as np

from patsy import dmatrices

from statsmodels.stats.outliers_influence import variance_inflation_factor

import statsmodels.api as sm

import scipy.stats as stats

from sklearn.metrics import mean_squared_error

import seaborn as sns

import matplotlib.pyplot as plt

import matplotlib.mlab as mlab
import scipy.io
plt.rcParams['font.sans-serif']=['SimHei']
plt.rcParams['axes.unicode_minus'] = False
# 数据读取

# #ccpp = pd.read_excel( 'CCPP.xlsx')ccpp.describe()
# data = scipy.io.loadmat('ENCDATA-2hp.mat') # 读取mat文件
# # path = scio.loadmat('FFT-2hp.mat')['FFT-2hp']
# print(data)
# train1=data['train3hp']
# test1=data['test3hp']
# train_y=data['train_y3hp']
# test_y=data['test_y3hp']
# data1=train1[:,1]
# print("train1",train1.shape)
# print("data",data1.shape)
# # sns.pairplot(data)
# # plt.show()
# #y, X = dmatrices( data1, data = t