Pierre Raybaut | 1 Nov 20:00

[ Python(x,y) ] New release : 2.1.4

Hi all,

Release 2.1.4 is now available on http://www.pythonxy.com.
(Full Edition, Basic Edition, Light Edition, Custom Edition and Update)

Changes history
Version 2.1.4 (11-01-2008)

     * Added:
           o ReportLab 2.2, the PDF generation library
           o Windows explorer integration: added a "Run in interactive 
mode" (python -i) option on Python files contextual menu
     * Updated:
           o The *-components listed below are not included in 
Python(x,y) Update 2.1.4 because of the huge size of this update 
installer which is very close to the Google Code per-file size limit 
(100MB). Even though these updates are quite minor, note that you can 
download them individually on here.
           o NumPy 1.2.1
           o SciPy 0.6.0.2 (minor update regarding deprecation warnings 
with NumPy 1.2.x)
           o matplotlib 0.98.3.3 (new 660-pages PDF documentation)
           o Enthought Tool Suite 3.0.2.3
           o VTK 5.2.0.4
           o *ITK 3.8.0.2
           o *wxPython 2.8.9.1
           o MDP 2.4
           o PySQlite 2.5.0
           o *Eclipse 3.4.1
           o Pydev 1.3.23
(Continue reading)

Frank Lagor | 2 Nov 00:43
Favicon

guassian_kde and kernel regression

This question is probably for Robert Kern, because I believe the he wrote the gaussian_kde class in scipy.stats.kde, however I would very much appreciate a response from anyone else who could help.  My question is:  Is there currently any way to perform weighted kernel density estimation using the gaussian_kde class?   If not, what needs to be done, and how do I get started?

Just for clarity sake-- by weighted KDE  I mean that I have more than just the distribution of points for the density estimate.  I also have an associated probability with each point.  In this case, I believe it becomes a regression problem and I think is referred to as kernel regression.  I would very much like to use the class to perform both KDE and wKDE.

Thanks in advance,
Frank

_______________________________________________
SciPy-user mailing list
SciPy-user <at> scipy.org
http://projects.scipy.org/mailman/listinfo/scipy-user
Anne Archibald | 2 Nov 00:51
Picon

Re: guassian_kde and kernel regression

2008/11/1 Frank Lagor <dfranci <at> seas.upenn.edu>:
> This question is probably for Robert Kern, because I believe the he wrote
> the gaussian_kde class in scipy.stats.kde, however I would very much
> appreciate a response from anyone else who could help.  My question is:  Is
> there currently any way to perform weighted kernel density estimation using
> the gaussian_kde class?   If not, what needs to be done, and how do I get
> started?
>
> Just for clarity sake-- by weighted KDE  I mean that I have more than just
> the distribution of points for the density estimate.  I also have an
> associated probability with each point.  In this case, I believe it becomes
> a regression problem and I think is referred to as kernel regression.  I
> would very much like to use the class to perform both KDE and wKDE.

The class does not support weights right now, but I don't think it
would be very difficult to add them to most parts of the code,
essentially just adding a "weights" optional argument. The automatic
covariance selection would need some rethinking; you'd need to hunt
down some research papers. (That method is only really appropriate for
unimodal distributions anyway.) But it does seem valuable.

Anne
Robert Kern | 2 Nov 09:18
Picon
Gravatar

Re: guassian_kde and kernel regression

On Sat, Nov 1, 2008 at 18:51, Anne Archibald <aarchiba <at> physics.mcgill.ca> wrote:
> 2008/11/1 Frank Lagor <dfranci <at> seas.upenn.edu>:
>> This question is probably for Robert Kern, because I believe the he wrote
>> the gaussian_kde class in scipy.stats.kde, however I would very much
>> appreciate a response from anyone else who could help.  My question is:  Is
>> there currently any way to perform weighted kernel density estimation using
>> the gaussian_kde class?   If not, what needs to be done, and how do I get
>> started?
>>
>> Just for clarity sake-- by weighted KDE  I mean that I have more than just
>> the distribution of points for the density estimate.  I also have an
>> associated probability with each point.  In this case, I believe it becomes
>> a regression problem and I think is referred to as kernel regression.  I
>> would very much like to use the class to perform both KDE and wKDE.
>
> The class does not support weights right now, but I don't think it
> would be very difficult to add them to most parts of the code,
> essentially just adding a "weights" optional argument. The automatic
> covariance selection would need some rethinking; you'd need to hunt
> down some research papers. (That method is only really appropriate for
> unimodal distributions anyway.) But it does seem valuable.

What Anne said.

--

-- 
Robert Kern

"I have come to believe that the whole world is an enigma, a harmless
enigma that is made terrible by our own mad attempt to interpret it as
though it had an underlying truth."
  -- Umberto Eco
Anand Patil | 2 Nov 15:33
Picon
Gravatar

Problem with mkl 10.0.2: undefined symbol: mkl_blas_xdgemm_1_thr_htn

Hi all,


I'm trying to build numpy from svn on Ubuntu with mkl 10.0.2.018, as it looks like mkl 10.0.3 and above won't work with numpy (I got the i_free error). I get the following problem, which I can't find on Google:

In [1]: from numpy import *

In [2]: A=asmatrix(eye(6000))

In [3]: A=A*A
/opt/intel/mkl/10.0.2.018/lib/em64t/libmkl_mc.so: undefined symbol: mkl_blas_xdgemm_1_thr_htn
 /opt/intel/mkl/10.0.2.018/lib/em64t/libmkl_mc.so: undefined symbol: mkl_blas_xdgemm_1_thr_htn

More confusingly, I can't find where the symbol is defined:

$ nm /opt/intel/mkl/10.0.2.018/lib/em64t/*.so | grep mkl_blas_xdgemm_1_thr_htn
nm: /opt/intel/mkl/10.0.2.018/lib/em64t/libmkl.so: File format not recognized

$

Anyone know how I can fix this?

Thanks,
Anand
_______________________________________________
SciPy-user mailing list
SciPy-user <at> scipy.org
http://projects.scipy.org/mailman/listinfo/scipy-user
Matthieu Brucher | 2 Nov 15:35
Picon

Re: Problem with mkl 10.0.2: undefined symbol: mkl_blas_xdgemm_1_thr_htn

Hi,

There is a problem with the MKL, but not with version prior to 10.0.2.
What libraries did you link on ?

Matthieu

2008/11/2 Anand Patil <anand.prabhakar.patil <at> gmail.com>:
> Hi all,
> I'm trying to build numpy from svn on Ubuntu with mkl 10.0.2.018, as it
> looks like mkl 10.0.3 and above won't work with numpy (I got the i_free
> error). I get the following problem, which I can't find on Google:
> In [1]: from numpy import *
> In [2]: A=asmatrix(eye(6000))
> In [3]: A=A*A
> /opt/intel/mkl/10.0.2.018/lib/em64t/libmkl_mc.so: undefined symbol:
> mkl_blas_xdgemm_1_thr_htn
>  /opt/intel/mkl/10.0.2.018/lib/em64t/libmkl_mc.so: undefined symbol:
> mkl_blas_xdgemm_1_thr_htn
> More confusingly, I can't find where the symbol is defined:
> $ nm /opt/intel/mkl/10.0.2.018/lib/em64t/*.so | grep
> mkl_blas_xdgemm_1_thr_htn
> nm: /opt/intel/mkl/10.0.2.018/lib/em64t/libmkl.so: File format not
> recognized
> $
> Anyone know how I can fix this?
> Thanks,
> Anand
> _______________________________________________
> SciPy-user mailing list
> SciPy-user <at> scipy.org
> http://projects.scipy.org/mailman/listinfo/scipy-user
>
>

--

-- 
Information System Engineer, Ph.D.
Website: http://matthieu-brucher.developpez.com/
Blogs: http://matt.eifelle.com and http://blog.developpez.com/?blog=92
LinkedIn: http://www.linkedin.com/in/matthieubrucher
Frank Lagor | 2 Nov 18:43
Favicon

Re: guassian_kde and kernel regression



On Sun, Nov 2, 2008 at 3:18 AM, Robert Kern <robert.kern <at> gmail.com> wrote:
On Sat, Nov 1, 2008 at 18:51, Anne Archibald <aarchiba <at> physics.mcgill.ca> wrote:
> 2008/11/1 Frank Lagor <dfranci <at> seas.upenn.edu>:
>> This question is probably for Robert Kern, because I believe the he wrote
>> the gaussian_kde class in scipy.stats.kde, however I would very much
>> appreciate a response from anyone else who could help.  My question is:  Is
>> there currently any way to perform weighted kernel density estimation using
>> the gaussian_kde class?   If not, what needs to be done, and how do I get
>> started?
>>
>> Just for clarity sake-- by weighted KDE  I mean that I have more than just
>> the distribution of points for the density estimate.  I also have an
>> associated probability with each point.  In this case, I believe it becomes
>> a regression problem and I think is referred to as kernel regression.  I
>> would very much like to use the class to perform both KDE and wKDE.
>
> The class does not support weights right now, but I don't think it
> would be very difficult to add them to most parts of the code,
> essentially just adding a "weights" optional argument. The automatic
> covariance selection would need some rethinking; you'd need to hunt
> down some research papers. (That method is only really appropriate for
> unimodal distributions anyway.) But it does seem valuable.

What Anne said.

--
Robert Kern

"I have come to believe that the whole world is an enigma, a harmless
enigma that is made terrible by our own mad attempt to interpret it as
though it had an underlying truth."
 -- Umberto Eco
_______________________________________________
SciPy-user mailing list
SciPy-user <at> scipy.org
http://projects.scipy.org/mailman/listinfo/scipy-user

Wonderful.  Thank you both very much for your responses.  I will soon get started working on it.
Take care,
Frank

_______________________________________________
SciPy-user mailing list
SciPy-user <at> scipy.org
http://projects.scipy.org/mailman/listinfo/scipy-user
Nils Wagner | 3 Nov 11:40
Picon
Favicon

odeint parameter

Hi all,

I have a question wrt.

odeint(func, y0, t, args=(), Dfun=None, col_deriv=0, 
full_output=0, ml=None, mu=None, rtol=None, atol=None, 
tcrit=None, h0=0.0, hmax=0.0, hmin=0.0, ixpr=0, mxstep=0, 
mxhn
il=0, mxordn=12, mxords=5, printmessg=0)

     h0 : float, (0: solver-determined)
         The step size to be attempted on the first step.
     hmax : float, (0: solver-determined)
         The maximum absolute step size allowed.
     hmin : float, (0: solver-determined)
         The minimum absolute step size allowed.

Is it really useful to start with the default values
h0=0.0, hmax= 0.0 and hmin=0.0 ?

Nils
Anand Patil | 3 Nov 13:10
Picon
Gravatar

Re: Problem with mkl 10.0.2: undefined symbol: mkl_blas_xdgemm_1_thr_htn

Just about every combination of mkl_core, mkl_def, mkl, guide and mkl_intel_thread .


Thanks,
Anand

On Sun, Nov 2, 2008 at 2:35 PM, Matthieu Brucher <matthieu.brucher <at> gmail.com> wrote:
Hi,

There is a problem with the MKL, but not with version prior to 10.0.2.
What libraries did you link on ?

Matthieu

2008/11/2 Anand Patil <anand.prabhakar.patil <at> gmail.com>:
> Hi all,
> I'm trying to build numpy from svn on Ubuntu with mkl 10.0.2.018, as it
> looks like mkl 10.0.3 and above won't work with numpy (I got the i_free
> error). I get the following problem, which I can't find on Google:
> In [1]: from numpy import *
> In [2]: A=asmatrix(eye(6000))
> In [3]: A=A*A
> /opt/intel/mkl/10.0.2.018/lib/em64t/libmkl_mc.so: undefined symbol:
> mkl_blas_xdgemm_1_thr_htn
>  /opt/intel/mkl/10.0.2.018/lib/em64t/libmkl_mc.so: undefined symbol:
> mkl_blas_xdgemm_1_thr_htn
> More confusingly, I can't find where the symbol is defined:
> $ nm /opt/intel/mkl/10.0.2.018/lib/em64t/*.so | grep
> mkl_blas_xdgemm_1_thr_htn
> nm: /opt/intel/mkl/10.0.2.018/lib/em64t/libmkl.so: File format not
> recognized
> $
> Anyone know how I can fix this?
> Thanks,
> Anand
> _______________________________________________
> SciPy-user mailing list
> SciPy-user <at> scipy.org
> http://projects.scipy.org/mailman/listinfo/scipy-user
>
>



--
Information System Engineer, Ph.D.
Website: http://matthieu-brucher.developpez.com/
Blogs: http://matt.eifelle.com and http://blog.developpez.com/?blog=92
LinkedIn: http://www.linkedin.com/in/matthieubrucher
_______________________________________________
SciPy-user mailing list
SciPy-user <at> scipy.org
http://projects.scipy.org/mailman/listinfo/scipy-user

_______________________________________________
SciPy-user mailing list
SciPy-user <at> scipy.org
http://projects.scipy.org/mailman/listinfo/scipy-user
Matthieu Brucher | 3 Nov 13:34
Picon

Re: Problem with mkl 10.0.2: undefined symbol: mkl_blas_xdgemm_1_thr_htn

Hi,

It seems you may have to link against libmkl_sequential.so. I don't
now why, because it depends on the parameters you used, but it may do
the trick.

Matthieu

2008/11/3 Anand Patil <anand.prabhakar.patil <at> gmail.com>:
> Just about every combination of mkl_core, mkl_def, mkl, guide and
> mkl_intel_thread .
> Thanks,
> Anand
>
> On Sun, Nov 2, 2008 at 2:35 PM, Matthieu Brucher
> <matthieu.brucher <at> gmail.com> wrote:
>>
>> Hi,
>>
>> There is a problem with the MKL, but not with version prior to 10.0.2.
>> What libraries did you link on ?
>>
>> Matthieu
>>
>> 2008/11/2 Anand Patil <anand.prabhakar.patil <at> gmail.com>:
>> > Hi all,
>> > I'm trying to build numpy from svn on Ubuntu with mkl 10.0.2.018, as it
>> > looks like mkl 10.0.3 and above won't work with numpy (I got the i_free
>> > error). I get the following problem, which I can't find on Google:
>> > In [1]: from numpy import *
>> > In [2]: A=asmatrix(eye(6000))
>> > In [3]: A=A*A
>> > /opt/intel/mkl/10.0.2.018/lib/em64t/libmkl_mc.so: undefined symbol:
>> > mkl_blas_xdgemm_1_thr_htn
>> >  /opt/intel/mkl/10.0.2.018/lib/em64t/libmkl_mc.so: undefined symbol:
>> > mkl_blas_xdgemm_1_thr_htn
>> > More confusingly, I can't find where the symbol is defined:
>> > $ nm /opt/intel/mkl/10.0.2.018/lib/em64t/*.so | grep
>> > mkl_blas_xdgemm_1_thr_htn
>> > nm: /opt/intel/mkl/10.0.2.018/lib/em64t/libmkl.so: File format not
>> > recognized
>> > $
>> > Anyone know how I can fix this?
>> > Thanks,
>> > Anand
>> > _______________________________________________
>> > SciPy-user mailing list
>> > SciPy-user <at> scipy.org
>> > http://projects.scipy.org/mailman/listinfo/scipy-user
>> >
>> >
>>
>>
>>
>> --
>> Information System Engineer, Ph.D.
>> Website: http://matthieu-brucher.developpez.com/
>> Blogs: http://matt.eifelle.com and http://blog.developpez.com/?blog=92
>> LinkedIn: http://www.linkedin.com/in/matthieubrucher
>> _______________________________________________
>> SciPy-user mailing list
>> SciPy-user <at> scipy.org
>> http://projects.scipy.org/mailman/listinfo/scipy-user
>
>
> _______________________________________________
> SciPy-user mailing list
> SciPy-user <at> scipy.org
> http://projects.scipy.org/mailman/listinfo/scipy-user
>
>

--

-- 
Information System Engineer, Ph.D.
Website: http://matthieu-brucher.developpez.com/
Blogs: http://matt.eifelle.com and http://blog.developpez.com/?blog=92
LinkedIn: http://www.linkedin.com/in/matthieubrucher

Gmane