#NEWLINE# def _compute_a( self ) : #NEWLINE# #INDENT# 'fixed effects parameters\n\n        Display (3.1) of\n        Laird, Lange, Stram (see help(Mixed)).\n\n        ' #NEWLINE# for unit in self . units : #NEWLINE# #INDENT# unit . fit ( self . a , self . D , self . sigma ) #NEWLINE# #UNINDENT# S = sum ( [ unit . compute_xtwx ( ) for unit in self . units ] ) #NEWLINE# Y = sum ( [ unit . compute_xtwy ( ) for unit in self . units ] ) #NEWLINE# self . Sinv = L . pinv ( S ) #NEWLINE# self . a = np . dot ( self . Sinv , Y )