function - Solving for two parameters using `optim()` in R? -


i'm trying find shape1 , shape2 in dbeta() such answer dbeta() 2 input values: .6 , .8 become 3.

i'm using below optim() don't exact result, expect getting values shape1 , shape2 when used .6 , .8 give 3 , 3, don't, why? 

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 expect `3` `.6` , `.8` don't.

i had brief @ this. don't think there's problem fit: here's picture of likelihood surface:

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()

enter image description here

filled circle fitted value, open circle minimum of grid; think difference numerical fuzz (i zoomed in few times make sure). in case, there's no evidence of weird multiple optima.

i believe issue you've set set of pair of points can't simultaneously matched any beta distribution; optim() giving best possible fit ...

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()

enter image description here


Comments

Post a Comment

Popular posts from this blog

PHP and MySQL WP -

i trying accomplish 2 things here code. writing php directly wp page using insert php plugin. 1) query mysql database , return results 1 column 4 columns. 2) make each result own hyperlink directs user new wp page runs new query on page. i able query show results , turn results dynamic hyperlink, cannot formatting down. right returns result 1 column. here have: [insert_php] $servername = "localhost"; $username = "login"; //edited privacy $password = "password"; //edited privacy $database = "database"; //edited privacy // create connection $conn = new mysqli($servername, $username, $password, $database); // check connection if ($conn->connect_error) { die("connection failed: " . $conn->connect_error); } echo "connected successfully"; //get results database $result = mysqli_query($conn,"select mine mines order mine"); $all_property = array(); //declare array saving property //showing proper...
Read more

android - InAppBilling registering BroadcastReceiver in AndroidManifest -

i'm implementing inapp billing in android app. works well, however, i'm trying decouple broadcast receiver activity manifest. suggestion in android's trivialdrive sample: // important: dynamically register broadcast messages updated purchases. // register receiver here instead of <receiver> in manifest // because call getpurchases() @ startup, therefore can ignore // broadcasts sent while app isn't running. // note: registering listener in activity bad idea, done here // because sample. regardless, receiver must registered after // iabhelper setup, before first call getpurchases(). currently there's class extends broadcastreceiver : public class iabbroadcastreceiver extends broadcastreceiver { /** * intent action receiver should filter for. */ public static final string action = "com.android.vending.billing.purchases_updated"; private final iabbroadcastlistener mlistener; public iabbroadcastreceiver(iabbroadcastlistener listener) { m...
Read more

nginx - phpPgAdmin - log in works but I have to login again after clicking on any links -

i have setup postgres 9.6 , phppgadmin on centos 7 server nginx, i can log in phppgadmin, postgresql server on right of page remains crossed out , link try browse (databases, roles etc) requires log in again. i'm thinking seems sessions problem, unsure how resolve it. any pointers on how fix or should looking @ in order track down problem? solved. turns out owner of /var/lib/php/session apache - needed nginx. sorted changing ownership nginx below: chown -r nginx:nginx /var/lib/php/session/
Read more