1、問(wèn)題描述

一套 1000 平方英尺 (sqft) 的房屋售價(jià)為300,000美元，一套 2000 平方英尺的房屋售價(jià)為500,000美元。這兩點(diǎn)將構(gòu)成我們的數(shù)據(jù)或訓(xùn)練集。面積單位為 1000 平方英尺，價(jià)格單位為 1000 美元。

希望通過(guò)這兩個(gè)點(diǎn)擬合線性回歸模型，以便可以預(yù)測(cè)其他房屋的價(jià)格。例如，面積為 1200 平方英尺的房屋價(jià)格是多少。

首先導(dǎo)入所需要的庫(kù)

import numpy as np
import matplotlib.pyplot as plt
plt.style.use('./deeplearning.mplstyle')

以下代碼來(lái)創(chuàng)建x_train和y_train變量。數(shù)據(jù)存儲(chǔ)在一維 NumPy 數(shù)組中。

# x_train is the input variable (size in 1000 square feet)
# y_train is the target (price in 1000s of dollars)
x_train = np.array([1.0, 2.0])
y_train = np.array([300.0, 500.0])
print(f"x_train = {x_train}")
print(f"y_train = {y_train}")

2、表示說(shuō)明

使用 m 來(lái)表示訓(xùn)練樣本的數(shù)量。 (x${}^{(i)}$, y${}^{(i)}$) 表示第 i 個(gè)訓(xùn)練樣本。由于 Python 是零索引的，(x${}^{(0)}$, y${}^{(0)}$) 是 (1.0, 300.0) ， (x${}^{(1)}$, y${}^{(1)}$) 是 (2.0, 500.0).

3、數(shù)據(jù)繪圖

使用 matplotlib 庫(kù)中的scatter()函數(shù)繪制這兩個(gè)點(diǎn)。 其中，函數(shù)參數(shù)marker 和 c 將點(diǎn)顯示為紅叉（默認(rèn)為藍(lán)點(diǎn)）。使用matplotlib庫(kù)中的其他函數(shù)來(lái)設(shè)置要顯示的標(biāo)題和標(biāo)簽。

# Plot the data points
plt.scatter(x_train, y_train, marker='x', c='r')
# Set the title
plt.title("Housing Prices")
# Set the y-axis label
plt.ylabel('Price (in 1000s of dollars)')
# Set the x-axis label
plt.xlabel('Size (1000 sqft)')
plt.show()

4、模型函數(shù)

線性回歸的模型函數(shù)（這是一個(gè)從 x 映射到 y 的函數(shù)）可以表示為 $$f_{w,b}(x^{(i)})=w{x}^{(i)}+b\text{\hspace{1em}(1)}$$

計(jì)算$f_{w,b}(x^{(i)})$ 的值，可以將每個(gè)數(shù)據(jù)點(diǎn)顯示地寫(xiě)為：

對(duì)于 $x^{(0)}$, f_wb = w * x[0] + b
對(duì)于 $x^{(1)}$, f_wb = w * x[1] + b

對(duì)于大量的數(shù)據(jù)點(diǎn)，這可能會(huì)變得笨拙且重復(fù)。 因此，可以在for 循環(huán)中計(jì)算輸出，如下面的函數(shù)compute_model_output 所示。

def compute_model_output(x, w, b):
    """
    Computes the prediction of a linear model
    Args:
      x (ndarray (m,)): Data, m examples 
      w,b (scalar)    : model parameters  
    Returns
      y (ndarray (m,)): target values
    """
    m = x.shape[0]
    f_wb = np.zeros(m)
    for i in range(m):
        f_wb[i] = w * x[i] + b
        
    return f_wb

調(diào)用 compute_model_output 函數(shù)并繪制輸出

w = 100
b = 100

tmp_f_wb = compute_model_output(x_train, w, b,)

# Plot our model prediction
plt.plot(x_train, tmp_f_wb, c='b',label='Our Prediction')

# Plot the data points
plt.scatter(x_train, y_train, marker='x', c='r',label='Actual Values')

# Set the title
plt.title("Housing Prices")
# Set the y-axis label
plt.ylabel('Price (in 1000s of dollars)')
# Set the x-axis label
plt.xlabel('Size (1000 sqft)')
plt.legend()
plt.show()

很明顯， $w=100$ 和 $b=100$ 不會(huì)產(chǎn)生適合數(shù)據(jù)的直線。

根據(jù)學(xué)過(guò)的數(shù)學(xué)知識(shí)，可以容易求出 $w=200$ 和 $b=100$

5、預(yù)測(cè)

現(xiàn)在我們已經(jīng)有了一個(gè)模型，可以用它來(lái)做出房屋價(jià)格的預(yù)測(cè)。來(lái)預(yù)測(cè)一下 1200 平方英尺的房子的價(jià)格。由于面積單位為 1000 平方英尺，所以 $x$ 是1.2。

w = 200                         
b = 100    
x_i = 1.2
cost_1200sqft = w * x_i + b    

print(f"${cost_1200sqft:.0f} thousand dollars")

輸出的結(jié)果是：$340 thousand dollars

總結(jié)

• 線性回歸建立一個(gè)特征和目標(biāo)之間關(guān)系的模型
  • 在上面的例子中，特征是房屋面積，目標(biāo)是房?jī)r(jià)。
  • 對(duì)于簡(jiǎn)單線性回歸，模型有兩個(gè)參數(shù) $w$ 和 $b$ ，其值使用訓(xùn)練數(shù)據(jù)進(jìn)行擬合。
  • 一旦確定了模型的參數(shù)，該模型就可以用于對(duì)新數(shù)據(jù)進(jìn)行預(yù)測(cè)。

