search for: grr

...lt;- x[2] ((1/(1+exp(-b0+b1))+(1/(1+exp(-b0)))+(1/(1+exp(-b0-b1)))))/3 } As you can see, my objective function is ((1/(1+exp(-b0+b1))+(1/(1+exp(-b0)))+(1/(1+exp(-b0-b1)))))/3 and I would like to solve for both b0 and b1. If I were to use optim then I would derive the gradient of the function (grr) as follows: fr2 <- expression(((1/(1+exp(-b0+b1))+(1/(1+exp(-b0)))+(1/(1+exp(-b0-b1)))))/3) grr <- deriv(fr2,c("b0","b1"), func=TRUE) and then simply use optim via optim(c(-5.2,0.22), fr1, grr) My problem is that I wish to place constraints (b0>=-0.2 and b1>= 0....

...1+exp(-b0)))+(1/(1+exp(-b0-b1)))))/3 > } > > As you can see, my objective function is > ((1/(1+exp(-b0+b1))+(1/(1+exp(-b0)))+(1/(1+exp(-b0-b1)))))/3 and I would > like to solve for both b0 and b1. > > If I were to use optim then I would derive the gradient of the function > (grr) as follows: > > fr2 <- > expression(((1/(1+exp(-b0+b1))+(1/(1+exp(-b0)))+(1/(1+exp(-b0-b1)))))/3) > grr <- deriv(fr2,c("b0","b1"), func=TRUE) > > and then simply use optim via > > optim(c(-5.2,0.22), fr1, grr) > > My problem is that I wish to...

Dear Friends. I found something very puzzling with constOptim(). When I change the parameters for ConstrOptim, the error messages do not seem to be consistent with each other: > constrOptim(c(0.5,0.3,0.5), f=fit.error, gr=fit.error.grr, ui=ui,ci=ci) Error in constrOptim(c(0.5, 0.3, 0.5), f = fit.error, gr = fit.error.grr, : initial value not feasible > constrOptim(c(0.5,0.9,0.5), f=fit.error, gr=fit.error.grr, ui=ui,ci=ci) Error in constrOptim(c(0.5, 0.9, 0.5), f = fit.error, gr = fit.error.grr, : initial val...

...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 by setting my equalities to inequalities x1+x2+x3 >= 1 x1+x2+x3 <= 1.001 However, I get a "Error in constrOptim(c(1.00004, 0, 0), fr, grr, ui = t(A), ci = b) :   initial value not feasible Execution halted" error. The values 1.00004,0,0 are well within the constraints domain. fr = function(x){   x1 = x[1]   x2 = x[2]   x3 = x[3]   x1*(1-x1) + x2*(1-x2) + x3*(1-x3) } grr = function(x){   x1 = x[1]   x2 = x[2]   x3 = x[3]   c(1...

...ange in the 'Rosenbrock Banana function' (taken from the Examples section of those help pages) adding extra (useless) parameters: fr <- function(x,extra1,extra2) { ## Rosenbrock Banana function x1 <- x[1] x2 <- x[2] extra1*extra2*100 * (x2 - x1 * x1)^2 + (1 - x1)^2 } grr <- function(x,extra1,extra2) { ## Gradient of 'fr' x1 <- x[1] x2 <- x[2] extra1*extra2*c(-400 * x1 * (x2 - x1 * x1) - 2 * (1 - x1), 200 * (x2 - x1 * x1)) } now, call it: optim(c(-1.2,1), fr, grr, extra1=c(1,2), extra2=c(3,4)) ok! but now: constrOptim(c...

...know another way to solve my problem of maximization under constraints? THANKS! My R instructions n=c(10,20,3,5) n1=n[1] n2=n [2] n3=n[3] n4=n[4] logfr=function(x) { ##function to maximize x1= x [1] x2= x[2] x3= x[3] x4= x[4] log(factorial(sum(n)))-sum(log (factorial(n)))+sum(n*log(x)) } grr.log <- function(x) { ## Gradient of 'log fr' x1=x[1] x2=x[2] x3=x[3] x4=x[4] return(n/x) } par.start= c(.19999999,.15,.4,.25) constr.coeff = rbind(diag(1,4,4),c(-1,1,1,1),c (0,-1,1,1),c(-1,-1,-1,-1), c(1,1,1,1)) constr.tn= c(0,0,0,0,0,0,-1,. 9999999) min= constrOptim(par.start, logf...

Hmm... i think it's time to go home now... forgot attachment. //T-? >Hi all. > >I am submitting a patch that adds DCE authentication and credential > >attaching. Fully legit, no ugly hacks ;) > >I looked (some) at Paul Henson's <henson at acm.org> patch for ssh-1.2.20 >while coding. > >The patch [attached] is against pre14, but DCE and AIX patches

...okay if (obj > obj.old) break } ###Here is an example modified from the examples of the help page of constrOptim: > fr <- function(x) { ## Rosenbrock Banana function + x1 <- x[1] + x2 <- x[2] + 100 * (x2 - x1 * x1)^2 + (1 - x1)^2 + } > grr <- function(x) { ## Gradient of 'fr' + x1 <- x[1] + x2 <- x[2] + c(-400 * x1 * (x2 - x1 * x1) - 2 * (1 - x1), + 200 * (x2 - x1 * x1)) + } > constrOptim(c(-1.2,0.9), fr, grr, ui=rbind(c(-1,0),c(0,-1)), ci=c(-1,-1)) $par [1] 0.9999761 0.9999522 $value [1]...

First, Thanks for any and all help!!!! I can't seem to figure out what I need to do, I've been fighting this for a month and am now beating my head off my desk with no solution to be found. I've read others having this issue but they were all older versions.. I am using 3.5.4,, Please read over and give me some input.. Every 7 days winbindd fails on the trust secret. The only way

Dear friends. I am using function constrOptim(c(0.5,0.3,0.5), fit.error, fit.error.grr, ui=-1*ui,ci=-1*ci) and I am confronted with error message "initial value not feasible" I plug in the initial value of (0.5,0.3,0.5) to function fit.error and fit.error.grr and have pretty reasonable result. I inequality "ui %*% theta - ci >= 0" as suggested in the R manual...

Hi, i have a function: zz<- (constrOptim(c(.5,0), fr, grr, ui=rbind(c(-1,0),c(1,-1)), ci=c(-0.9,0.1))) i can get the result by using command (for example): zz$par now if i can use multicore: zz<-parallel(constrOptim(c(.5,0), fr, grr, ui=rbind(c(-1,0),c(1,-1)), ci=c(-0.9,0.1))) result < collect(zz) i cant get my the result: result$par because mult...

Hi I am trying to install "rms" package but while installing it shows following error package 'survival' 2.36-2 is loaded, but >= 2.36.3 is required by 'rms' what to do? i am using linux OS I have tried by updated r-base-core but it didnt work regards GRR [[alternative HTML version deleted]]

Apparently nobody has a suggestion to why smbd process all of a sudden starts sucking CPU like crazy? Restarted samba and then restarting the entire server does not fix it. Searching the archives shows I'm not the only one with this problem. The archives do not show what to do about it yet. Grr... Redhat 9.0 Samba 3.04 ( I have 2.2.7 servers that have done this to me also) - Gary

...r*(tot/sh*1000.0)))/(sigmae*ss/(ss+lbar*(tot/sh*1000.0 ))))*(sigmae*ss/(ss+lbar*(tot/sh*1000.0 )))*sqrt((lambda*lambda/(sigmae*ss/(ss+lbar*(tot/sh*1000.0 )))/(sigmae*ss/(ss+lbar*(tot/sh*1000.0))))))))))) fit.error=function(rec,lambda,lbar) {sum((eval(apple)*1000-orange)^2/(orange^2)) } fit.error.grr=function(rec,lambda, lbar) {drec=sum(eval(D(apple,'rec'))*(eval(apple)*1000-orange)/(orange^2)) dlambda=sum(eval(D(apple,'lambda'))*(eval(apple)*1000-orange)/(orange^2)) dlbar=sum(eval(D(apple,'lbar'))*(eval(apple)*1000-orange)/(orange^2)) c(drec,dlambda,dlbar) } rr=c(4.33...

...39;wx/keyword_ctors'' require ''wx/keyword_defs'' These files do appear in the svn sources. If I copy the files from the svn sources, I can now run the examples. Unfortunately, I''ve done too much mucking around to know exactly where the segfault issue went away. Grr. _______________________________________________ wxruby-users mailing list wxruby-users@rubyforge.org http://rubyforge.org/mailman/listinfo/wxruby-users

.../log/samba/log.%m max log size = 1000 dns proxy = No interfaces = eth1 lo bind interfaces only = yes log level = 2 msdfs:8 auth:5 winbind:5 idmap:5 acls:3 panic action = /usr/share/samba/panic-action %d valid users = @"GREENREEFERS\grr" vfs objects = acl_xattr map acl inherit = Yes store dos attributes = Yes Check that /etc/krb5.conf exists and looks like this: [libdefaults] default_realm = GREENREEFERS.NO dns_lookup_realm = false dns_lookup_kdc = true Check that /etc/resolv.co...

Hi, I am new to R, and am working on some optimization problems - I was wondering if there was a way that I could show all the iterations in R -i.e. showing for each iteration, what the iteration is, how much the function is, what the norm of the gradient would be... Any help is greatly appreciated! thank you!! -- View this message in context:

Hi All. As I mentioned earlier, I've been having weird failures with both sshd and sftp. The sshd one turned out the be my fault (misplaced "}", grr) but the sftp one doesn't appear to be. The sftp regression tests fail on the current portable tree on Solaris and AIX (with my own mods to auth.c and regress/, but I don't think that's related). The test that fails is: $ ./sftp -P sftp-server -B 5 -R 1 -b ./batch Attaching to sftp-...

...charge me a monthly fee (I don't send enough faxes to justify it). T.38 is a requirement; I need to know that a fax has gone through at the time I send it (store-and-forward, efax, etc aren't good enough). I don't need a DID. Broadvox offers T.38, but don't offer pay-as-you-go. Grr.... - a -- "I didn't see it then, but it turned out that getting fired was the best thing that could have ever happened to me. The heaviness of being successful was replaced by the lightness of being a beginner again, less sure about everything. It freed me to enter one of the mos...

This time, it seems that browser() wins. Grr ... -k =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-devel-request@...