Hướng dẫn code for normal distribution in python - mã cho phân phối chuẩn trong python

Ngay cả khi bạn không ở trong lĩnh vực thống kê, bạn cũng phải bắt gặp thuật ngữ phân phối bình thường.Normal Distribution”.

Phân phối xác suất là một hàm thống kê mô tả khả năng thu được các giá trị có thể mà một biến ngẫu nhiên có thể thực hiện. Bằng cách này, chúng tôi có nghĩa là phạm vi của các giá trị mà một tham số có thể thực hiện khi chúng tôi chọn ngẫu nhiên các giá trị từ nó.

Một phân phối xác suất có thể rời rạc hoặc liên tục.

Giả sử ở một thành phố, chúng ta có chiều cao của người lớn trong độ tuổi từ 20-30 tuổi, từ 4,5 ft đến 7 ft.

Nếu chúng tôi được yêu cầu nhận 1 người lớn ngẫu nhiên và hỏi những gì anh ấy/cô ấy (giả sử giới tính không ảnh hưởng đến chiều cao) sẽ là gì? Không có cách nào để biết chiều cao sẽ là gì. Nhưng nếu chúng ta có sự phân phối chiều cao của người lớn trong thành phố, chúng ta có thể đặt cược vào kết quả có thể xảy ra nhất.

Phân phối bình thường là gì?

Một phân phối bình thường còn được gọi là phân phối Gaussian hoặc đường cong chuông nổi tiếng. Mọi người sử dụng cả hai từ có thể thay thế cho nhau, nhưng nó có nghĩa là điều tương tự. Đó là một phân phối xác suất liên tục.Normal Distribution is also known as a Gaussian distribution or famously Bell Curve. People use both words interchangeably, but it means the same thing. It is a continuous probability distribution.

Hàm mật độ xác suất (PDF) cho phân phối bình thường:

Trong đó, = trung bình, = độ lệch chuẩn, x = giá trị đầu vào.

Terminology:

• Trung bình - giá trị trung bình là trung bình thông thường. Tổng của tổng số điểm chia cho tổng số điểm. – The mean is the usual average. The sum of total points divided by the total number of points.
• Độ lệch chuẩn - Độ lệch chuẩn cho chúng ta biết cách truyền bá dữ liệu của Google như thế nào. Đó là thước đo giá trị quan sát được bao xa từ giá trị trung bình. – Standard deviation tells us how “spread out” the data is. It is a measure of how far each observed value is from the mean.

Có vẻ khó khăn, phải không? Nhưng nó rất đơn giản.

1. Ví dụ thực hiện phân phối bình thường

Hãy để một cái nhìn vào mã bên dưới. Chúng tôi sẽ sử dụng Numpy và Matplotlib cho cuộc biểu tình này:

# Importing required libraries

import numpy as np
import matplotlib.pyplot as plt

# Creating a series of data of in range of 1-50.
x = np.linspace(1,50,200)

#Creating a Function.
def normal_dist(x , mean , sd):
    prob_density = (np.pi*sd) * np.exp(-0.5*((x-mean)/sd)**2)
    return prob_density

#Calculate mean and Standard deviation.
mean = np.mean(x)
sd = np.std(x)

#Apply function to the data.
pdf = normal_dist(x,mean,sd)

#Plotting the Results
plt.plot(x,pdf , color = 'red')
plt.xlabel('Data points')
plt.ylabel('Probability Density')

2. Thuộc tính phân phối bình thường

Hàm mật độ phân phối bình thường chỉ đơn giản là chấp nhận một điểm dữ liệu cùng với giá trị trung bình và độ lệch chuẩn và ném một giá trị mà chúng ta gọi là mật độ xác suất.probability density.

Chúng ta có thể thay đổi hình dạng của đường cong chuông bằng cách thay đổi độ lệch trung bình và độ lệch chuẩn.

