当前位置: 首页 > news >正文

上海市做网站今日头条最新版

上海市做网站,今日头条最新版,logo制作免费版,wordpress函数语言今日是代码对sigmoid函数的实现和运用 #linear_model线性回归 #名字虽然叫逻辑回归,作用于分类 #分类:类别 #回归:预测 from sklearn.linear_model import LogisticRegression 实现函数 import numpy as np import matplotlib.pyplot as pl…

今日是代码对sigmoid函数的实现和运用

#linear_model线性回归
#名字虽然叫逻辑回归,作用于分类
#分类:类别
#回归:预测
from sklearn.linear_model import LogisticRegression

实现函数

import numpy as np
import matplotlib.pyplot as pltdef sigmoid(x):return 1/(1+np.exp(-x))x = np.linspace(-5,5,100)y = sigmoid(x)plt.plot(x,y,color='green')

损失函数

from sklearn import datasets
from sklearn.linear_model import LogisticRegression
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
#Z-score归一化
from sklearn.preprocessing import scale,StandardScaler

加载数据

X,y=datasets.load_breast_cancer(return_X_y=True)X=X[:,:2]#切片两个特征
display(X.shape)
display(y.shape)

建模

model = LogisticRegression()#训练,和之前线性回归,类似
#后面其他方法,算法,类似
model.fit(X,y)

逻辑回归线性方程拿出来:系数
w1 = model.coef_[0,0]
w2 = model.coef_[0,1]
b = model.intercept_
print('方程系数',w1,w2)
print('截距',b)

sigmoid函数

def sigmoid(X,w1,w2,b):z = w1*X[0] + w2*X[1] + b#方程表示return 1/(1+ np.exp(-z))

损失函数

def loss_function(X,y,w1,w2,b):loss = 0for X_i,y_i in zip(X,y):p = sigmoid(X_i,w1,w2,b)#概率p=np.clip(p,0.0001,0.999)#裁剪loss+= -y_i * np.log(p) +(1-y_i)* np.log(1-p)return loss

定义参数w1,w2取值空间

w1_space= np.linspace (w1 - 2,w1 +2,100)w2_space = np.linspace(w2 - 2,w2 +2,100)

损失计算

loss1_ = np.array([loss_function(X,y,i,w2,b) for i in w1_space])
loss1_

loss2_ = np.array([loss_function(X,y,w1,i,b) for i in w2_space])
loss2_

可视化¶

fig1 = plt.figure(figsize=(12,9))plt.subplot(2,2,1)
plt.plot(w1_space,loss1_,color='green')plt.subplot(2,2,2)
plt.plot(w1_space,loss1_,color='red')

逻辑回归代码实现

import numpy as np
from sklearn import datasets
from sklearn.linear_model import LogisticRegression
#将数据拆分
from sklearn.model_selection import train_test_split
X,y = datasets.load_iris(return_X_y=True)
cood = y!=2#过滤数据:类别是2,过滤掉
X=X[cood]
y=y[cood]
y

加载数据并拆分
#将调练数据测试数据:80% 训练数据,保留20%,测试数据
X_train,X_test,y_train,y_test = train_test_split(X,y,test_size=0.2)display(X_train.shape,X_test.shape)
display(y_train.shape,y_test.shape)

训练
model = LogisticRegression()model.fit(X_train,y_train)y_pred = model.predict(X_test)
print('预测结果是:',y_pred)proba_ = model.predict_proba(X_test)
print('预测概率是:\n',proba_)

y_pred

proba_.argmax(axis=1)

概率手动计算
def sigmoid(x):return 1/(1+np.exp(-z))#方程系数和截距
w=model.coef_
b=model.intercept_#求解线性方程
z=X_test.dot(w.reshape(-1))+bp=sigmoid(z)#列合并
#np.column_stack([1-p,p])
np.concatenate([(1-p).reshape(-1,1),p.reshape(-1,1)],axis = 1)[:5]

model.predict_proba(X_test)

http://www.yidumall.com/news/58534.html

相关文章:

  • 网站模板设计开发奇零seo赚钱培训
  • 企业集团网站建设方案网站优化服务
  • 小企业网站服务器seo技术培训
  • ...温岭做网站品牌推广策划
  • 怎么找网站优化官网咨询
  • 构建动态网站嘉兴百度快照优化排名
  • 用老域名做新网站蜘蛛搜索引擎
  • 外贸怎么做网站深圳网站页面设计
  • 手机黄山网站成都网站建设创新互联
  • 东莞网站建设做网站sem优化是什么意思
  • 海淀网站制作营销活动策划
  • 网站 广州谷歌google下载
  • wordpress 在线聊天插件福建seo排名培训
  • 龙华营销型网站制作哪家好查询网站注册信息
  • 网站制作公司 全贵州社会新闻最新消息
  • 网站加载很慢营业推广策略
  • 易旅游网站建设企业网搭建
  • wordpress教程下载网站主题搜索引擎优化排名品牌
  • 北海网站开发网络营销策略案例
  • 高端网站建设 司法搜索优化整站优化
  • 做网站的客户需求百度帐号个人中心
  • 网站的数据库丢失百度搜索关键词
  • 单位网站建设的优势免费行情软件app网站下载大全
  • 一个网站锚文本可以做几个市场调研报告ppt模板
  • 做网站的网络公司好看的web网页
  • 网站建设字体颜色代码seo教程seo入门讲解
  • 没网站做推广网站快速收录软件
  • 政府网站集约化建设模式网站友链查询源码
  • 淮南北京网站建设itmc平台seo优化关键词个数
  • 广州一起做的网站网络营销的成功案例有哪些