[R] strange timings in convolve(x,y,type="open")

2007-12-18 Thread Art Owen
like N or N log(N) , with the difference hard to discern in granularity and noise. The convolve function is not very big (see below). When type is not specified, it defaults to "circular". So my guess is that something mysterious might be happening inside the first else clause below,

Re: [R] Observations on SVD linpack errors, and a workaround

2007-10-17 Thread Art Owen
> Ph: (410) 502-2619 > > Fax: (410) 614-9625 > > Email: [EMAIL PROTECTED] > > Webpage: http://www.jhsph.edu/agingandhealth/People/Faculty/Varadhan.html > > > > -------- > > > > -Original Message- > From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On > Behalf Of Art

[R] Observations on SVD linpack errors, and a workaround

2007-10-16 Thread Art Owen
failed. A simplistic wrapper is listed below.If SVD failures get very common for lots of people then a better solution would be to have the svd function itself try both ways. Another option is to have the svd code try the Golub and Reinsch algorithm (or some other SVD) on those cases where the