Thay đổi giá trị trung bình sẽ chuyển đường cong theo giá trị trung bình đó, điều này có nghĩa là chúng ta có thể thay đổi vị trí của đường cong bằng cách thay đổi giá trị trung bình trong khi hình dạng của đường cong vẫn còn nguyên vẹn.

Hình dạng của đường cong có thể được kiểm soát bởi giá trị của độ lệch chuẩn. Một độ lệch chuẩn nhỏ hơn sẽ dẫn đến một đường cong giới hạn chặt chẽ trong khi giá trị cao sẽ dẫn đến một đường cong trải rộng hơn.

Một số thuộc tính tuyệt vời của phân phối bình thường:

• Trung bình, chế độ và trung bình đều bằng nhau.
• Tổng diện tích dưới đường cong bằng 1.
• Đường cong là đối xứng xung quanh trung bình.

Quy tắc thực nghiệm cho chúng ta biết rằng:

• 68% dữ liệu nằm trong một độ lệch chuẩn của giá trị trung bình.
• 95% dữ liệu nằm trong hai độ lệch chuẩn của giá trị trung bình.
• 99,7% dữ liệu nằm trong ba độ lệch chuẩn của giá trị trung bình.

Đó là một trong những phân phối quan trọng nhất trong tất cả các số liệu thống kê. Sự phân phối bình thường là kỳ diệu vì hầu hết các hiện tượng xảy ra tự nhiên tuân theo phân phối bình thường. Ví dụ, huyết áp, điểm IQ, chiều cao theo phân phối bình thường.

Tính toán xác suất với phân phối bình thường

Để tìm xác suất của một giá trị xảy ra trong một phạm vi trong phân phối bình thường, chúng ta chỉ cần tìm khu vực dưới đường cong trong phạm vi đó. tức là chúng ta cần tích hợp hàm mật độ.

Vì phân phối bình thường là một phân phối liên tục, khu vực dưới đường cong thể hiện xác suất.

Trước khi nhận được chi tiết trước tiên, hãy để Lừa chỉ biết phân phối bình thường tiêu chuẩn là gì.

Một phân phối bình thường tiêu chuẩn tương tự như phân phối bình thường với trung bình = 0 và độ lệch chuẩn = 1. is just similar to a normal distribution with mean = 0 and standard deviation = 1.

Z = (x-μ)/ σ

Giá trị Z ở trên còn được gọi là điểm Z. Một điểm Z cho bạn một ý tưởng về mức độ trung bình của một điểm dữ liệu. z-score. A z-score gives you an idea of how far from the mean a data point is.

Nếu chúng ta dự định tính toán xác suất theo cách thủ công, chúng ta sẽ cần tra cứu giá trị Z của chúng ta trong bảng Z để xem giá trị tỷ lệ phần trăm tích lũy. Python cung cấp cho chúng tôi các mô -đun để thực hiện công việc này cho chúng tôi. Hãy để vào đó.

1. Tạo đường cong bình thường

Chúng tôi sẽ sử dụng chức năng lớp

# import required libraries
from scipy.stats import norm
import numpy as np
import matplotlib.pyplot as plt
import seaborn as sb

# Creating the distribution
data = np.arange(1,10,0.01)
pdf = norm.pdf(data , loc = 5.3 , scale = 1 )

#Visualizing the distribution

sb.set_style('whitegrid')
sb.lineplot(data, pdf , color = 'black')
plt.xlabel('Heights')
plt.ylabel('Probability Density')

Giả sử chúng ta có dữ liệu về chiều cao của người lớn trong một thị trấn và dữ liệu theo phân phối bình thường, chúng ta có kích thước mẫu đủ với trung bình bằng 5.3 và độ lệch chuẩn là 1.

Thông tin này là đủ để tạo một đường cong bình thường.

# import required libraries
from scipy.stats import norm
import numpy as np
import matplotlib.pyplot as plt
import seaborn as sb

