联立方程式 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去做 估计出来的参数也一样
所以没什么意义 可能是我理解错误
不知道这个函数的用意是什么
谢谢