聯立方程式 or system equation
用MLE去做 我目前只會做一條equation的MLE
我目前上網找的只能做到一條equation的MLE估計
目前code如下:
#==========================================
library(stats4)
set.seed(1001)
N <- 100
x <- runif(N)
y <- 5 * x + 3 + rnorm(N)
LL <- function(beta0, beta1, mu, sigma) {
# Find residuals
#
R = y - x * beta1 - beta0
#
# Calculate the likelihood for the residuals (with mu and sigma as parameters)
#
R = suppressWarnings(dnorm(R, mu, sigma))
#
# Sum the log likelihoods for all of the data points
#
-sum(log(R))
}
fit <- mle(LL, start = list(beta0 = 3, beta1 = 1, mu = 0, sigma=1))
summary(fit)
#==========================================
R = y - x * beta1 - beta0
這是我的一條equation
我目前只會做到一條
如果要做一條以上，要如何做?
#=========================================
另外我也有找到這個package "systemfit"
它可以一次做兩個以上的迴歸
不過 我的理解是 他是分開做的
因為我分別用lm去做 估計出來的參數也一樣
所以沒什麼意義 可能是我理解錯誤
不知道這個函數的用意是什麼
謝謝