- 博客(7)
- 收藏
- 关注
RSS订阅原创 WLS算法
data_new['e']=data_new['y']-model75.fittedvaluessigma1=np.round((np.sum((data_new[data_new['region']==1]['e'])**2))/(len(data_new[data_new['region']==1]['e'])-1),2)sigma2=np.round((np.sum((data_new[data_new['region']==2]['e'])**2))/(len(data_new[data_new
2021-12-07 19:48:50
1370
原创 logist回归基本代码
data=pd.read_csv("C:/Users/Administrator/Desktop/P3.csv")data['constant']=1#表12.2m1 = sm.Logit(data['y'], data[['constant','x1','x2','x3']]).fit()print(m1.summary())==========================================================================Dep. Vari
2021-12-06 19:26:43
824
原创 特征值特征向量
cor=data.iloc[:,2:].corr()lamda=np.linalg.eig(cor)[0]#0特征值，1特征向量lamda=np.round(lamda,3)table=PrettyTable()table.add_column("Variable",["1",'2','3','4','5','6'])table.add_column("特征值",lamda)print(table)
2021-11-29 21:42:56
405
原创 sklearn岭回归
from sklearn.linear_model import Ridgek=np.linspace(0,1,29)#(start, end, num=num_points)使k属于[0,1]y=aaa['import']x=aaa.iloc[:,1:4]θ1=[]θ2=[]θ3=[]SSE=[]for i in range(len(k)): ridge = Ridge(k[i]) ridge =ridge .fit(x, y) θ1.append(ridge....
2021-11-27 20:35:46
965
原创 岭迹图的制作
# 设置matplotlib参数，正确显示中文和'-'符号。plt.rcParams.update( { 'text.usetex': False, 'font.family': 'stixgeneral', 'mathtext.fontset': 'stix', })alpha=np.linspace(0,1,220)intercept_dict = {}#回归系数保存在这里for a in alpha: result = ...
2021-11-23 22:34:49
1761
原创 定义pca
def PCA(dataSet): #计算平均值 n=len(dataSet) meanVal=np.mean(dataSet,axis=0) #axis=0 , Along the direction of the row #中心化及尺度变换 CenteringMat=dataSet-meanVal StdVal=np.std(dataSet,axis=0) StdVal=(StdVal*np.sqrt(n))/np.sqrt(n-1) Sta...
2021-11-20 18:00:43
1211
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人