拓端tecdat|R語言最優(yōu)化問題中的共軛函數(shù)
原文鏈接:http://tecdat.cn/?p=18993
原文出處:拓端數(shù)據(jù)部落公眾號
在回歸模型研究中,我們將討論優(yōu)化,而經(jīng)典工具就是所謂的共軛。給定函數(shù)f:Rp→R,其共軛值為函數(shù)f ?:Rp→R使得

可視化考慮一個簡單的拋物線函數(shù)(在維度1中)f(x)= x ^ 2 / 2,然后f ?(2)是線x?2x與函數(shù)f(x)之間的最大距離。
f = function(x) x^2/2
fstar = function(y) max(y*x-vf)
我們可以在下圖上看到。
polygon(c(x[idx2],rev(x[idx2])),c(vf[idx2],rev(x0*x[idx2],col=rgb(0,1,0,.3,border=NA)
abline(a=0,b=x0,col="red")
segments(x[i],x0*x[i],x[i],f(x[i]),lwd=3,col="red")


在這種情況下,我們實際上可以計算f?,因為

一階條件是x?= y,因此

實際上,對于?p的共軛,我們可以使用以下代碼對其進(jìn)行可視化
f = function(x) abs(x)^p/p
fstar = function(y) max(y*x-vf)
vi(1.5)

f = function(x) abs(x)^p/p
fstar = function(y) max(y*x-vf)
vi(1, YL=c(0,10))

在那種情況下,如果f(x)= ∣x∣則

另一種情況是?

我們可以在下面看到
f = function(x) exp(x)
fstar = function(y) max(y*x-vf)
vi(1,YL=c(-3,3))


最受歡迎的見解
1.用R語言模擬混合制排隊隨機服務(wù)排隊系統(tǒng)
2.R語言中使用排隊論預(yù)測等待時間
3.R語言中實現(xiàn)馬爾可夫鏈蒙特卡羅MCMC模型
4.R語言中的馬爾科夫機制轉(zhuǎn)換(Markov regime switching)模型
5.matlab貝葉斯隱馬爾可夫hmm模型
6.用R語言模擬混合制排隊隨機服務(wù)排隊系統(tǒng)
7.Python基于粒子群優(yōu)化的投資組合優(yōu)化
8.R語言馬爾可夫轉(zhuǎn)換模型研究交通傷亡人數(shù)事故預(yù)測
9.用機器學(xué)習(xí)識別不斷變化的股市狀況——隱馬爾可夫模型的應(yīng)用