R 函数 optim()、nlm()、optimize() 可以用来求函数极值,因此可以用来计算最大似然估计。optimize() 只
能求一元函数极值
正态分布最大似然估计有解析表达式。作为示例,用 R 函数进行数值优化求解。
对数似然函数为:
定义 R 的优化目标函数为上述对数似然函数去掉常数项以后乘以 −2,求其最小值点。目标函数为:
objf.norm1 <- function(theta, x){
mu <
R 函数 optim()、nlm()、optimize() 可以用来求函数极值,因此可以用来计算最大似然估计。optimize() 只
能求一元函数极值
正态分布最大似然估计有解析表达式。作为示例,用 R 函数进行数值优化求解。
对数似然函数为:
定义 R 的优化目标函数为上述对数似然函数去掉常数项以后乘以 −2,求其最小值点。目标函数为:
objf.norm1 <- function(theta, x){
mu <