Trong thống kê, chúng tôi nói rằng hồi quy là tuyến tính khi nó tuyến tính trong các tham số. Khớp các mô hình tuyến tính là một nhiệm vụ dễ dàng, chúng ta có thể sử dụng phương pháp bình phương nhỏ nhất và lấy các tham số tối ưu cho mô hình của mình. Trong Python, bạn có thể đạt được điều này bằng cách sử dụng một loạt các thư viện như
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
Bây giờ, trước khi chúng ta đi vào các mô hình, điều quan trọng cần lưu ý là bối cảnh về lý do và khi nào bạn sẽ sử dụng các loại mô hình này. Chúng không nhất thiết phải thay thế cho tổ hợp của bạn [XGBoost, LightGBM, Random Forest, Neural Network, v.v. ] mô hình học máy. Tốt nhất là nghĩ về các mô hình phi tuyến tính này trong bối cảnh khớp đường cong, tôi. e. bạn muốn mô hình hóa một hiện tượng bằng cách sử dụng một đường cong mà sau này bạn có thể giải thích. Vì vậy, đây là một số cân nhắc về lý do tại sao bạn sẽ sử dụng chúng
khả năng diễn giải. Nếu bạn lập mô hình một sự kiện bằng cách sử dụng các đường cong, thì kết quả của bạn sẽ dễ hiểu hơn nhiều, vì thông thường phương trình toán học có ý nghĩa đối với từng tham số. Quan trọng. điều này không có nghĩa là bạn mất đi sức mạnh của mô hình máy học phức tạp của mình, vì bạn vẫn có thể sử dụng mô hình máy học cùng với đường cong mà bạn đang điều chỉnh [e. g. phân đoạn dữ liệu của bạn bằng mô hình, sau đó điều chỉnh các đường cong có thể hiểu được trên từng phân đoạn]. Điều này sẽ giúp ích không chỉ cho người xây dựng mô hình trong tương lai mà còn cho các bên liên quan trong kinh doanh. Sẽ dễ dàng hơn nhiều để giải thích mô hình của bạn đang làm gì nếu đầu ra có thể hiểu được
Thuộc tính toán học. Đôi khi, điều thực sự quan trọng là đảm bảo một số thuộc tính toán học trên đầu ra mô hình của bạn. Ví dụ: nếu bạn đang lập mô hình
4 của một số loài côn trùng thì có thể dự đoán chiều cao của bạn sẽ tăng một cách đơn điệu khi tuổi tăng lên. Việc mô tả những hiện tượng này bằng các phương trình toán học sẽ cho bạn khả năng áp dụng ràng buộc này trong giai đoạn lập mô hình. Bạn vừa khớp với một hàm có thuộc tính đơn điệu nàyfrom sympy import * a, b, c, x = symbols["a b c x"] expr = a*x**2 + b*x + c diff[expr, a], diff[expr, b], diff[expr, c]
Tuyến tính x Phi tuyến
Điều quan trọng là phải phân biệt rằng khi tôi nói phi tuyến tính, ý tôi là phi tuyến tính trong các tham số. Vì vậy, ví dụ, một phương trình bậc hai $y = ax^2 + bx + c$ là tuyến tính, bởi vì mô hình này là tuyến tính trong các tham số $a, b, c$. Mô hình có thể không tuyến tính trong $x$, nhưng nó vẫn có thể tuyến tính trong các tham số
Để hiểu rõ hơn về các mô hình tuyến tính và phi tuyến tính, hãy xem xét các ví dụ sau
$$ \begin{equation} y = \beta_0 + \beta_1 x \end{equation} $$
$$ \begin{equation} y = \beta_0 [1 + \beta_1]^x \end{equation} $$
$$ \begin{equation} y = \ \beta_0\cdot\sin\left[x^{\beta_1}\right]\ +\ \beta_2\cdot\cos\left[e^{x\beta_3}\right]
Phương trình $[1]$ là một dòng đơn giản và các tham số $\beta_0, \beta_1$ là tuyến tính trên $y$, vì vậy đây là một ví dụ về mô hình tuyến tính. Phương trình $[2]$ là cái mà chúng tôi gọi là công thức Tăng trưởng theo cấp số nhân, trong đó $\beta_0$ thường biểu thị điểm bắt đầu, $x$ là thước đo thời gian và $\beta_1$ [thường được gọi là $r$] biểu thị tốc độ tăng trưởng. Trong trường hợp này, chúng tôi muốn có được điểm xuất phát và tốc độ tăng trưởng;
$$ y = \beta_0 [1 + \beta_1]^x \\\
log[y] = log[\beta_0] + x\cdot log[1 + \
Và sau đó gọi $y'=log[y]$, $\beta_0’ = log[\beta_0]$, $\beta_1'=log[1 + \beta_1]$, bạn có thể viết lại Tăng trưởng theo cấp số nhân thành
$$y’ = \beta_0’ + \beta_1’x$$
Và điều chỉnh OLS [Bình phương nhỏ nhất thông thường] bằng cách sử dụng công thức này, vì đây là mô hình tuyến tính [đây được gọi là mô hình log-tuyến tính]. Cuối cùng, phương trình $[3]$ là một mô hình phi tuyến tính, vì hồi quy là phi tuyến tính trong các tham số $\beta_0, \beta_1, \beta_2, \beta_3, \beta_4$ [Tôi hy vọng vậy, khi tôi tạo nhanh phương trình này
Trong bài viết này, trước tiên tôi sẽ trình bày các công cụ chính mà chúng ta sẽ sử dụng để điều chỉnh các mô hình, sau đó giải thích một loạt các mô hình phi tuyến hữu ích + mã + đồ thị của mô hình cho bất cứ khi nào bạn cần điều chỉnh các phương trình kỳ diệu như vậy. Tôi cũng sẽ không đi sâu hơn vào mục đích và cách sử dụng từng model vì nó không phải trọng tâm của bài viết.
Lắp đường cong trong Python
Một phần lớn của bài viết này dựa trên một blogpost khác, nhằm vào người dùng R và sử dụng các thư viện R với các phương trình được tạo sẵn và các quy trình tích hợp sẵn để tạo điều kiện thuận lợi cho cuộc sống của nhà khoa học dữ liệu. Mặc dù đôi khi tôi sử dụng R, nhưng hầu hết thời gian tôi sử dụng Python, vì vậy thật tuyệt khi giữ cơ sở mã dự án của tôi với số lượng ngôn ngữ lập trình riêng biệt tối thiểu. Đáng buồn thay, không có thư viện dễ sử dụng/phổ biến nào trong Python tích hợp sẵn các mô hình này [nếu bạn biết bất kỳ thư viện nào, vui lòng cho tôi biết. ], nên chúng ta sẽ phải thực hiện các mô hình và sáng tạo trong quá trình điều chỉnh đường cong khi cần thiết, đặc biệt là khi chọn các tham số ban đầu.
Bên cạnh bộ ba thần thánh [Pandas, Matplotlib và Numpy], đây là những thư viện chính mà chúng tôi sẽ sử dụng, cho các mục đích phù hợp
- scipy
- Sympy
Chúng tôi chủ yếu sẽ sử dụng các quy trình tối ưu hóa từ
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
6. Chức năng gọi khi tối ưu hóa, tôi. e. mô hình của bạn. Nó phải có chữ ký tương tự nhưfrom sympy import * a, b, c, x = symbols["a b c x"] expr = a*x**2 + b*x + c diff[expr, a], diff[expr, b], diff[expr, c]
7, trong đófrom sympy import * a, b, c, x = symbols["a b c x"] expr = a*x**2 + b*x + c diff[expr, a], diff[expr, b], diff[expr, c]
8 làfrom sympy import * a, b, c, x = symbols["a b c x"] expr = a*x**2 + b*x + c diff[expr, a], diff[expr, b], diff[expr, c]
8 của bạn vàfrom sympy import * a, b, c, x = symbols["a b c x"] expr = a*x**2 + b*x + c diff[expr, a], diff[expr, b], diff[expr, c]
00 là các tham số khác. Có những chữ ký hợp lệ khác, nhưng hãy bắt đầu với chữ ký đơn giản nàyfrom sympy import * a, b, c, x = symbols["a b c x"] expr = a*x**2 + b*x + c diff[expr, a], diff[expr, b], diff[expr, c]
01. Không bắt buộc. Hàm trả về Jacobian của hàm của bạnfrom sympy import * a, b, c, x = symbols["a b c x"] expr = a*x**2 + b*x + c diff[expr, a], diff[expr, b], diff[expr, c]
6. Đôi khi, nó làm cho việc tối ưu hóa trở nên dễ dàng hơn vì trình tối ưu hóa sẽ không phải tính toán gần đúng độ dốc của hàm của bạn bằng sốfrom sympy import * a, b, c, x = symbols["a b c x"] expr = a*x**2 + b*x + c diff[expr, a], diff[expr, b], diff[expr, c]
03. Các giá trịfrom sympy import * a, b, c, x = symbols["a b c x"] expr = a*x**2 + b*x + c diff[expr, a], diff[expr, b], diff[expr, c]
8from sympy import * a, b, c, x = symbols["a b c x"] expr = a*x**2 + b*x + c diff[expr, a], diff[expr, b], diff[expr, c]
05. Các giá trịfrom sympy import * a, b, c, x = symbols["a b c x"] expr = a*x**2 + b*x + c diff[expr, a], diff[expr, b], diff[expr, c]
06from sympy import * a, b, c, x = symbols["a b c x"] expr = a*x**2 + b*x + c diff[expr, a], diff[expr, b], diff[expr, c]
07. Tùy chọn [nhưng rất khuyến khích cho các mô hình phi tuyến tính]. Dự đoán ban đầu cho các tham sốfrom sympy import * a, b, c, x = symbols["a b c x"] expr = a*x**2 + b*x + c diff[expr, a], diff[expr, b], diff[expr, c]
08. Không bắt buộc. Các ràng buộc đối với các giá trị có thể có của từng tham số [hữu ích để kiểm soát hành vi chung của quá trình tối ưu hóa]from sympy import * a, b, c, x = symbols["a b c x"] expr = a*x**2 + b*x + c diff[expr, a], diff[expr, b], diff[expr, c]
09. Không bắt buộc. Số lần đánh giá chức năng tối đa [về cơ bản, bao nhiêu lần lặp lại]from sympy import * a, b, c, x = symbols["a b c x"] expr = a*x**2 + b*x + c diff[expr, a], diff[expr, b], diff[expr, c]
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
Các mô hình này cũng yêu cầu kiến thức toán học khá tốt [phép tính và hàm cơ bản], để bạn biết mình đang làm gì. Cuối cùng, tôi cũng sẽ sử dụng rộng rãi công cụ tuyệt vời Desmos để vẽ các phương trình của chúng ta. Mã đầy đủ sẽ có sẵn trong repo Github được liên kết ở cuối bài viết
Đoán ban đầu và Jacobian
Đối với các mô hình phi tuyến tính, đôi khi điều quan trọng là chúng tôi cung cấp cho trình tối ưu hóa dự đoán ban đầu tốt [
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
Đây là một ví dụ về sự không phù hợp, khi bạn không chuyển tham số
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
Trong đoạn mã bên dưới, tôi cũng sẽ cung cấp hàm Jacobian cho từng mô hình. Chúng không cần thiết cho sự phù hợp, vì
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
Đây là một ví dụ về cách tính đạo hàm riêng của một hàm bằng cách sử dụng
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
$[x^2, x, 1]$
Mô hình lồi/lõm
QUAN TRỌNG. Không phải tất cả các mô hình này đều phi tuyến tính;
Phân rã theo cấp số nhân
Các mô hình phân rã hàm mũ có nhiều ứng dụng khác nhau cho hóa học [phân rã chất], sinh học, kinh tế lượng, v.v.
Tham số hóa chung
$$Y = a\cdot e^{k\cdot X}$$
giải thích thông số
- $a$ là giá trị của $Y$ khi $X = 0$
- $k$ là mức tăng hoặc giảm tương đối của $Y$ khi tăng một đơn vị $X$
Ở đây chúng tôi sử dụng dữ liệu về phản ứng hóa học phân hủy chất, tôi. e. khi thời gian trôi qua nồng độ của một số chất giảm khi nó xuống cấp. Chúng tôi muốn tìm tốc độ phân hủy chất theo thời gian. Như chúng ta đã thấy trước đó, nếu chúng ta chỉ chuyển hàm cho
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
Đây là lý do tại sao điều quan trọng là bạn phải biết chính xác những gì bạn đang cố gắng lập mô hình và những thông số đại diện cho điều gì. Phù hợp thì Dự đoán sẽ không hoạt động =]
Từ việc xem xét dữ liệu và tham số hóa của chúng tôi, chúng tôi có thể đoán rằng
- giá trị của $a$ phải gần với
79from sympy import * a, b, c, x = symbols["a b c x"] expr = a*x**2 + b*x + c diff[expr, a], diff[expr, b], diff[expr, c] - giá trị của
20 phải là số âm, vì dữ liệu cho thấy hành vi phân rã theo cấp số nhânfrom sympy import * a, b, c, x = symbols["a b c x"] expr = a*x**2 + b*x + c diff[expr, a], diff[expr, b], diff[expr, c]
Hãy nghĩ về
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
Sau đó, hãy thử khớp với p0=[100, -1e-3] và xem điều gì sẽ xảy ra
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
Tốt hơn nhiều
Phân rã theo cấp số nhân với tiệm cận thấp hơn
Tham số hóa chung
$$Y = b + [a - b]\cdot e^{k\cdot X}$$
giải thích thông số
- $a$ là giá trị của $Y$ khi $X = 0$
- $b$ là đường tiệm cận dưới nơi phân rã sẽ ổn định
- $k$ là mức tăng hoặc giảm tương đối của $Y$ khi tăng một đơn vị $X$
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
Mô hình tiệm cận [Số mũ âm]
Về cơ bản ngược lại với mô hình phân rã theo cấp số nhân. Nó có thể được sử dụng để mô tả các hiện tượng trong đó $Y$ tăng trưởng có giới hạn khi $X$ tiến đến vô cùng
Ví dụ: nếu $Y$ là chiều dài của một con cá và $X$ đại diện cho tuổi của nó, thì chúng tôi mong đợi sự tăng trưởng [tăng chiều dài của cá] khi nó già đi, nhưng chiều dài của nó sẽ không tăng liên tục, nó có thể sẽ ổn định . Ví dụ về loài cá này dựa trên định luật von Bertalanffy. Ngoài ra, chúng tôi không mong đợi con cá có chiều dài $0$ khi nó được sinh ra, vì vậy tham số hóa chung của loại mô hình này cũng có một tham số để kiểm soát điều này.
Tham số hóa chung
$$Y = b - [b - a]\cdot e^{-k\cdot X}$$
giải thích thông số
- $a$ là giá trị của $Y$ khi $X = 0$
- $b$ là đường tiệm cận trên nơi tăng trưởng sẽ ổn định [giá trị tối đa có thể đạt được cho $Y$]
- $k$ là mức tăng hoặc giảm tương đối của $Y$ khi tăng một đơn vị $X$
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
Mô hình tiệm cận [bị ràng buộc. bắt đầu từ 0]
Đôi khi mô hình của bạn cần bắt đầu từ 0 [i. e. $a = 0$ từ phương trình trước]. Điều này dẫn đến việc tham số hóa đơn giản hơn, mặc dù có cùng hành vi
Tham số hóa chung
$$Y = b\cdot [1 - e^{-k\cdot X}]$$
giải thích thông số
- $b$ là đường tiệm cận trên nơi tăng trưởng sẽ ổn định [giá trị tối đa có thể đạt được cho $Y$]
- $k$ là mức tăng hoặc giảm tương đối của $Y$ khi tăng một đơn vị $X$
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
hồi quy sức mạnh
Tham số hóa chung
$$Y = a\cdot X^{b}$$
giải thích thông số
- $a$ là giá trị của $Y$ khi $X = 0$
- $b$ là sức mạnh [kiểm soát cách $Y$ tăng so với $X$
Hồi quy lũy thừa tương đương với một đường cong hàm mũ nhưng với logarit của $X$
$$aX^b = a\cdot e^{log[X^b]} = a\cdot e^{b\cdot log[X]}$$
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
Ở đây, tôi đang điều chỉnh mô hình Sức mạnh cho bộ dữ liệu có số lượng loài thực vật theo khu vực lấy mẫu của một số thí nghiệm
- diện tích. tổng diện tích của không gian lấy mẫu để đo thực vật
- số loài. số loài được tìm thấy
Một thuộc tính hữu ích của tham số Power là bạn có thể ràng buộc nó sao cho phù hợp tùy thuộc vào hành vi tăng trưởng mà bạn đang cố gắng mô hình hóa
- $0 < b < 1$. Hình lồi, $Y$ tăng khi $X$ tăng
- $b < 0$. Hình lõm, $Y$ giảm khi $X$ tăng
- $b > 1$. Hình lõm, $Y$ tăng khi $X$ tăng
Dưới đây là một số ví dụ
Đường cong sigmoidal
Đây là hai đường cong chính để kiểm tra khi bạn cần lập mô hình các sự kiện có hình dạng sigmoidal
đường cong hậu cần
Đường cong logistic phải quen thuộc với bất kỳ nhà khoa học dữ liệu nào. Nó bắt nguồn từ hàm phân phối logistic tích lũy và đối xứng quanh một điểm uốn
Tham số hóa chung
$$Y = a + \frac{c - a}{1 + e^{b[X - d]}}$$
giải thích thông số
- $c$ là tiệm cận trên
- $a$ là tiệm cận dưới
- $d$ kiểm soát vị trí của điểm uốn [so với X]
- $b$ độ dốc xung quanh điểm uốn và chúng ta có thể sử dụng nó để kiểm soát hành vi 'tăng trưởng'. Nếu $b < 0$ thì $Y$ tăng khi $X$ tăng và nếu $b > 0$ thì $Y$ giảm khi $X$ tăng
Tham số hóa ở trên là chung, nhưng tùy thuộc vào những gì bạn đang cố gắng điều chỉnh, nó có thể được giảm bớt. Ví dụ: nếu chúng ta thấy $a = 0$, chúng ta có một phương trình đơn giản hơn. Chúng ta cũng có thể buộc $c = 1$ hơn nữa, làm cho phương trình trở nên đơn giản hơn. Tất cả phụ thuộc vào những ràng buộc của vấn đề của bạn
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
Trong ví dụ này, chúng tôi điều chỉnh Đường cong Logistic để lập mô hình Tăng trưởng Thực vật, trong trường hợp này là Củ cải đường. Các biến là
- dae. ngày sau khi cây mọc
- không trọng lượng. Đo trọng lượng [gam trên mét vuông chất khô]
Bộ dữ liệu đã được xuất bản trong bài báo này. Scott, R. K. , Wilcockson, S. J. , & Moisey, F. r. [1979]. Ảnh hưởng của thời gian trừ cỏ đến sinh trưởng và năng suất củ cải đường
Chức năng Gompertz
Hay cụ thể hơn, Đường cong Gompertz. Không phải lúc nào vấn đề chúng tôi đang lập mô hình cũng đối xứng và hàm Gompertz có thể mô hình hóa các “tốc độ tăng trưởng” khác nhau xung quanh điểm uốn. Tôi khuyên bạn nên xem trang Wikipedia cho chức năng này, đây là một số ví dụ về nơi nó có thể được sử dụng
- Mức độ phổ biến của điện thoại di động, trong đó chi phí ban đầu cao [do đó mức độ hấp thụ chậm], sau đó là giai đoạn tăng trưởng nhanh, tiếp theo là mức tiêu thụ chậm lại khi đạt đến mức bão hòa
- Dân số trong một không gian hạn chế, khi tỷ lệ sinh đầu tiên tăng và sau đó chậm lại khi đạt đến giới hạn tài nguyên
- Mô hình hóa sự phát triển của khối u
- Mô hình hóa tác động thị trường trong tài chính và các khoản vay địa phương tổng hợp năng động
Tham số hóa chung
$$Y = a + [c - a]\cdot e^{-e^{b[X - d]}}$$
Các tham số có cách hiểu tương tự như trong Đường cong Logistic
giải thích thông số
- $c$ là tiệm cận trên
- $a$ là tiệm cận dưới
- $d$ kiểm soát vị trí của điểm uốn
- $b$ độ dốc xung quanh điểm uốn và chúng ta có thể sử dụng nó để kiểm soát hành vi 'tăng trưởng'. Nếu $b < 0$ thì $Y$ tăng khi $X$ tăng và nếu $b > 0$ thì $Y$ giảm khi $X$ tăng
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
Đường cong Gompertz có một tham số hóa thay thế đảo ngược mô hình tăng trưởng so với mô hình chung
$$Y = a + [c - a]\cdot \left[1 - e^{-e^{-b[X - d]}}\right]$$
Tôi sẽ không phù hợp với bất cứ thứ gì vì nó khá giống với cái trước
Dưới đây là một ví dụ để hình dung rõ hơn về sự khác biệt giữa các đường cong sygmoidal. Các thông số ở đây là.
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]
Kết luận + Mã
Tôi hy vọng những đoạn mã và ví dụ mã này sẽ hữu ích cho bạn [và cho chính tôi trong tương lai]. Có rất nhiều mô hình phi tuyến tính khác mà tôi không đưa vào, một số trong số đó là
- phương trình Michaelis-Menten
- Đường cong giảm năng suất/mật độ
- đường cong log-hậu cần
- đường cong Weibull
Như đã đề cập trước đó, bài viết này dựa trên cơ sở đó bao gồm nhiều chi tiết hơn về các phương trình trên, mặc dù tất cả mã phù hợp đều nằm trong R. Một số bộ dữ liệu tôi đã sử dụng cũng được lấy từ gói
from sympy import *
a, b, c, x = symbols["a b c x"]
expr = a*x**2 + b*x + c
diff[expr, a], diff[expr, b], diff[expr, c]