Question

我想估计 R 中的模型。

其中一部分是由两个 OLS 组成的有限混合模型。

作为R的新生，我不知道如何在R中编写此概率密度函数。

我想知道您是否可以提供帮助。

概率密度函数如下：

f(y|x)=(p/σ1)*φ(y-x*b1/σ1)+((1-p)/σ2)*φ(y-x*b2/σ2)

我已经使用stata编写了一个示例：

 
   gen double f1'=normalden($ML_y1,xb1',exp(lns1'))
   gen doublef2'=normalden($ML_y1,xb2',exp(lns2'))
   tempvar p
   gen double p'=exp(lp')/(1+exp(lp'))
   replacelnf'=ln(p'*f1'+(1-p')*f2')

我想知道您是否可以告诉我如何在R中编写此函数。

非常感谢，期待您的帮助

Answer 1

请参见功能FLXMRglm。密度用dnorm

估算

library(flexmix)
FLXMRglm

# your case 


if (family == "gaussian") {
        z@defineComponent <- function(para) {
            predict <- function(x, ...) {
                dotarg = list(...)
                if ("offset" %in% names(dotarg)) 
                  offset <- dotarg$offset
                p <- x %*% para$coef
                if (!is.null(offset)) 
                  p <- p + offset
                p
            }
            logLik <- function(x, y, ...) dnorm(y, mean = predict(x, 
                ...), sd = para$sigma, log = TRUE)
            new("FLXcomponent", parameters = list(coef = para$coef, 
                sigma = para$sigma), logLik = logLik, predict = predict, 
                df = para$df)
        }
        z@fit <- function(x, y, w, component) {
            fit <- lm.wfit(x, y, w = w, offset = offset)
            z@defineComponent(para = list(coef = coef(fit), df = ncol(x) + 
                1, sigma = sqrt(sum(fit$weights * fit$residuals^2/mean(fit$weights))/(nrow(x) - 
                fit$rank))))
        }
    }

如何在R中为混合模型编写此函数？

1 个答案: