# p r 5   l i n e a r   r e g r e s s i o n   a n d   M S E , M A E , r s m 
 
 
 
 i m p o r t   p a n d a s   a s   p d                                 
 
 i m p o r t   n u m p y   a s   n p                                   
 
 i m p o r t   m a t p l o t l i b . p y p l o t   a s   p l t         
 
 f r o m   s k l e a r n . l i n e a r _ m o d e l   i m p o r t   L i n e a r R e g r e s s i o n       
 
 f r o m   s k l e a r n . m e t r i c s   i m p o r t   m e a n _ s q u a r e d _ e r r o r ,   m e a n _ a b s o l u t e _ e r r o r ,   r 2 _ s c o r e     
 
 
 
 
 
 d f   =   p d . r e a d _ c s v ( " t e m p e r a t u r e s . c s v " )       
 
 d a t a   =   d f [ [ ' Y E A R ' ,   ' A N N U A L ' ] ]                   
 
 d a t a   =   d a t a . d r o p n a ( )                                       
 
 p r i n t ( d a t a ) 
 
 
 
 
 
 X   =   d a t a [ ' Y E A R ' ] . v a l u e s . r e s h a p e ( - 1 ,   1 )     
 
 y   =   d a t a [ ' A N N U A L ' ] . v a l u e s     
 
 
 
 
 
 m o d e l   =   L i n e a r R e g r e s s i o n ( )                             
 
 m o d e l . f i t ( X ,   y )                                                 
 
 y _ p r e d   =   m o d e l . p r e d i c t ( X )     
 
 
 
 
 
 m s e   =   m e a n _ s q u a r e d _ e r r o r ( y ,   y _ p r e d )         
 
 m a e   =   m e a n _ a b s o l u t e _ e r r o r ( y ,   y _ p r e d )       
 
 r 2   =   r 2 _ s c o r e ( y ,   y _ p r e d )                                 
 
 p r i n t ( " M S E   = " ,   m s e ) 
 
 p r i n t ( " M A E   = " ,   m a e ) 
 
 p r i n t ( " R    = " ,   r 2 ) 
 
 
 
 
 
 p l t . f i g u r e ( f i g s i z e = ( 1 0 ,   6 ) ) 
 
 p l t . s c a t t e r ( X ,   y ,   c o l o r = ' b l u e ' ,   l a b e l = ' A c t u a l   D a t a ' )         
 
 p l t . p l o t ( X ,   y _ p r e d ,   c o l o r = ' r e d ' ,   l a b e l = ' R e g r e s s i o n   L i n e ' )   
 
 p l t . x l a b e l ( ' Y E A R ' ) 
 
 p l t . y l a b e l ( ' A n n u a l   T e m p e r a t u r e   (  C ) ' ) 
 
 p l t . t i t l e ( ' Y e a r   v s   A n n u a l   T e m p e r a t u r e   ( L i n e a r   R e g r e s s i o n ) ' ) 
 
 p l t . l e g e n d ( ) 
 
 p l t . s h o w ( ) 