groupstats.r

来自「这个包里面含有利用极大似然估计方法来拟合成组数据和条件数据的有限混合分布模型的函」· R 代码 · 共 15 行

## last modified June 2002

groupstats <- function(mixdat) 
{
    m <- nrow(mixdat)
    n <- sum(mixdat[, 2])
    amid <- (c(3 * mixdat[1, 1] - 2 * mixdat[2, 1], mixdat[-m, 
        1]) + c(mixdat[-m, 1], 3 * mixdat[m - 1, 1] - 2 * mixdat[m - 
        2, 1]))/2
    amid[1] <- ifelse(amid[1] <= 0 & mixdat[1, 1] > 0, mixdat[1, 
        1]/2, amid[1])
    data.frame(pi = 1, mu = sum(amid * mixdat[, 2])/n, sigma = sqrt(var(rep(amid, 
        mixdat[, 2]))))
}