similar to: constrOptim (PR#8158)

Greetings, I have a nagging problem with contrasts and I can't seem to resolve it. A factor exists with four levels (lib1, lib2, con1, con2) and when I check the contrasts or set the contrasts to any of the prespecified ones, I do not get the exact contrasts necessary to test the theoretically relevant ones. I need orthogonal contrasts that look just like this matrix: con1 con2

Dear r helpers I run the following code for nested logit and got a message that Error in optim(c(0, 0, 0, 0, 0.1, -2, -0.2), fr, hessian = TRUE, method = "BFGS") : initial value in 'vmmin' is not finite What does this mean? and how can I correct it? Thank you June > yogurt = read.table("yogurtnp.csv", header=F,sep=",")> attach(yogurt)>

Dear Rhelpers I have the problem with initial values, could you please tell me how to solve it? Thank you June > p = summary(maxLik(fr,start=c(0,0,0,1,0,-25,-0.2))) Error in maxRoutine(fn = logLik, grad = grad, hess = hess, start = start, : NA in the initial gradient > p = summary(maxLik(fr,start=c(0,0,0,1,0,-25,-0.2),method="BFGS")) Error in optim(start, func, gr =

I'm doing some string manipulation on a vector of file names, and noticed something curious. When I strsplit the vector, I get a list of character vectors. The list is numbered, as lists are. When I cast that list as a data frame with 'as.data.frame()', the resulting columns have names derived from the original filenames. Example code is below. My question is, where are these names

Hi R users, I am new to R and am trying to merge data frames in the following way. Suppose I have n data frames each with two fields. Field 1 is common among data frames but may have different entries. Field 2 is different. Data frame 1: Src Target1 1 aaa 1 bbb 1 ccc 2 aaa 3 ddd Data frame 2: Src Target2 2 aaaa 3 dddd 4 bbbb 4

My decision tree grows only with one split and based on what I see in E-Miner it should split on more variables. How can I adjust splitting criteria in R? Also is there way to indicate that some variables are binary, like variable Info_G is binary so in the results would be nice to see "2) Info_G=0" instead of "2) Info_G<0.5". Thank you in advance! And thanks for Eric who

Hi when I use con1 <- socketConnection(...) in R and want to send text from another application written in Delphi to R, do I just have to send the text or do I have to implement more control characters and so on? Is con1 <- socketConnection(port=6011, server=TRUE) writeLines("plot(rnorm(100))", con1) just sending the text in "plot(rnorm(100))" to the socket or is

I have an application called dynaterm which I can run successfully under wine but to achieve 110% success I need to pass this application a parameter which must contain double quotes E.G. C:/Program Files/dynaterm8/dynaterm.exe -connect:"con1" So I've tried escaping the quotes and setting the launcher command to wine "C:/Program Files/dynaterm8/dynaterm.exe

Please help me, I use R-1.4.1 and everything works fine. But with R-1.5.1 ist does'nt work: I have this put into a file called "socket.R" while(TRUE){ con1 <- socketConnection(host="localhost",port=6021,server=TRUE); source(con1); close(con1); } So I can use R with: source("socket.R") and R answers me when I send a command to it via socket 6021. But in

Using the package /RMySQL/ I have established a connection ('con1') to a MySQL-Database called 'mc'. This database contains a table called 'mc_2000' which has 200000 observations on several variables. One of these variables is called 'B5' (type 'double'). When I generate a numeric vector in R called 'B6' that has the same length as 'B5': How

Dear All, I am using this query it returns id : id <- dbGetQuery(con1,"SELECT id FROM tenants WHERE name LIKE '%consim%'") But In my case the string "consim" is there in another variable(it is coming from configuration file); > str <- "consim" I am trying to replace the string some this like, but it not working: > id <-

Hi Thanks Gabor for your suggestion. I am posting the code that worked for me. dataframe1 = data.frame(cbind(Src = c(1,1,1,2,3), Target1 = c('aaa','bbb','ccc','aaa','ddd'))); #must be data frame dataframe2 = data.frame(cbind(Src = c(2,3,4,4,4), Target2 = c('aaaa','dddd','bbbb','eeee','ffff'))); dataframe3 =

Hi, Is there any way to get the hessian matrix from the "constrOptim" function without supplying gradient function? Thanks. --------------------------------- Bring words and photos together (easily) with [[alternative HTML version deleted]]

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

Hello, I'm trying to run constrOptim. It returns to me an error about the fact that constraints arguments (ui and ci) are non compatibles: > optout= constrOptim(startparams, f=ImpulseSS, grad=grImpulse, ui=UI, ci=CI, data=gexp[k,], t=t) Error in ui %*% theta : non-conformable arguments I would like to point out that I can calculate that product in the command line: > UI %*%

Hi, I'm having trouble using the constrOptim function to generate the 9-component vector argmin of the function ELfsds: ELfsds <- function(pvechat){ LG=0 for(i in 1:9){ LG=LG+log(pvechat[i]) } return(-LG) } with accompanying gradient function: gradfunc <- function(thetavec){ g=1/(9*thetavec) return(g) } The constraints on the optimization problem are: 1 - components of

I have a function myFunction(beta,x) where beta is a vector of coefficients and x is a data frame (think of it as a matrix). I want to optimize the function myFunction() by ONLY changing beta, i.e. x stays constant, with 4 constraints. I have the following code (with a separate source file for the function): rm(list=ls()) source('mySourceFile')

Hi everybody, please, could you give me help ? I scanned the help archives and didn't found hints... I want to solve a large sparse linear system subjected to an inequality constrains (all solutions positive) and an equality constrain (all solutions sum to 1), thus I tried to fool constrOptim using: x[1] + 0 + ... + 0 >= 0 ... 0 + 0 + ... + x[n] >= 0 x[1] + x[2] + ... +

Hello, I got a question with the R function constrOptim. >From the R help, it says that the return values of "constrOptim" are the same as "optim". For the return value "convergence" of the function "optim", the values should be 0, 1, 10, 51 and 52. See http://www.maths.lth.se/help/R/.R/library/stats/html/optim.html When I use constrOptim, I get

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