找到方程系統的根

我想用multiroot命令在rootSolve包中找到b和m。代碼如下。我嘗試了不同的起始值，但結果是NaN或NaNs produced。找到方程系統的根

n <- 23 
model <- function(theta){ 
    b <- theta[1] 
    m <- theta[2] 
    power <- exp(-(x-m)*b) 
    a <- -n/sum(log(1-power)) 
    betat <- apply(x,1,function(x) (x-m)*power/(1-power)) 
    mut <- apply(x,1, function(x) power/(1-power)) 
    F1 <- n/b-sum(x)+n*m+(a-1)*sum(betat) 
    F2 <- n*b-b*(a-1)*sum(mut) 
    c(F1=F1,F2=F2) 
} 
multiroot(f = model, start = c(.5, .5))

那麼有人可以解釋我哪裏出錯了嗎？

來源

2016-01-20 mert

你肯定有辦法解決嗎？除了這個非線性問題通常很難。他們可能對起始值，收斂標準非常明智... – jogo

我添加了系統並報告瞭解決方案 – mert

發佈時有錯誤嗎？ – mert

library(rootSolve) 
x<- c(17.88,28.92,33,41.52,42.12,45.6,48.40,51.84,51.96,54.12,55.56,67.80, 
    68.64,68.64,68.88,84.12,93.12, 98.64,105.12,105.84,127.92,128.04,173.4) 
n <- length(x) 

model <- function(theta){ 
    b <- theta[1] 
    m <- theta[2] 
    power <- exp(-(x-m)*b) 
    a <- -n/sum(log(1-power)) 
    F1 <- n/b-sum(x-m) + (a-1)*sum((x-m)*power/(1-power)) 
    F2 <- n*b - b*(a-1)*sum(power/(1-power)) 
    c(F1=F1,F2=F2) 
} 
# model(c(b = 0.031, m = 4.748)) 
multiroot(f = model, start = c(.03, 5))

所以結果是：

> multiroot(f = model, start = c(.03, 5)) 
$root 
[1] 0.03140027 4.55976021 

$f.root 
      F1   F2 
-2.046363e-12 -6.217249e-15 

$iter 
[1] 5 

$estim.precis 
[1] 1.02629e-12

來源

2016-01-20 15:30:44 jogo

謝謝。但仍然存在問題。初始值的小變化給了NaN。其實我檢查了這個結果，並將它們與原始紙張進行比較。這個結果比那些更好。但開始的價值觀確實是一個大問題。 – mert

如果沒有其他幫助必須對起始值進行網格搜索。 – jogo

找到方程系統的根

回答

相關問題