Displaying 4 results from an estimated 4 matches for "bestval".
Did you mean:
testval
2006 Sep 23
1
Fitdistr() versus nls()
...----- BEGINNING OF CODE
----------------------------------------------------------------
cdf.all=read.table("all_failures.cdf", header=FALSE, col.names=c
("ttr", "cdf"), sep=":" )
allvals.x=array(t(cdf.all[1]))
allvals.y=array(t(cdf.all[2]))
library(MASS)
bestval.exp.nls=bestval.exp.fit=-1
plot(allvals.x, allvals.y)
for(it in 1:100){
#extract random samples
random=sort(sample(1:length(allvals.x), 15))
somevals.x=allvals.x[c(random)]
somevals.y=allvals.y[c(random)]
#fit with nls and fitdistr
fit.exp = fitdistr(somevals.y, "exponential")
nls...
2007 Jan 05
1
help for memory problem with 64-bit machines
...t;)
*** caught segfault ***
address 0x47cd37c, cause 'memory not mapped'
Traceback:
1: .Fortran("anneal", as.integer(criterio), as.integer(p),
as.double(as.vector(mat)), as.integer(kmin), as.integer(kmax),
as.double(valores), as.integer(vars), as.double(bestval), as.integer(bestvar),
as.integer(nexclude), as.integer(exc), as.integer(ninclude), as.integer(inc),
as.integer(nsol), as.integer(niter), as.logical(improvement),
as.double(cooling), as.double(temp), as.integer(coolfreq),
as.integer(length(pcindices)), as.integer(pcindices), as.logical(esp),
a...
2012 Feb 10
3
Schwefel Function Optimization
All,
I am looking for an optimization library that does well on something as chaotic as the Schwefel function:
schwefel <- function(x) sum(-x * sin(sqrt(abs(x))))
With these guys, not much luck:
> optim(c(1,1), schwefel)$value
[1] -7.890603
> optim(c(1,1), schwefel, method="SANN", control=list(maxit=10000))$value
[1] -28.02825
> optim(c(1,1), schwefel, lower=c(-500,-500),
2012 Feb 10
3
Schwefel Function Optimization
All,
I am looking for an optimization library that does well on something as chaotic as the Schwefel function:
schwefel <- function(x) sum(-x * sin(sqrt(abs(x))))
With these guys, not much luck:
> optim(c(1,1), schwefel)$value
[1] -7.890603
> optim(c(1,1), schwefel, method="SANN", control=list(maxit=10000))$value
[1] -28.02825
> optim(c(1,1), schwefel, lower=c(-500,-500),