similar to: R-alpha: R 0.50.a2 shortest segfault ...

Here's a trivial patch (Martin!). -k --- src/library/base/man/Arithmetic.orig Tue Aug 5 22:33:16 1997 +++ src/library/base/man/Arithmetic Tue Aug 5 22:35:51 1997 @@ -15,6 +15,7 @@ ALIAS(^) ALIAS(%%) ALIAS(%/%) +ALIAS(Arithmetic) VALUE( These binary operators perform arithmetic on vector objects. They return numeric vectors containing the result of the element

Dear List, I am using constrOptim to solve the following fr1 <- function(x) { b0 <- x[1] b1 <- 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

Dear List, I have the code below, where I am using the constrained optimisation package, 'constrOptim.nl' to find the values of two values, b0 and b1. I have no problems when I enter further variable information DIRECTLY into the functions, fn, and heq. In this instance I require fn to have -0.0075 appended to it, and in the case of heq, h[1] has -0.2. library(alabama)

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),

Hi, I think that there something wrong with the 'constrOptim' max/minimization function because she doesn?t send extra arguments to 'optim' call. Fact: When I use optim in a f(x,theta)-like function, everything goes ok. But using constrOptim with the same function leads to error... Proof: Make a small change in the 'Rosenbrock Banana function' (taken from the Examples

Dear R users, I?m a graduate students and in my master thesis I must obtain the values of the parameters x_i which maximize this Multinomial log?likelihood function log(n!)-sum_{i=1]^4 log(n_i!)+sum_ {i=1}^4 n_i log(x_i) under the following constraints: a) sum_i x_i=1, x_i>=0, b) x_1<=x_2+x_3+x_4 c)x_2<=x_3+x_4 I have been using the ?ConstrOptim? R-function with the instructions

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

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 ...

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 and it is

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 multicore add process id. for example:

Hi, Attached is a patch to 1.2.1pre18 that includes a number of portability changes, and more a more complete HPUX10.20 implementation. Thanks to Ben Taylor's utmpx patch, utmpx and wtmpx support are now in for HP, with a few caveats. It compiles cleanly on HPUX10.20 with gcc, cleanly on Linux, and with a few minor warnings on Solaris 2.6. Most compat changes are made via autoconf which

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

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]]

Dear friends. I use ConstrOptim( ) and got error message "initial value not available". My understanding of "initial value not available" is that one of the following 3 cases happens: 1.The objective function is not well defined at the point of the initial value. 2. The differentiation of the objective function is not well defined at the point of the initial value. 3. The

Sorry to be such a nuisance but after asking for help with compiling R-1.1.0 on an Intel-based machine I now have another problem on the Alpha DS-10 server which my students use. Compilation and instalation proceed perfectly, but when I invoke R I get: [etdx01:bin]$ ./R Fatal error: The X11 shared library could not be loaded. The error was /home/etptupaf/R/lib/R/bin/R_X11.so: undefined

I was going through the example below which is taken from the example section in the R documentation for accessing the function call stack. I am confused and I have 3 questions that I was hoping someone could answer. 1) why is y equal to zero even though the call was done with gg(3) 2) what does parents are 0,1,2,0,4,5,6,7 mean ? I understand what a parent frame is but how do the #'s relate

Not sure if I already reported this ... it seems that as of 0.50.a2, the adj graphics parameter is completely ignored. Try e.g. plot(1:10) title(main = "This is a title", adj = 0) title(main = "This is a title", adj = 1) No difference ... -k =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=- r-devel mailing list -- Read

A set of patches for R-0.50-a1 is now available as ftp://stat.auckland.ac.nz/pub/R/R-0.50-a1.patch1.gz The patches mainly fix problems reported since R-0.50-a1 but some older problems are also fixed. Here is the list of changes. Ross o Many subsetting and mutation problems with the new "expression" type have now been fixed. o When ask=T is set in par() the user is instructed

The files: wx/accessors.rb wx/keyword_ctors.rb wx/keyword_defs.rb do not exist in the tarball: % pwd /home/rwa/Distributions/test/wxruby-1.9.2 % find . -name accessors.rb % but are required in wx.rb: % tail wx.rb Dir.glob(class_files) do | class_file | require ''wx/classes/'' + class_file[/\w+\.rb$/] end # Load in syntax sweetner require ''wx/accessors''

Hi I want to extract columns from a data frame using a vector with the desired column names. This short example uses the select argument in the subset function to accomplish what I am trying to do. Is there a better solution? #names of desired columns colnames <- c("col1","col3") #my data data <-