我想使用FinAsym packagage估算PIN码。我有两个数据集,结构和类都相同。 这个工作正常:
dput(ms)
structure(c(4, 5, 4), .Dim = c(3L, 1L))
par0 <- c(0.5, 0.5, 0.5, 0.5)
param_optim <- optim(par0, pin_likelihood, gr=NULL, ms)
alph <- c( param_optim$par[3] )
miu <- c( param_optim$par[2] )
delt <- c( param_optim$par[4] )
epsi <- c( param_optim$par[1] )
PSI <- c( ( 1-miu )*epsi + miu )
pin <- c( miu/PSI )
并提供所需的输出
dput(pin)
5.85127760950505e-07
但是这个不起作用
dput(m2)
structure(c(580, 2163, 2185), .Dim = c(3L, 1L))
par0 <- c(0.5, 0.5, 0.5, 0.5)
param_optim <- optim(par0, pin_likelihood, gr=NULL, m2)
alph <- c( param_optim$par[3] )
miu <- c( param_optim$par[2] )
delt <- c( param_optim$par[4] )
epsi <- c( param_optim$par[1] )
PSI <- c( ( 1-miu )*epsi + miu )
pin <- c( miu/PSI )
并给出了以下错误:
Error in optim(par0, pin_likelihood, gr = NULL, m2) :
function cannot be evaluated at initial parameters
我很困惑,因为两者在结构和类别上都相同,那么为什么我会收到此错误?任何帮助表示赞赏。