0

我正在尝试使用 depmixS4 来估计面板数据的 HMM。

这是我的代码:

rModels <- list(
    list(
        GLMresponse(formula=total_sa_attributed_orders~1,data=dfsmallsorted,family=gaussian())
    ),
    list(
        GLMresponse(formula=total_sa_attributed_orders~1,data=dfsmallsorted,family=gaussian())
    )
)

transition <- list()
transition[[1]] <- transInit(~weekly_total_email_open_count,nstates=2,data=dfsmallsorted)
transition[[2]] <- transInit(~weekly_total_email_open_count,nstates=2,data=dfsmallsorted)

inMod <- transInit(~1,ns=2,data=data.frame(rep(1,1887)),family=multinomial("identity"))
mod <- makeDepmix(response=rModels,transition=transition,prior=inMod,
ntimes=ntimevector,homogeneous=FALSE)

set.seed(123)
fm1 <- fit(mod)

但是我在提供起始值时出错。我不知道如何在我的情况下提供起始值。

4

0 回答 0