similar to: Optimizing nested function with nlminb()

On Tue, Jan 10, 2017 at 08:13:21PM -0600, Dirk Eddelbuettel wrote: > > On 10 January 2017 at 17:48, frederik at ofb.net wrote: > | Hi R Devel, > | > | I just ran into a corner case with 'strptime'. Recall that the "%OS" > | conversion accepts fractional seconds: > | > | > strptime("17_35_14.01234.mp3","%H_%M_%OS.mp3")$sec > |

Hi R Devel, I just ran into a corner case with 'strptime'. Recall that the "%OS" conversion accepts fractional seconds: > strptime("17_35_14.01234.mp3","%H_%M_%OS.mp3")$sec [1] 14.01234 Unfortunately for my application it seems to be "greedy", in that it tries to parse a decimal point which might belong to the rest of the format: >

1. lbeta and beta do not work properly: lbeta returns its first argument and beta gives the lbeta result. In names.c lines 245-6, the codes for these should be 2 and 3 instead of 1 and 2 2. crtl-C does not work (except the first time) on Red Hat elf Linux (which has many many other problems as well) nor on the previous version of Linux for Amiga. It worked on Slackware aout Linux and now works

1. lbeta and beta do not work properly: lbeta returns its first argument and beta gives the lbeta result. In names.c lines 245-6, the codes for these should be 2 and 3 instead of 1 and 2 2. crtl-C does not work (except the first time) on Red Hat elf Linux (which has many many other problems as well) nor on the previous version of Linux for Amiga. It worked on Slackware aout Linux and now works

Hi, How to solve this problem. The following is the code. betabinexch0=function(theta,data) + { + eta=theta[,1] + K=theta[,2] + y=data[,1]; n=data[,2] + N=length(y) + val=0*K; + for (i in 1:N) + val=val+lbeta(K*eta+y[i],K*(1-eta)+n[i]-y[i]) + val=val-N*lbeta(K*eta,K*(1-eta)) + val=val-2*log(1+K)-log(eta)-log(1-eta) + return(val) + } > data(cancermortality) >

Dear community, I am currently trying to fit an ordinal logistic regression model with the polr function. I often get the same error message : "attempt to find suitable starting values failed", for example with : require(MASS) data(iris) polr(Species~Sepal.Length+Sepal.Width+Petal.Length+Petal.Width,iris) (I know the response variable Species should be nominal but I do as levels

Dear friends i am trying to fit an Ornstein-Uhlenbeck process by MAXIMUM LIKELYHOOD method. i found these formulas on http://www.sitmo.com/article/calibrating-the-ornstein-uhlenbeck-model/ this is the mean-reverting process http://r.789695.n4.nabble.com/file/n4637271/process.txt process.txt and this is the script that i am using....... ouFit.ML=function(spread) { n=length(spread)

What is the likelyhood that major inexpensive webhosts like godaddy, bluehost, etc. will upgrade to RoR 1.1? Is this going to be like PHP 5 where it has to percolate for a year or more before it becomes widly available? Your thoughts? Along the same lines... is it possible to adopt some of the new improved Ajax / javascript capabilities without actually upgrading the ruby installation?

(1) ...need to speed up a monte-carlo sampling...any suggestions about how I can get R to use all 8 cores of a mac pro would be most useful and very appreciated... (2) spent the last few hours trying to get pnmath to compile under os- x 10.5.4... using gcc version 4.2.1 (Apple Inc. build 5553) as downloaded from CRAN, xcode 3.0... ...xcode 3.1 installed over top of above after

NAME xx - twice as dirty SYNOPSIS ~ > gem install "double x" require "xx" include XX::XHTML doc = xhtml_{ html_{ head_{ title_{ " go xx! " } } body_{ " one more and it would be illegal " } } } URI http://rubyforge.org/frs/?group_id=1024 DESCRIPTION xx is a library designed to extend ruby objects

Hi everyone, I am wondering what the likelyhood would be to get a "delete machine script" implemented in the samba core that would delete a machine upon successful un-join of a samba domain. Alternatively, why is this not a good idea? tom.

A friend runs 6.8 but would like to add postscreen to his email server. It is provided with postfix version 2.8 and up. The newest postfix for 6.8 is 2.6.x. Any likelyhood it will be further updated? Or is there another source for newer postfix versions for CentOS 6.8? Similar query for the clamav db, now at 0.99.1 with 0.99.2 out for quite some time. Jon -- Jon H. LaBadie

Hi, I tried to sync two files (50 MB and 100 MB) with my webdav folder using rsync 3.0.5 with Mac OS X (10.4.11) Terminal and X11. With the Terminal application the sync always failed with both files. With X11 I once was successful synchronizing the 50 MB file. The error message is always like this (also when I set the timeout to eg. 7200 seconds): io timeout after 1003 seconds --

Dear all, I'm trying to estimate the parameters of a special case of a poisson model, where the specified equation has an integral and several fixed parameters. I think that the MLE command in STATS4 package could be a good choice, but it's a little complicated. I've got some problems with the offset and I don't understand some of the functions. Do you know where can I find some

Begging indulgence of the list - asking here because due to the higher likelyhood of getting an answer... I have an application that often loses connection from one machine to another when one (or more) of the machines has particular brands of 1gig nics, but which runs rock solid when on 10meg nics and some 1 gig nics. The application senses (falsely) that the connection has been shut down

As some of you R-devel readers may know, the plot() method for "lm" objects is based in large parts on contributions by John Maindonald, subsequently "massaged" by me and other R-core members. In the statistics litterature on applied regression, people have had diverse oppinions on what (and how many!) plots should be used for goodness-of-fit / residual diagnostics, and to my

I am testing the likelyhood of updating a couple RH9 systems to centos 5 386. I found at document that talked about updating to centos 4. steps being: 1.) backup 2.) rpm --import http://pubmirrors.reflected.net/centos/5/os/i386/RPM-GPG-KEY-centOS-5 rpm -Uvh --nodeps http://pubmirrors.reflected.net/centos/5/os/i386/CentOS/yum-3.0.5-1.el5.centos.5.noarch.rpm

Works for me: > strptime("17_35_14.01234.mp3","%H_%M_%OS")$sec [1] 14.01234 > strptime("17_35_14.mp3","%H_%M_%OS")$sec [1] 14 Just leave off the ".mp3" in your time pattern. Relevant section from the help ("Details") for strptime: strptime converts character vectors to class "POSIXlt": its input x is first converted by

[I send this to R-help, since the bug may affect users..] The beta(.) and ln(beta(..)) functions are not giving proper results in R-0.49. This does NOT affect usual beta-distribution calculations, since the internal beta(.) and lbeta(.) are okay. beta(1,2) [1] -0.6931472 exp(beta(1,2)) [1] 0.5 lbeta(3,4) [1] 3 ------- For R-devel'ers : Using gdb clearly shows that beta(.) really

Linking my VC++ application with the standalone Rmath library yields the following; ------ Build started: Project: Complex plugin, Configuration: Debug Win32 ------ Linking... Creating library .\../Debug/complex_plugin.lib and object .\../Debug/complex_plugin.exp libRmath.a(mlutils.o) : warning LNK4217: locally defined symbol __iob imported in function _REprintf libRmath.a(dbeta.o) :