在R中使用`optim（）`解決兩個參數？

Question

我試圖找到shape1和shape2在dbeta()使得從dbeta()兩個輸入值答案：.6和.8成爲3。

我使用下面的optim()但沒有得到確切的結果，這是我希望得到值shape1和shape2與.6和.8使用時給3和3，但他們不這樣做，爲什麼呢？

f <- function(x) { 
y <- c(3, 3) - dbeta(c(.6, .8), shape1 = x[1], shape2 = x[2]) 
} 

AA = optim(c(1, 1), function(x) sum(f(x)^2), control = list(reltol = (.Machine$double.eps))) 

parms = unname(AA$par) 

dbeta(c(.6, .8), parms[1], parms[2]) # Here I expect to get `3` for `.6` and `.8` but I don't. 

Answer 1

我簡要介紹了一下。我不認爲有與合適的任何問題：這裏的可能性表面的照片：

library(emdbook) 
cc <- curve3d(g(c(x,y)),xlim=c(1,20),ylim=c(1,20), 
       sys3d="none") 
pp <- which(cc$z==min(cc$z),arr.ind=TRUE) 
png("betasurf.png") 
with(cc,image(x,y,z)) 
points(parms[1],parms[2],pch=16) 
points(cc$x[pp[1]],cc$y[pp[2]],pch=1) 
dev.off() 

實心圓是擬合值，空心圓是電網的最低;我認爲差異只是數字模糊（我放大了幾次以確保）。無論如何，沒有任何證據表明像多個optima那樣奇怪。

我認爲這個問題是您與他們建立了一套對點的那不能同時通過匹配任何 Beta分佈; optim()是給你最佳的契合......

png("betatmp.png") 
curve(dbeta(x,parms[1],parms[2]),from=0,to=1) 
points(c(0.6,0.8),c(3,3),pch=16) 
dev.off()