similar to: Increasing precision of rgenoud solutions

I would like to limit the significant figures of the calibrated parameters determined by genoud() in the Rgenoud package. Below is some example output, where column 1 is model run number, columns 2-7 are the parameter values, and columns 8-12 are model fit statistics. I would like genoud to internally limit parameters to 4 decimal places as shown in this output. It is clear that the function is

Dear everyone, I am trying to run rgenoud on several chips simultaneusly. I used the instructions provided on Jasjeet Sekhon's Homepage (http://sekhon.berkeley.edu/rgenoud/multiple_cpus.html). However, I have the newer version of R (R 2.12) installed - for a 64-bit machine. So, when I tried to install the special version of "snow" from a zip file provided by Jasjeet on his page, R

Dear R user fellows, I would like to ask you about the package 'rgenoud' which is a genetic optimization tool. I ran the function 'genoud' with two variables to be minimized by the following command. result<-genoud(fn,nvars=2,starting.values=c(0.5,0), pop.size=1000, max.generations=10, wait.generations=3) Then, I had the following error message. Error in

Hi R users, "genoud" function of "rgenoud" package will optimize my function. If opt = genoud(fn,2,max=TRUE,starting.value=c(1,10),........) opt$value will give the optimized value of the function, "fn". My problem is from the same opt, can I get the value of the function at the initial parameter values? I need the initial value of the function for

I am trying to use the optimizing function genoud() with the snow package on a couple of i686 machines running Redhat Linux WS4 . I don't know anything about PVM or MPI, so I just followed the directions in snow and rgenoud for the simplest method and started a socket cluster. My function fn for genoud involves updating an input file for a separate numerical model with the latest parameter

I've been using RGenoud for a while and it worked smoothly so far. However I came across a strange problem lately (for me at least...). It fails after the first individual and I get the following error message when I set MemoryMatrix = FALSE : Error in genoud(fn, nvars = 8, max = TRUE, pop.size = 10, : REAL() can only be applied to a 'numeric', not a 'character' when I

Hello! I need to purchase a new box, which I would like to optimize for good R performance. For the record, I will run Fedora Core 5 as and OS, and I wanted to know if anyone has experience with how the following affects R performance: - Is there a big advantage to having a 64-bit CPU over having a 32-bit? - Does an Opteron offer any advantages over an Athlon, and if yes, does it justify an

Hello - and sorry for a possibly stupid question, I'm just starting to learn rgenoud. I am defining a function with 5 parameters (p1, p2, p3, p4a, and p4b) and then want to optimize it using genoud. But I am doing something wrong. Before genoud is even able to run it says: "Error in p2 + 1.2 : 'p2' is missing". I assume I did not specify it right. My code is below. The task

I have been using the rgenoud package for a nonlinear least-squares problem with lots of local minima, and it works very well but takes lots of time. According to the article refrenced in the documentation, the original GENOUD-software by the same authors seems to allow for user-supplied analytical derivatives instead of numerical approximations, which would probably save some time. Does anybody

Dear R group, since I do only have a moderately fast MacBook and I would like to get my results faster than within 72h per run ( :-((( ), I tried Amazon AWS, which offers pretty fast computers via remote access. I don't want to post any code because its several hundreds lines of code, I am not looking for the "optimal" answer, but maybe some suggestions from you if faced similar

Hello list, I am dealing with a noisy function (gradient,hessian not available) with simple boundary constraints (x_i>0). I've tried constrOptim() using nelder mead to minimize it but it is way too slow and the returned results are not satisfying. simulated annealing is so hard to tune and it always crashes R program in my case. I wonder if there are any packages or functions can do

Hello, I am trying to do model autocalibration using the snow and rgenoud packages. The function I want to run in task-parallel fashion across multiple machines is one that pre- and post-processes data and runs an external model code. My problem is that external file I/O is happening only in the master node and not in the slaves. I have followed Jasjeet Sekhon's suggestion to test the

I am trying to use rgenoud and snow with an external model and file I/O. Unlike the typical application of rgenoud and snow, I need to run an external executable and do pre- and post-processing of input and output files for each parameter set generated by genoud(). I'm hoping that someone can suggest improvements to my approach and a work-around for file I/O problems I've encountered when

I see something new and unexpected here. > update.packages() trying URL `http://cran.r-project.org/src/contrib/PACKAGES' Content type `text/plain; charset=iso-8859-1' length 163467 bytes opened URL .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... ......... downloaded

Hello! I need to compare 2 datasets whether they come from the same distribution. I use function ks.boot{Matching}. And what is the confidence level of the p-value, returned by ks.boot function? The code is: set=read.table("http://stella.sai.msu.ru:8080/~gala/data/testsets.csv", header=T,sep=',') set1=set[!is.na(set$set1),'set1']

I am trying to determine whether two samples are identical or not. I'm aware that somebody can use the Kolmogorov-Smirnoff test to compare empirical distributions, but since my samples have ties I'm not sure if I'm getting the right p-values for the comparison. Can the Kolmogorov-Smirnoff test be adjusted for the case when ties exists and are there any functions that already

Does anyone else find that using the Var.calc option (for heteroscedasticity consistent std. errors) in Match() (from the Matching library) slows down computation of the matching estimator by a lot? I don't really understand why when I use this option it slows down so much, but for me it does significantly. I want to use the heteroscedasticity consistent std. errors in my project, but as long

I'm interested in doing multinomial logistic regression with equality constraints on some of the parameter values. For example, with categorical outcomes Y_1 (baseline), Y_2, and Y_3, and covariates X_1 and X_2, I might want to impose the equality constraint that \beta_{2,1} = \beta_{3,2} that is, that the effect of X_1 on the logit of Y_2 is the same as the effect of X_2 on the

Dear group, does anybody know how to export libraries/functions to all nodes when launching snow? I want to use a function from fBasics (dstable) for a rGenoud optimization routine, but I fail "making the function accessible" to the nodes created. I know how it works for variables, I know how it works in snowfall(which cant be used in that case), but I dont know how it culd work in

Hello R users, I am trying to minimize two functions with genoud. It is actually one function with two sets of data, each of them having two unknown variables (called Vcmax and gi) which have the same value in each of the function. They are called f.1 and f.2 in the code below. My objective to minimize the functions in order to get the two variables equal in each of the functions. Furthermore, I