Macc0927 2021-01-22 00:16
浏览 7

R语言 EM 伽马 建模时报错了 请教各位大神,问题出在那里?

df是包含TP1NP一个变量的数据框

tsh.hist <- hist(df1$TP1NP, breaks = seq(0,max(df1$TP1NP) + 0.1, 0.1), plot = FALSE)
tsh.df <- data.frame(tsh = tsh.hist$breaks[2:length(tsh.hist$breaks)],freq = tsh.hist$counts)
tsh.df <- as.mixdata(tsh.df)
fitg <- mix(tsh.df,
            mixparam(c(0.1,2,6),c(0.1,1,2)),
            constr = list(
              conpi = "NONE",
              conmu = "MFX",
              consigma = "SFX",
              fixpi = NULL,
              fixmu = c(FALSE,FALSE, FALSE),
              fixsigma = c(FALSE,FALSE, FALSE)),
            "gamma", emsteps = 1)

fitg.coef <- coef(fitg, natpar = TRUE) -> tshg

报错:Error in nlm(mixlike, lmixdat = mixdat, lmixpar = fitpar, ldist = dist,  :    参数里不能有遗漏值 

  • 写回答

0条回答 默认 最新

    报告相同问题?

    悬赏问题

    • ¥15 shape_predictor_68_face_landmarks.dat
    • ¥15 slam rangenet++配置
    • ¥15 有没有研究水声通信方面的帮我改俩matlab代码
    • ¥15 对于相关问题的求解与代码
    • ¥15 ubuntu子系统密码忘记
    • ¥15 信号傅里叶变换在matlab上遇到的小问题请求帮助
    • ¥15 保护模式-系统加载-段寄存器
    • ¥15 电脑桌面设定一个区域禁止鼠标操作
    • ¥15 求NPF226060磁芯的详细资料
    • ¥15 使用R语言marginaleffects包进行边际效应图绘制