# Creating the distribution
data = np.arange(1,10,0.01)
pdf = norm.pdf(data , loc = 5.3 , scale = 1 )

#Visualizing the distribution

sb.set_style('whitegrid')
sb.lineplot(data, pdf , color = 'black')
plt.xlabel('Heights')
plt.ylabel('Probability Density')

Phương pháp lớp

# import required libraries
from scipy.stats import norm
import numpy as np
import matplotlib.pyplot as plt
import seaborn as sb

# Creating the distribution
data = np.arange(1,10,0.01)
pdf = norm.pdf(data , loc = 5.3 , scale = 1 )

#Visualizing the distribution

sb.set_style('whitegrid')
sb.lineplot(data, pdf , color = 'black')
plt.xlabel('Heights')
plt.ylabel('Probability Density')

# import required libraries
from scipy.stats import norm
import numpy as np
import matplotlib.pyplot as plt
import seaborn as sb

# Creating the distribution
data = np.arange(1,10,0.01)
pdf = norm.pdf(data , loc = 5.3 , scale = 1 )

#Visualizing the distribution

sb.set_style('whitegrid')
sb.lineplot(data, pdf , color = 'black')
plt.xlabel('Heights')
plt.ylabel('Probability Density')

# import required libraries
from scipy.stats import norm
import numpy as np
import matplotlib.pyplot as plt
import seaborn as sb

# Creating the distribution
data = np.arange(1,10,0.01)
pdf = norm.pdf(data , loc = 5.3 , scale = 1 )

#Visualizing the distribution

sb.set_style('whitegrid')
sb.lineplot(data, pdf , color = 'black')
plt.xlabel('Heights')
plt.ylabel('Probability Density')

# import required libraries
from scipy.stats import norm
import numpy as np
import matplotlib.pyplot as plt
import seaborn as sb

# Creating the distribution
data = np.arange(1,10,0.01)
pdf = norm.pdf(data , loc = 5.3 , scale = 1 )

#Visualizing the distribution

sb.set_style('whitegrid')
sb.lineplot(data, pdf , color = 'black')
plt.xlabel('Heights')
plt.ylabel('Probability Density')

# import required libraries
from scipy.stats import norm
import numpy as np
import matplotlib.pyplot as plt
import seaborn as sb

# Creating the distribution
data = np.arange(1,10,0.01)
pdf = norm.pdf(data , loc = 5.3 , scale = 1 )

#Visualizing the distribution

sb.set_style('whitegrid')
sb.lineplot(data, pdf , color = 'black')
plt.xlabel('Heights')
plt.ylabel('Probability Density')

2. Tính xác suất của dữ liệu cụ thể xảy ra

Bây giờ, nếu chúng tôi được yêu cầu chọn một người ngẫu nhiên từ phân phối này, thì xác suất chiều cao của người sẽ nhỏ hơn 4,5 ft là bao nhiêu?

Khu vực dưới đường cong như trong hình trên sẽ là xác suất chiều cao của người sẽ nhỏ hơn 4,5 ft nếu được chọn ngẫu nhiên từ phân phối. Hãy để xem cách chúng ta có thể tính toán điều này trong Python.

Vùng dưới đường cong không là gì ngoài việc tích hợp hàm mật độ với giới hạn bằng -∞ đến 4,5.

norm(loc = 5.3 , scale = 1).cdf(4.5)

0.211855 or 21.185 %

Dòng mã ở trên cho thấy xác suất có 21,18% cơ hội là nếu một người được chọn ngẫu nhiên từ phân phối bình thường với giá trị trung bình là 5,3 và độ lệch chuẩn là 1, thì chiều cao của người sẽ dưới 4,5 ft .

Chúng tôi khởi tạo đối tượng của lớp

# import required libraries
from scipy.stats import norm
import numpy as np
import matplotlib.pyplot as plt
import seaborn as sb

