我正在尝试使用 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)
但是我在提供起始值时出错。我不知道如何在我的情况下提供起始值。