mathematical optimization - R DEoptim() function: How to select parameters to be opimised? -

-

i wish estimate parameters following example:

but how can make deoptim optimise 2 of 3 parameters? - there direct method this?

rm(list=ls()) t <- seq(0.1,20,length=100)  hobs <- 20 + 8*exp(-0.05*t) hsim <- function(p,t) {p[1] + p[2]*exp(-p[3]*t)} upper <- c(30,10,1) lower <- -upper  resfun <- function(p, t, hobs) {           r <- hobs - hsim(p,t)           return(t(r)%*%r) }  deoptim(resfun, lower, upper, hobs = hobs, t = t,                        deoptim.control(np = 80, itermax = 200, f = 1.2, cr = 0.7, trace =false))

i found workaround cf below. there more elegant way of doing this?

rm(list=ls()) t <- seq(0.1,20,length=100)  hobs <- 20 + 8*exp(-0.05*t) #hkorr <- rnorm(100,0,0.2) hsim <- function(p,t) {p[1] + p[2]*exp(-p[3]*t)} upper <- c(20,10,1) lower <- c(1,1,0.001) sel = c(0,1,1) ini = c(20,na,na) # correct upper , lower selected parameters upper <- upper[which(sel==1)] lower <- lower[which(sel==1)]  resfun <- function(par,t,hobs, sel, ini) {           p <- rep(na,3)           p[which(sel == 0)] <- ini[which(sel==0)]           p[which(sel == 1)] <- par           r <- hobs - hsim(p,t)           return(t(r)%*%r) } deoptim(resfun, lower, upper, hobs = hobs, t = t, sel = sel, ini = ini,         deoptim.control(np = 80, itermax = 200, f = 1.2, cr = 0.7, trace =false))

you can set upper , lower parameter value equal whichever parameter want hold constant.

# hold first parameter constant @ 20 opt <- deoptim(resfun, c(20, -10, -1), c(20, 10 ,1), hobs=hobs, t=t,                deoptim.control(np=80, itermax=200, f=1.2, cr=0.7))

Comments

Post a Comment

Popular posts from this blog

php - Magento - Deleted Base url key -

-
i working on cleaning url: http://mydomain.com/shop/shop/deals i going through settings in configuration tab in admin panel on magento , change secure , unsecure base urls remove first base url "shop" web page read: http://mydomain.com/shop/deals so changed unsecure base urls mydomain.com/shop mydomain.com , secure 12.345.67.89/shop 12.345.67.89. wrong thing do. i cannot reload page or admin console. how access these files through ftp , change correctly there?? when make changes in admin panel changes made in database not in file. have make changes in database. 1.go database , select core_config_data table under table find web/unsecure/base_url , web/secure/base_url . make changes here. don't forget clear cache in var/cache folder
Read more

javascript - Tooltipster plugin not firing jquery function when button or any click even occur -

-
i want use tooltipster plugin click element in body show tooltip , want have close button inside tooltip close it. however use jquery click() function handle won't fired. tried solution in post. worked when tooltip triggered hover event. tooltipster plugin not firing jquery function when link in in tooltip clicked the original solution using hover in jsfiddle trigger: 'hover', instead using click show tooltip trigger: 'click', http://jsfiddle.net/bcqyl/7/ it won't fire click event in jquery block. can captured normal javascript function only. checked should code in tooltipster plugin locked , captured "click()" event inside tips when in click trigger mode. i tried other event onchange event of radio button fired when using code in original solution @evan right. tooltip doesn't exist in dom , can't bind/delegate action (even if go way $(document) after domready). if you're after close button inside of too...
Read more

java - WrongTypeOfReturnValue exception thrown when unit testing using mockito -

-
my test list<person> mylist; @test public void testisvalidperson() { mylist = new arraylist<person>(); mylist.add(new person("tom")); when(persondao.get(person)).thenreturn(mylist); when((persondao.get(person)).isempty()).thenreturn(false);//------exception thrown boolean result = service.isvalid("tom"); assertfalse(result); } method tested: public boolean isvalid(string person){ persondao = new persondao(); person personobj = new person(person); return (persondao.get(person).isempty())?false : true; } exception thrown: org.mockito.exceptions.misusing.wrongtypeofreturnvalue: boolean cannot returned get() get() should return list *** if you're unsure why you're getting above error read on. due nature of syntax above problem might occur because: 1. exception *might* occur in wrongly written multi-threaded tests. please refer mockito faq on limitations of concurrency testing. 2. sp...
Read more