similar to: Parameters setting in functions optimization

Dear R users I am trying to use OPTIMX(OPTIM) for nonlinear optimization. There is no error in my code but the results are so weird (see below). When I ran via OPTIM, the results are that Initial values are that theta0 = 0.6 1.6 0.6 1.6 0.7. (In fact true vales are 0.5,1.0,0.8,1.2, 0.6.) -------------------------------------------------------------------------------------------- >

Your function is giving NaN's during the optimization. The R-forge version of optimx() has functionality specifically intended to deal with this. NOTE: the CRAN version does not, and the R-forge version still has some glitches! However, I easily ran the code you supplied by changing optim to optimx in the penultimate line. Here's the final output. KKT condition testing Number of

Dear R users When I use OPTIMX with BFGS, I've got the following error message. ----------------------------------------------------------------- > optimx(par=theta0, fn=obj.fy, gr=gr.fy, method="BFGS") Error: Gradient function might be wrong - check it! ----------------------------------------------------------------- So, I checked and checked my gradient function line by

Hello! I am trying to create an R optimization routine for a task that's currently being done using Excel (lots of tables, formulas, and Solver). However, otpim seems to be finding a local minimum. Example data, functions, and comparison with the solution found in Excel are below. I am not experienced in optimizations so thanks a lot for your advice! Dimitri ### 2 Inputs:

I am still getting these in my maillog: Jul 13 07:39:38 pop5.cfu.net dovecot: IMAP(breu at cfu.net): Corrupted index cache file /var/spool/mail/filer/storage//cfu.net/b/r/breu/dovecot.index.cache: invalid record size Jul 13 09:35:03 pop5.cfu.net dovecot: IMAP(breu at cfu.net): Corrupted index cache file /var/spool/mail/filer/storage//cfu.net/b/r/breu/dovecot.index.cache: invalid record size

Hi Please let me know if we can store the linear model object in the data base and retrive the object and output from them Data<- read.csv("C:/FE and RE.csv") Formula="Y~X2+X3+X4 lmobject = lm(formula=Formula,data=Data) can i store the lm object in the database and and is it possible to retrive it and get the summary information -- View this message in context:

Dear R-helpers, I need to show a linear fit through a subset of the data within each combination of levels of two factors. So I prepared an xyplot with different panels for each level of one of the factors, and different symbols within each panel for the levels of the second factor. My problem is selecting the subset of each combination through which the line should be fit for subsequent

Hello, I'm trying to build a frequency table for a vector, broken down by the combination of factors. For further analyses, I need to have the result arranged in a new data frame, with the upper limit of the histogram's breaks, the per bin count, and the factors identifying each record. My problem is including the latter: ---<---------------cut

Hello: The development version of Ecdat on R-Forge contains a vignette in which optim(?, method=?L-BFGS-B?) stops with an error message while violating the lower bound. To see all the details, try the following: install.packages("Ecdat", repos="http://R-Forge.R-project.org") Then do "help(pac=Ecdat)" -> "User guides, package

Hi people, I have a text file like this one posted: snp_id gene chromosome distance_from_gene_center position pop1 pop2 pop3 pop4 pop5 pop6 pop7 rs2129081 RAPT2 3 -129993 "upstream" 0.439009 1.169210 NA 0.233020 0.093042 NA -0.902596 rs1202698 RAPT2 3 -128695 "upstream" NA

Hello, I am attempting to utilize the 'alabama' package to solve a constrained nonlinear optimization problem. The problem has both equality and inequality constraints (heq and hin functions are used). All constraints are smooth, i.e. I can differentiate easily to produce heq.jac and hin.jac functions. My initial solution is feasible; I am attempting to maximize a function, phi. As

http://en.wikipedia.org/wiki/Binomial_distribution Bino Hi there, How can I get the legend in the probability density plot in that http://en.wikipedia.org/wiki/File:Binomial_distribution_pmf.svg wiki page ? and the cumalative plot as well? ==================================== # Binomial Distribution x=0:40

Still trying to get a handle on c, cat, paste, and lists. Here is my current problem. I have a good number of vectors that I want to selectively combine into a dataframe. The exact name of the vectors changes in a looping arrangement but can easily be grepped. The problem is: how can I use a list produced by a grep in a data.frame command: i.e. > ls() [1] "a"

Hi there, is there a possibility to count the number of appearances of an element in a vector ? i mean of any given element.. deliver all elements which are exactly xtimes in this vector ? thx in advance !!

This message is in MIME format --_=XFMail.1.1.p0.Linux:970825095458:252=_ Content-Type: text/plain; charset=us-ascii Here is a patch which fixes some problems with time series functions. Some examples of what goes wrong... > x <- ts(rnorm(100),start=1,deltat=2) > start(x) Error in ts[1] : object is not subsettable > end(x) Error in ts[2] : object is not subsettable > y <-

I'm using Dovceot/IMAP on Linux and Outlook clients on WIN7 workstations. Mail on Linux is stored in Maildir format. I'm searching for where Outook keeps its information on color categories in IMAP. According to Diane Poremsky at slipstick.com, "Outlook stores it in the mapi properties of each message. If you use MFCMAPI to viuw the messages, you'll see the properties."

Hi, Is it possible to create a ts() object, whose data is daily based BUT measured only on working days? In other words, suppose I have a data set with 255 observations, measured from 29 June 1959 to 30 June 1960. How would I create such a data? I tried something like: ts(c(...), start(1959, 180)) but I'm not sure what to use for frequency. In other words I don't know how to

Running samba version 3.2 on AIX 6.1. Mapping share on windows 7 workstation. I can map the share using \\<ip<file:///\\%3cip> address>\testarea. If I attempt to map using using \\servername\testarea<file:///\\servername\testarea> I receive an error. smbclient -L //<servername>/testarea on the unix box is successful When I run a net view command I receive the

Hi, I'm using Cepstral as a TTS Engine for Asterisk with Swift application. It works fine when I have just 1 voice installed. Now I have 2 voices in the same language installed but I can't seem to find the way to select which voice to use in Swift's application in Asterisk. Does anyone know?? Thank you, Pim

I've run across a situation where I have a directory of .gz files. I do want to rsync this directory. I've tried without the --compress (because they are already compressed) but the files never get rsynced. Interactively, I get the following error: :rsync error: some files could not be transferred (code 23) at main.c(1045) In the /etc/rsyncd.log, for each of the files, I get: