similar to: constrained optimization in R

Thanks Thomas and Marc that is what I was looking for. -Jason The DESCRIPTION file for a package lists the license; this is also given in the CRAN entry for each package. -thomas Marc Schwartz <MSchwartz@MedAnalytics.com> 03/30/2004 09:02 AM Please respond to MSchwartz To: Jason.L.Higbee@stls.frb.org cc: R-Help

R: Is there a reason why there isn't a Windows Binary version of RQuantlib on CRAN? Usually when there is no binary, I just source the source code, but this one appears to have various calls and methods and things like that so I'm hesitant to do so. I know there has been a big discussion on why Rmetrics doesn't have source for unix/linux, but that isn't on CRAN. Through

When I load R-WinEdt (>library(RWinEdt), I get the warning: Warning message: multi-argument returns are deprecated in: return(InstallRoot, RWinEdtInstalled) I have upgraded to R 1. 8 on Windows, by copying non-base libraries into the 1.8 library folder and updating the help. I also reinstalled R-WinEdt from the zip file as detailed in RWinEdt ReadMe file using the recommend (A)

R: This stems from my curiosity about the previous thread about a request for glm.nb code. The issue of package licenses was brought up and I was hoping for some clarification on that. Using the function license() or licence() gives info on the license for R, but something like license(MASS) does not give info on the license for the MASS package (perhaps it might be good to expand the

R: I got a warning message when running the cor.test function using both Spearman and Kendall rank correlations saying that the p-value may be incorrect due to ties in the data. My data has 35 obs and one series has 6 pairs of ties. Does anyone know if this would likely have a great effect on the p-values calculated.. The values look good; tau = -0.68 with p-value = 8e-9 and rho = =0.84

R: I have a program that runs a For loop for days and I need to (if possible) pause the program. Any ideas on how to do that? If I use stop, how certain can I be that all the statements executed in the previous iteration of the for loop? Thanks, Jason Higbee Research Associate Federal Reserve Bank of St. Louis E: jason.l.higbee@stls.frb.org [[alternative HTML version deleted]]

R-listers: I may be asking too much from R, but is there a way to use time indexing on a time series object. For instance: > tsobject <- ts(1:12, start =1999, freq = 4) > tsobject Qtr1 Qtr2 Qtr3 Qtr4 1999 1 2 3 4 2000 5 6 7 8 2001 9 10 11 12 > tsobject[1999,Qtr4] Error in NextMethod("[") : Object "Qtr4" not found I would

R-Listers: I am doing a quasi-maximum likelihood estimation and I get a "subscript out of bound" error message, Typically I would think this means that a subscript used in the function is literally out of bounds however I don't think this is the case. All I change in the code is a constant, that is hard-wired in (not data dependent and not parameter dependent), furthermore,

R UseRs: I can't figure out how to append intermediate terms inside a For loop. My use of append(), various indexing, and use of data frames, vectors, matrices has been fruitless. Here's a simplified example of what I'm talking about: i <- 1 for(i in 10) { v[i] <- i/10 } > v [1] 1 NA NA NA NA NA NA NA NA 1 I would like: [1] 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9

_______________________________________________________________________________________ Hi, I was just wondering if there was any available R code that could handle general constrained optimisation problems. At the moment I'm using nlminb and optim, both of which allow box constraints on the parameters, but ideally I'd like to be able to specify more general constraints on the solution

in order to know which production the custumer most like,i design a question as follow : Q:there are six production listed below.according to your preference,the production you like most is_____,the production you secondly like is ____,and the third is_____. productionA productionB productionC productionD productionE productionF when the data is collected. i type in a

Note that package.descripton() is deprecated in R-1.9.0. You have to use packageDescription() instead. Jim James W. MacDonald Affymetrix and cDNA Microarray Core University of Michigan Cancer Center 1500 E. Medical Center Drive 7410 CCGC Ann Arbor MI 48109 734-647-5623 >>> Sundar Dorai-Raj <sundar.dorai-raj at PDF.COM> 03/30/04 10:51AM >>> You can access this from

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible Url : https://stat.ethz.ch/pipermail/r-help/attachments/20080305/80536e8c/attachment.pl

Hi Everyone, I'm trying to minimize a function using constrOptim with the simulated annealing method SANN. If I understand constrOptim well, it basically passes most of its arguments to optim while somehow enforcing the constraints. My problem is, that since SANN does not need gradients, when using optim with SANN, the gr argument of optim is used to specify a function to create the next

Hi All, I'm trying to do the following constrained optimization example. Maximize x1*(1-x1) + x2*(1-x2) + x3*(1-x3) s.t. x1 + x2 + x3 = 1 x1 >= 0 and x1 <= 1 x2 >= 0 and x2 <= 1 x3 >= 0 and x3 <= 1 which are the constraints. I'm expecting the answer x1=x2=x3 = 1/3. I tried the "constrOptim" function in R and I'm running into some issues. I first start off

i want to estimate parameters with maximum likelihood method with contraints (contant numbers). for example sum(Ai)=0 and sum(Bi)=0 i have done it without the constraints but i realised that i have to use the contraints. Without constraints(just a part-not complete): skellamreg_LL=function(parameters,z,design) { n=length(z); mu=parameters[1]; H=parameters[2]; Apar=parameters[3:10];

i have a function like 1+sin(a+bx) where -pi/2<=a+bx<=pi/2 i made a progrom using constrOptim() function but it is not giving good result. it depends on the initial value. but when i am doing simulation it is impossible of find the best initial value for every step. also i am not exactly sure how to input the restriction though i have read the help file of the function. here x is a set of

Hi, I'm very interested in using the constrOptim() function currently in the R-devel sources. In particular, I'm trying to fit point process conditional intensity models via maximum likelihood. However, I noticed that the gradient of the objective function must be supplied for all but the Nelder-Mead method. I was wondering why this was because optim() itself does not require a gradient

Hi. I need some advises on how to use R to find pi (i is the index) with the following objective function and constraint: max (sum i)[ f(ai, bi, pi) * g(ci, di, pi) * Di ] s.t. (sum i)[ f(ai, bi, pi) * Di * pi] / (sum i)[ f(ai, bi, pi) * Di ] <= constant f and g are diffentiable. So, I am thinking of optim with method = "BFGS"? But wonder how to include the

Dear all, There appears to be a bug in how constrOptim handles ... arguments that are suppose to be passed to optim, according to the documentation. This means you can't get the hessian to be returned, for example (so this is a real problem, and not just a question of mistaken documentation). Looking at the code, it appears that a call to the user-defined f includes the ..., when the ...