# Creating the distribution
data = np.arange(1,10,0.01)
pdf = norm.pdf(data , loc = 5.3 , scale = 1 )

#Visualizing the distribution

sb.set_style('whitegrid')
sb.lineplot(data, pdf , color = 'black')
plt.xlabel('Heights')
plt.ylabel('Probability Density')

# import required libraries
from scipy.stats import norm
import numpy as np
import matplotlib.pyplot as plt
import seaborn as sb

# Creating the distribution
data = np.arange(1,10,0.01)
pdf = norm.pdf(data , loc = 5.3 , scale = 1 )

#Visualizing the distribution

sb.set_style('whitegrid')
sb.lineplot(data, pdf , color = 'black')
plt.xlabel('Heights')
plt.ylabel('Probability Density')

Giá trị xác suất tích lũy từ -∞ đến ∞ sẽ bằng 1.

Bây giờ, một lần nữa chúng tôi được yêu cầu chọn một người ngẫu nhiên từ phân phối này, thì xác suất chiều cao của người sẽ ở trong khoảng từ 6,5 đến 4,5 ft là bao nhiêu?

cdf_upper_limit = norm(loc = 5.3 , scale = 1).cdf(6.5)
cdf_lower_limit = norm(loc = 5.3 , scale = 1).cdf(4.5)

prob = cdf_upper_limit - cdf_lower_limit
print(prob)

0.673074 or 67.30 %

Mã trên đầu tiên tính toán giá trị xác suất tích lũy từ -∞ đến 6,5 và sau đó là giá trị xác suất tích lũy từ -∞ đến 4.5. Nếu chúng ta trừ CDF là 4,5 từ CDF là 6.5, kết quả chúng ta nhận được là diện tích dưới đường cong giữa giới hạn 6,5 và 4,5.

Bây giờ, điều gì sẽ xảy ra nếu chúng ta được hỏi về xác suất rằng chiều cao của một người được chọn ngẫu nhiên sẽ trên 6,5ft?

Thật đơn giản, như chúng ta biết tổng diện tích dưới đường cong bằng 1 và nếu chúng ta tính toán giá trị xác suất tích lũy từ -∞ đến 6,5 và trừ nó từ 1, kết quả sẽ là xác suất mà chiều cao của một người được chọn ngẫu nhiên sẽ là Trên 6,5ft.

cdf_value = norm(loc = 5.3 , scale = 1).cdf(6.5)
prob = 1- cdf_value
print(prob)

0.115069 or 11.50 %.

Đó là rất nhiều thứ để chìm vào, nhưng tôi khuyến khích tất cả tiếp tục thực hành khái niệm thiết yếu này cùng với việc thực hiện bằng cách sử dụng Python.

Mã hoàn chỉnh từ việc triển khai trên:

# import required libraries
from scipy.stats import norm
import numpy as np
import matplotlib.pyplot as plt
import seaborn as sb

# Creating the distribution
data = np.arange(1,10,0.01)
pdf = norm.pdf(data , loc = 5.3 , scale = 1 )

#Probability of height to be under 4.5 ft.
prob_1 = norm(loc = 5.3 , scale = 1).cdf(4.5)
print(prob_1)

#probability that the height of the person will be between 6.5 and 4.5 ft.

cdf_upper_limit = norm(loc = 5.3 , scale = 1).cdf(6.5)
cdf_lower_limit = norm(loc = 5.3 , scale = 1).cdf(4.5)

prob_2 = cdf_upper_limit - cdf_lower_limit
print(prob_2)

#probability that the height of a person chosen randomly will be above 6.5ft

cdf_value = norm(loc = 5.3 , scale = 1).cdf(6.5)
prob_3 = 1- cdf_value
print(prob_3)

Sự kết luận

Trong bài viết này, chúng tôi đã có một số ý tưởng về phân phối bình thường, một đường cong bình thường trông như thế nào và quan trọng nhất là việc thực hiện nó trong Python.

Học hỏi!