襷 糾 螻襴 譯殊伎 曙＾蟇 伎 覈襯 豕る 覦 牛 企讌. 蠏覿 OLS 谿 螻燕 豕る 覦覯願, ML 磯襯 豕る 蟆企. 蠍一 谿螻燕企 磯螳 覈螳 螻 豕 豕螳 豕 企麹.
optim 企蟆 覈襯 豕貅譯朱 企. 螳 襯 危ｴ覲伎.
> x <- rnorm(100)
> y <- 3*x+2+rnorm(100)
> lm(y~x)
Call:
lm(formula = y ~ x)
Coefficients:
(Intercept) x
1.959 2.908
蠏 0, 譴ｸ谿 1 蠏覿襦覿 100螳 螳 れ豢伎 覲 x 麹. 蠏碁Μ螻 x 蠍一瑚鍵 3 螻燕螻 ｸ 2襯 れ れ姶襯 y襯 燕. OLS 蠏覿 れ 螳螻 蟇一 觜訣蟆 蠍一瑚鍵 2.908, ｸ 1.959 谿場. (れ豢襦 襷 一危一企襦 れ 襷 蟲豌伎 螳 譟郁 るゼ .)
願螻 螳 optim 伎伎 企慨. 一 OLS 襯 襷れ企慨. 覈 豢 朱誤磯れ 覯″磯 覦れ碁.
ols <- function(b){
t <- b[1] + b[2]*x
sum( (y - t)^2 )
}
伎 麹 豐蠍郁螻 蟷 覈襯 optim 蟆譯朱 . optim 覈襯 豕る襦 豕貅狩 蟆曙一 覈 -襯 覿譴. 豐蠍郁 覯″一 蠍語企 朱誤一 螳 螳 .
> optim(c(0,1), ols)
$par
[1] 1.958803 2.907654
$value
[1] 89.22844
$counts
function gradient
69 NA
$convergence
[1] 0
$message
NULL
par 豢 朱誤一 螳 碁. lm 襦 れ OLS 蠏覿螻 螳 蟆 . value 豕 覈 螳企. 谿螻燕 89.22. count 語 襯 企 69覯襷 豕螳 る 蟆 . convergence螳 0企 企る 蟆 碁. 蠍壱 誤 伎 襷 谿瑚.