附錄

deeplearning.mplstyle 源碼：文章來(lái)源地址http://www.zghlxwxcb.cn/news/detail-624117.html

# see https://matplotlib.org/stable/tutorials/introductory/customizing.html
lines.linewidth: 4
lines.solid_capstyle: butt

legend.fancybox: true

# Verdana" for non-math text,
# Cambria Math

#Blue (Crayon-Aqua) 0096FF
#Dark Red C00000
#Orange (Apple Orange) FF9300
#Black 000000
#Magenta FF40FF
#Purple 7030A0

axes.prop_cycle: cycler('color', ['0096FF', 'FF9300', 'FF40FF', '7030A0', 'C00000'])
#axes.facecolor: f0f0f0 # grey
axes.facecolor: ffffff  # white
axes.labelsize: large
axes.axisbelow: true
axes.grid: False
axes.edgecolor: f0f0f0
axes.linewidth: 3.0
axes.titlesize: x-large

patch.edgecolor: f0f0f0
patch.linewidth: 0.5

svg.fonttype: path

grid.linestyle: -
grid.linewidth: 1.0
grid.color: cbcbcb

xtick.major.size: 0
xtick.minor.size: 0
ytick.major.size: 0
ytick.minor.size: 0

savefig.edgecolor: f0f0f0
savefig.facecolor: f0f0f0

#figure.subplot.left: 0.08
#figure.subplot.right: 0.95
#figure.subplot.bottom: 0.07

#figure.facecolor: f0f0f0  # grey
figure.facecolor: ffffff  # white

## ***************************************************************************
## * FONT                                                                    *
## ***************************************************************************
## The font properties used by `text.Text`.
## See https://matplotlib.org/api/font_manager_api.html for more information
## on font properties.  The 6 font properties used for font matching are
## given below with their default values.
##
## The font.family property can take either a concrete font name (not supported
## when rendering text with usetex), or one of the following five generic
## values:
##     - 'serif' (e.g., Times),
##     - 'sans-serif' (e.g., Helvetica),
##     - 'cursive' (e.g., Zapf-Chancery),
##     - 'fantasy' (e.g., Western), and
##     - 'monospace' (e.g., Courier).
## Each of these values has a corresponding default list of font names
## (font.serif, etc.); the first available font in the list is used.  Note that
## for font.serif, font.sans-serif, and font.monospace, the first element of
## the list (a DejaVu font) will always be used because DejaVu is shipped with
## Matplotlib and is thus guaranteed to be available; the other entries are
## left as examples of other possible values.
##
## The font.style property has three values: normal (or roman), italic
## or oblique.  The oblique style will be used for italic, if it is not
## present.
##
## The font.variant property has two values: normal or small-caps.  For
## TrueType fonts, which are scalable fonts, small-caps is equivalent
## to using a font size of 'smaller', or about 83%% of the current font
## size.
##
## The font.weight property has effectively 13 values: normal, bold,
## bolder, lighter, 100, 200, 300, ..., 900.  Normal is the same as
## 400, and bold is 700.  bolder and lighter are relative values with
## respect to the current weight.
##
## The font.stretch property has 11 values: ultra-condensed,
## extra-condensed, condensed, semi-condensed, normal, semi-expanded,
## expanded, extra-expanded, ultra-expanded, wider, and narrower.  This
## property is not currently implemented.
##
## The font.size property is the default font size for text, given in points.
## 10 pt is the standard value.
##
## Note that font.size controls default text sizes.  To configure
## special text sizes tick labels, axes, labels, title, etc., see the rc
## settings for axes and ticks.  Special text sizes can be defined
## relative to font.size, using the following values: xx-small, x-small,
## small, medium, large, x-large, xx-large, larger, or smaller


font.family:  sans-serif
font.style:   normal
font.variant: normal
font.weight:  normal
font.stretch: normal
font.size:    8.0

font.serif:      DejaVu Serif, Bitstream Vera Serif, Computer Modern Roman, New Century Schoolbook, Century Schoolbook L, Utopia, ITC Bookman, Bookman, Nimbus Roman No9 L, Times New Roman, Times, Palatino, Charter, serif
font.sans-serif: Verdana, DejaVu Sans, Bitstream Vera Sans, Computer Modern Sans Serif, Lucida Grande, Geneva, Lucid, Arial, Helvetica, Avant Garde, sans-serif
font.cursive:    Apple Chancery, Textile, Zapf Chancery, Sand, Script MT, Felipa, Comic Neue, Comic Sans MS, cursive
font.fantasy:    Chicago, Charcoal, Impact, Western, Humor Sans, xkcd, fantasy
font.monospace:  DejaVu Sans Mono, Bitstream Vera Sans Mono, Computer Modern Typewriter, Andale Mono, Nimbus Mono L, Courier New, Courier, Fixed, Terminal, monospace


## ***************************************************************************
## * TEXT                                                                    *
## ***************************************************************************
## The text properties used by `text.Text`.
## See https://matplotlib.org/api/artist_api.html#module-matplotlib.text
## for more information on text properties
#text.color: black

到了這里，關(guān)于【機(jī)器學(xué)習(xí)】Linear Regression的文章就介紹完了。如果您還想了解更多內(nèi)容，請(qǐng)?jiān)谟疑辖撬阉鱐OY模板網(wǎng)以前的文章或繼續(xù)瀏覽下面的相關(guān)文章，希望大家以后多多支持TOY模板網(wǎng)！