similar to: Likelihood returning inf values to optim(L-BFGS-B) other

I have a function called within .First(), as in .First=function() { ...blah... BIG(partofblah) #BIG is my function, n=partofblah when called ...foo... } The partofblah component of blah is a number obtained from readline(), which is then an argument to BIG() BIG=function(n=10,removeask=T) { z <- sapply(ls(pos=1), function(x)object.size(get(x))) ...stuff... } When .First()

Dear R-help list, I am working on an optimization with R by evaluating a likelihood function that contains lots of Gamma calculations (BGNBD: Hardie Fader Lee 2005 Management Science). Since I am forced to implement lower bounds for the four parameters included in the model, I chose the optim() function mith L-BFGS-B as method. But the likelihood often returns inf-values which L-BFGS-B

Nurza, Try running a while loop steping out until you have a start and finish thats the function is opposite in sign. You need a "start" and "finish" where F is + and - on either side of the loop. Graphing F might help. step<-10 checkme<-F(start)*F(finish+step) while(checkme>0){ initialstep<-initialstep*2 checkme<-F(start)*F(finish+step) }

Hi all: Well, after making the change to make use of the CVS trees and switching to ices2 (as opposed to the libshout example streamer), the icecast server at bumpy.braille.uwo.ca:9000 is back up. I've put my old 22khz mono files (encoded with beta4) back up on test.ogg. In addition, we've implemented a vorbis subset of Kirk Reiser's egoplay MP3 stream. You can now hear all the

Hi all: Well, after making the change to make use of the CVS trees and switching to ices2 (as opposed to the libshout example streamer), the icecast server at bumpy.braille.uwo.ca:9000 is back up. I've put my old 22khz mono files (encoded with beta4) back up on test.ogg. In addition, we've implemented a vorbis subset of Kirk Reiser's egoplay MP3 stream. You can now hear all the

I have been dabbling with getting some of our Windows applications running under Linux using the CodeWeavers preview 20010629. I am making good progress and am very pleased with the results. However I am currently stuck on one problem. The Windows software is developed using Borland Delphi and the application includes a number of forms that are intended to be displayed full screen. The

Hi Devs, I'm developing a tool that relies on semantic information in bytecode labels (i.e. block names). I've discovered that clang is stripping these named labels (along with some virtual register names) when I run on a virtual machine. Well, I'm using VirtualBox, and have tried two different versions of Ubuntu and some different clang builds. Can someone point me in the right

Exception type exceptions.ValueError Exception Handler Information Traceback (most recent call last): File "/usr/share/rhn/rhn_applet/rhn_applet.py", line 521, in refresh self.model.refresh(force) File "/usr/share/rhn/rhn_applet/rhn_applet_model.py", line 457, in refresh source.refresh(force) File "/usr/share/rhn/rhn_applet/rhn_applet_source.py", line

Hello everybody, In order to predict income for different time points, I fitted a linear model with polynomial effects using BMA (bicreg(...)). It works fine, the results are consistent with what we are looking for. Now, we would like to predict income for a future time point t_next and of course draw the prediction interval around the estimated value for this point t_next. I've found the

Dear All, I am running Wine under Linux - Suse 10.1 (kernel 2.6.16). GUI apps work fine (eg "wine notepad", "wine matlab") and wcmd works too ("wine wcmd"), but any CUI program causes a crash, like this: /home/gcunning> wineconsole --backend=user "c:freedos\edit.com" Warning: unprotecting memory to allow real-mode calls. NULL pointer accesses

Dear R users, I am new to R. I would like to find *maximum likelihood estimators for psi and alpha* based on the following *log likelihood function*, c is consumption data comprising 148 entries: fn<-function(c,psi,alpha) { s1<-sum(for(i in 1:n){(c[i]-(psi^(-1/alpha)*(lag(c[i],-1))))^2* (lag(c[i],-1)^((-2)*(alpha+1)) )}); s2<- sum(for(m in 1:n){log(lag(c[m],-1)^(((2)*alpha)+2))});

Hello I have the following problem : The dataframe TEST has multiple lines for a same person because : there are differents values of Nom or differents values of Prenom but the values of Matricule or Sexe or Date.de.naissance are the same. TEST <- structure(list(Matricule = c(66L, 67L, 67L, 68L, 89L, 90L, 90L, 91L, 108L, 108L, 108L), Nom = structure(c(1L, 2L, 2L, 4L, 8L, 5L, 6L, 9L, 3L, 3L,

On Fri, 21 Mar 2003, Remco B. Brink wrote: > <quote who="Geoff Shang"> > > > Is there some XML magic I need to use now to get directory info > > published? > > Can you send over your existing configuration? Yeah. Ignore the comments, this began as the sample config file and the comments weren't removed or changed. <?xml

hey all, maybe i am dumb but it took me a few hours to work out which cvs was which and where ogg went etc etc...so i am completing some documentation (the icantcode site was down too when i checked last time, is it still alive?) anyway, I am going to make some good documentation for my site soon but heres a skeleton including some snippets from docs online already - appologies in advance for

Hi Folks: I have written a simple scheduler for ices2. It basically plays a background playlist in sequential, random or shuffle mode. Each time it is called by ices it checks to see if there is a file in /tmp and if there is builds it into a request file which plays immediately after the current piece finishes. You can use 'at' or cron to set up jobs to place playlists in /tmp to be

Hi Everyone, I'm having a problem with arima.sim. Namely specifying inital values for the series. If I generate a random walk > vs = rnorm(100,0,1) > xs = cumsum(vs) and fit an ARIMA(1,0,0) to it > xarima = arima(xs,order=c(1,0,0)) > xarima Call: arima(x = xs, order = c(1, 0, 0)) Coefficients: ar1 intercept 0.9895 8.6341 s.e. 0.0106 6.1869 I should

I'll add our init script for redhat to it then: note: it does work, but restarting is not always working, due to ports still being open :-( --------------------------------------------------------------------- #!/bin/sh # # description: starts/stops the ices audioencoder # # chkconfig: 345 99 20 # # Source function library. . /etc/rc.d/init.d/functions

Dear, I want to compute coefficient of determination (R-squared) to complement AIC for model selection of multivariable GLM. However, I found this is not a built-in function in glm. neither is it available through reviewing the question in the R-help archive. Please kindly help and thanks a lot. Elaine [[alternative HTML version deleted]]

Hello I'm new to the list. I went through recent archives and didn't find any mention of this, sorry if this has already been discussed. When listening to a stream (http://bumpy.braille.uwo.ca:9000/test.ogg) winamp crashed after a few minutes. When i tried going back in and went back to the stream it worked for a few minutes again and crashed. The times are different (8:27 and 6:20, and

Hey folks, I'd like to get a .deb of the unmodified_drivers put together. I just filed ITP, and I would love to have your help with it. http://wp.colliertech.org/cj/?p=653 I've never built a kernel package before, nor have I used DKMS, so it's going to be a bumpy ride for me. Do any of you have experience with these tools? Have you built the unmodified_drivers modules recently?