U
    ĂmśdĐ  ă                	   @   sÚ  d dl Zd dlmZ ddlmZmZ d\ZZ	dZ
dZe de
eˇdd…ejf Ze eˇZe ede e eˇˇ  dˇdd…ejf Zedej edˇ  Zd	Zedde…dd…f Zedde…dd…f Zejeded
 …dd…f eed
 d d…dd…f f Zejeded
 …dd…f eed
 d d…dd…f f Ze e ed
 ˇe ed
 d eˇfˇZeeeeddŤZe  eˇZ!e "ˇ  e #eedee!dˇ e $dˇ e "ˇ  e #ee %ˇ dede!dˇ e $dˇ dS )é    Né   )ÚGaussProcessÚkernel_euclid)é2   é   é   é
   g{®Gáz„?gš™™™™™ą?é   é   gü©ńŇMb@?)ZkernelZ
ridgecoeffZbozr.z4euclid kernel: true y versus noisy y and estimated yzbo-zgo-zr.-zJeuclid kernel: true (green), noisy (blue) and estimated (red) observations)&ÚnumpyÚnpZmatplotlib.pyplotZpyplotZpltZkernridgeregress_classr   r   ÚmÚkÚupperÚscaleZlinspaceZnewaxisZxsÚsinZxs1ÚsumÚsqrtÚabsZy1trueÚrandomZrandnÚy1ZstrideZxstrainZystrainZr_ZhstackZarangeÚindexZgp1ZpredictZyhatr1ZfigureZplotÚtitleZravel© r   r   úi/home/sam/Atlas/atlas_env/lib/python3.8/site-packages/statsmodels/sandbox/regression/example_kernridge.pyÚ<module>   s2   
.::(˙

