similar to: vectorizing sapply() code (Modified by Aaron J. Mackey)

I could use some help generating a time series for the Mackey-Glass equation: dx/dt = 0.2 x(t-tau)/(1 + x(t-tau)^10) - 0.1 x(t) for the case there tau = 17. I tried the "ddesolve" package but the dde(...) function seems to hang, not producing anything. Can someone show me the R script how to do this? - Mike Beddo [[alternative HTML version deleted]]

I have a MLE task that for a small number of parameters finishes in a reasonable amount of time, but for my "real" case (with 17 parameters to be estimated) either takes far too long (over a day), or fails with "computationally singular" errors. So a) are there any parallel implementations of optim() (in R or otherwise) and b) how can I make my function more robust?

http://tinyurl.com/36pw6k Aaron swartz attacs Reddit brings it down --~--~---------~--~----~------------~-------~--~----~ You received this message because you are subscribed to the Google Groups "Ruby on Rails: Talk" group. To post to this group, send email to rubyonrails-talk-/JYPxA39Uh5TLH3MbocFFw@public.gmane.org To unsubscribe from this group, send email to

Hi all, I'm a PhD student at sydney uni and am trying to run a non linear mixed model program to obtain estimates of parameters describing dairy cow lactation curves. At present, I have been able to get the data to converge using the S plus (S plus 2000) nlme function. However, when I put the same data into R (R 1.9.0), add in the nlme package and run the code, it does not converge by the

Hi I'm sorry about the triviality of my problem. I have a vector (v) of three columns (logA, logB, id). I want to compute (and plot) the correlation between logA and logB for different thresholds of id (e.g. >30, etc). So I tried: for(i in 1:100){ points(cor(v$logA[v$id>i], v$logB[v$id>i], use="complete.obs"), i)) } (i created a plot object already) but it comes with

In findInterval's help page, it says 'v[i[j]] <= x[j] < v[i[j] + 1]'. I'm wondering if there is a variant of it such that 'v[i[j]] < x[j] <= v[i[j] + 1]'.

Dear all, Suppose I have a linear mixed-effects model (from the package nlme) with nested random effects (see below); how would I present the results from the random effects part in a publication? Specifically, I?d like to know: (1) What is the total variance of the random effects at each level? (2) How can I test the significance of the variance components? (3) Is there something like an

For the life of me I couldn't work out what to searc I have an m*n numeric matrix x and a numeric vector y (of length m+n-1) How do I do a calculation like this? z[i,j] = x[i,j] * y[i+j] ? Well, one can write a pair of loops, or write a single loop within which we calculate a vector at a time, but ... is there a "neat" way to do it? tiny example:

I have not devoted time to setting up ROracle since binaries are not available and it seems to require some effort to compile (see http://cran.r-project.org/web/packages/ROracle/index.html). On the other hand, RODBC worked more or less magically once I set up the data sources. What is your success using ROracle and why would it be preferable to RODBC ? -Avram On Thursday, September 11,

Hi all, given this example #start a<-c(0,70,50,100,70,650,1300,6900,1780,4930,1120,700,190,940, 760,100,300,36270,5610,249680,1760,4040,164890,17230,75140,1870,22380,5890,2430) length(a) b<-c(0,0,10,30,50,440,1000,140,70,90,60,60,20,90,180,30,90, 3220,490,20790,290,740,5350,940,3910,0,640,850,260) length(b) out<-ks.test(log10(a+1),log10(b+1)) # max distance D

Is there an R-language interface to the R API C-language functions Rf_logspace_add() and Rf_logspace_sub()? I don't see one but I may not looking under the right name. Various packages have functions which do that same sort of thing (log(exp(x)+exp(y)) and log(exp(x)-exp(y)) without unnecessary floating point errors). They have names like matrixStats::logSumExp(lx, na.rm=FALSE, ...)

[sent this to ext2-devel but haven't seen any other traffic there, so cc-ing it here] [Last week] while configuring a new box I came across a site suggesting: chattr -R -S +j +A /var/spool/postifx to improve postfix's speed on ext3, in lieu of mounting a separate fs there w/ data=journal. So I tried it. Both on the new box (2.4.21-7um) and on my laptop (about three csets shy of the

Just now I ran tune2fs -j on the root filesystem of a box running 2.6.0-test8. Then I edited /etc/fstab and changed the FS type to from ext2 to ext3, saved the file, and invoked vim on the file again. A few moments after this, the box hung. Unfortunately X was running at the time, and so I don't have any messages to cite. Is this a known problem?

All -- I'm curious is if it safe or even wise to run the following command on a mounted filesystem, namely root (/)? tune2fs -j /dev/hda1 What about if someone goes into single user mode and runs this first? mount -o remount,ro / And then to enable it, runs this? mount -t ext2 -o remount,rw / I assume it is not safe to do so, but one user in my LUG assumes otherwise. Just curious,

Hi list, I just noticed in the chattr(1) man page that there exists a 'j' attribute which actually does the same thing that data=journal mount option would do, expect it only sets this property on the directories of my choice. Big Thanks to the one that came up with that idea :) It allows me to set +j on mail queues that must sit on the same partitions as the rest of the system (the guy

Hi All I have Googled the message Error in J[time] : invalid subscript type 'closure', but can't quite make sense from the results that I get. What does it mean? I'm trying to access an element in an array. Here is a short bit of code that produces the error: s_fun <- function(parms,time) { with(as.list(c(parms)),{ print(J[time]) }) } times = seq(from=0, to=100, by =

Hi: I have a query related to the J and Jcross functions in the SpatStat package. I use J to finding indications of clustering in my data, and Jcross to look for dependence between point patterns. I use the envelope function to do Monte Carlo tests to look for significance. So far so good. My question is how I can test to see if tests are significantly different. For example, if find J of

So, I'm building 3.0.23b for Slackware (since they don't have it out yet[1]), and I've noticed these two lines in the script that Slackware uses to build Samba 3.0.23 from source (which I'm modifying to build 3.0.23b): # -j options don't seem to work... make Anyone know why that comment might be there? Is it true that Samba can't be built with "-j2" or

I'm sorry if this has already been posted - I've only been subscribed for a month or two and I did search thoroughly through the archives. Basically I am using Syslinux to boot my flash disk with Linux (MOVIX2). I am using a Kingston Data Traveler 2.0. When syslinux -s /dev/sda is executed I get an error of 'ERROR 3246: Boot sector read failed' yet it seems that something has

https://bugzilla.netfilter.org/cgi-bin/bugzilla/show_bug.cgi?id=78 Summary: -m psd -j TARPIT returns all ports open from nmap Product: iptables userspace Version: unspecified Platform: i386 OS/Version: RedHat Linux Status: NEW Severity: normal Priority: P2 Component: unknown AssignedTo: