Re: [Numpy-discussion] Correlation filter

2009-11-20 Thread josef . pktd
On Fri, Nov 20, 2009 at 2:28 PM, Keith Goodman wrote: > On Fri, Nov 20, 2009 at 11:17 AM,   wrote: >> On Fri, Nov 20, 2009 at 1:51 PM, Keith Goodman wrote: >>> def corr3(x, y): >>>    x = x - x.mean() >>>    x /= x.std() >>>    nx = x.size >>>    one = np.ones(nx) >>>    xy = lfilter(x, 1, y) >>>

Re: [Numpy-discussion] Correlation filter

2009-11-20 Thread Keith Goodman
On Fri, Nov 20, 2009 at 11:17 AM, wrote: > On Fri, Nov 20, 2009 at 1:51 PM, Keith Goodman wrote: >> def corr3(x, y): >>    x = x - x.mean() >>    x /= x.std() >>    nx = x.size >>    one = np.ones(nx) >>    xy = lfilter(x, 1, y) >>    sy = lfilter(one, 1, y) >>    sy2 = lfilter(one, 1, y*y) >>  

Re: [Numpy-discussion] Correlation filter

2009-11-20 Thread josef . pktd
On Fri, Nov 20, 2009 at 1:51 PM, Keith Goodman wrote: > On Fri, Nov 20, 2009 at 8:53 AM,   wrote: >> scipy.signal.correlate  would be fast, but it will not be easy to >> subtract the correct moving mean. Subtracting a standard moving mean >> would subtract different values for each observation in

Re: [Numpy-discussion] Correlation filter

2009-11-20 Thread Keith Goodman
On Fri, Nov 20, 2009 at 8:53 AM, wrote: > scipy.signal.correlate  would be fast, but it will not be easy to > subtract the correct moving mean. Subtracting a standard moving mean > would subtract different values for each observation in the window. > > One possibility would be to look at a moving

Re: [Numpy-discussion] Fitting a curve on a log-normal distributed data

2009-11-20 Thread Gökhan Sever
On Thu, Nov 19, 2009 at 9:12 PM, Ian Mallett wrote: > Hello, > > My analysis shows that the exponential regression gives the best result > (r^2=87%)--power regression gives worse results (r^2=77%). Untransformed > data gives r^2=76%. > > I don't think you want lognorm. If I'm not mistaken, that

Re: [Numpy-discussion] Correlation filter

2009-11-20 Thread josef . pktd
On Fri, Nov 20, 2009 at 10:51 AM, Keith Goodman wrote: > I have a short 1d array x and a large 2d array y. I'd like to locate > the places in the y array that are most like (correlated to) the x > array. > > My first attempt, corr1, is too slow. My second attempt, corr2, is > faster but still slow

[Numpy-discussion] Correlation filter

2009-11-20 Thread Keith Goodman
I have a short 1d array x and a large 2d array y. I'd like to locate the places in the y array that are most like (correlated to) the x array. My first attempt, corr1, is too slow. My second attempt, corr2, is faster but still slow. I reuse the same y many times, so my third attempt will probably

Re: [Numpy-discussion] memmap limits

2009-11-20 Thread Francesc Alted
A Thursday 19 November 2009 00:48:13 Robert Kern escrigué: > On Wed, Nov 18, 2009 at 17:43, Mathew Yeates wrote: > > What limits are there on file size when using memmap? > > With a modern filesystem, usually you are only limited to the amount > of contiguous free space in your process's current

Re: [Numpy-discussion] Failure building pdf doc

2009-11-20 Thread Pauli Virtanen
Fri, 20 Nov 2009 10:55:35 +0900, David Cournapeau wrote: > While checking everything builds for the 1.4.0 release, I noticed a > problem with building the latex version: > > writing... done > processing numpy-user.tex... user/index user/introduction > user/whatisnumpy user/install user/howtofind u

Re: [Numpy-discussion] datetime64

2009-11-20 Thread Fernando Perez
Hi Travis, [...] On Thu, Nov 19, 2009 at 10:03 PM, Travis Oliphant wrote: > Again, the NEP has not been fully implemented yet.   What is implemented > works as far as I can tell, but could use more tests.   I would like to > finish the core functionality before 1.4.0 and will try to do that, but

Re: [Numpy-discussion] Numpy/Scipy for EC2

2009-11-20 Thread David Goldsmith
On Thu, Nov 19, 2009 at 5:45 PM, Dan Yamins wrote: > Hi all: > > I'm just writing to report on my experience using Starcluster, which > enables the use of NumPy and Scipy in the Amazon EC2 cloud computing > environment. The purpose of my email is to extol Starcluster's qualities, > and suggest t