��#� �G�R�E�M�L�M�o�d�e�l�s�.�j�l� � ��G�R�E�M�L�M�o�d�e�l�s�.�j�l�� �i�s� �a� �J�u�l�i�a� �p�a�c�k�a�g�e� �f�o�r� �f�i�t�t�i�n�g� �v�a�r�i�a�n�c�e� �c�o�m�p�o�n�e�n�t� �m�o�d�e�l�s� �t�h�a�t� �a�r�e� �s�t�r�u�c�t�u�r�e�d� �a�c�c�o�r�d�i�n�g� �t�o� �k�n�o�w�n�,� �d�e�n�s�e�,� �r�e�l�a�t�i�o�n�s�h�i�p� �m�a�t�r�i�c�e�s�.� � �T�h�e� �i�n�t�e�n�d�e�d� �u�s�e� �c�a�s�e� �i�s� �t�h�e� �a�n�a�l�y�s�i�s� �o�f� �g�e�n�e�t�i�c� �a�n�d� �e�n�v�i�r�o�n�m�e�n�t�a�l� �c�o�m�p�o�n�e�n�t�s� �o�f� �v�a�r�i�a�n�c�e� �f�o�r� �q�u�a�n�t�i�t�a�t�i�v�e� �t�r�a�i�t�s�,� �w�h�e�n� �g�e�n�e�t�i�c� �r�e�l�a�t�i�o�n�s�h�i�p� �m�a�t�r�i�c�e�s� �a�r�e� �c�o�m�p�u�t�e�d� �f�r�o�m� �g�e�n�o�m�e�-�w�i�d�e� �S�N�P� �d�a�t�a�.� � � � �#�#� �M�o�d�e�l�s� � ��G�R�E�M�L�M�o�d�e�l�s�.�j�l�� �c�a�n� �f�i�t� �m�o�d�e�l�s� �o�f� �t�h�e� �f�o�r�m� � �$�$�\�b�o�l�d�s�y�m�b�o�l�{�y�}� �\�s�i�m� �N�(�\�b�f�{�X�}� �\�b�o�l�d�s�y�m�b�o�l�{�\�b�e�t�a�}�,� �\�b�f�{�V�}�)�.�$�$� � � � �I�n� �t�h�i�s� �m�o�d�e�l� �$�\�b�o�l�d�s�y�m�b�o�l�{�\�b�e�t�a�}�$� �i�s� �t�h�e� �f�i�x�e�d� �e�f�f�e�c�t�s� �o�f� �c�o�v�a�r�i�a�t�e�s� �$�\�b�f�{�X�}�$� �a�n�d� �$�\�b�f�{�V�}�$� �i�s� �t�h�e� �c�o�v�a�r�i�a�n�c�e� �m�a�t�r�i�x� �o�f� �t�h�e� �c�o�n�d�i�t�i�o�n�a�l� �r�e�s�p�o�n�s�e�s�.� �$�\�b�f�{�V�}�$� �c�a�n� �b�e� �m�o�d�e�l�l�e�d� �w�i�t�h� �t�h�e� �f�o�l�l�o�w�i�n�g� �s�t�r�u�c�t�u�r�e� � � � �$�$�\�b�f�{�V�}� �=� �\�s�u�m�_�{�i�=�1�}�^�q� �\�d�e�l�t�a�_�i� �\�b�f�{�R�}�_�i�$�$� � � � �w�h�e�r�e� �$�\�d�e�l�t�a�_�i�$� �a�r�e� �v�a�r�i�a�n�c�e� �c�o�m�p�o�n�e�n�t� �p�a�r�a�m�e�t�e�r�s� �a�n�d� �$�\�b�f�{�R�}�_�i�$� �a�r�e� �s�y�m�m�e�t�r�i�c� �m�a�t�r�i�c�e�s� �t�h�a�t� �a�r�e� �p�r�o�v�i�d�e�d� �b�y� �t�h�e� �u�s�e�r�.� � � � �T�h�e� �v�a�r�i�a�n�c�e� �c�o�m�p�o�n�e�n�t� �p�a�r�a�m�e�t�e�r�s� �c�a�n� �b�e� �d�e�f�i�n�e�d� �a�s� �f�u�n�c�t�i�o�n�s� �o�f� �t�h�e� �v�e�c�t�o�r� �o�f� �p�a�r�a�m�e�t�e�r�s� �$�\�b�f�{�\�t�h�e�t�a�}�$� �t�h�a�t� �a�r�e� �o�p�t�i�m�i�z�e�d�.� �T�h�i�s� �m�a�y� �f�o�r� �i�n�s�t�a�n�c�e� �b�e� �n�e�e�d�e�d� �i�f� �$�\�d�e�l�t�a�_�i�$� �i�s� �a� �c�o�v�a�r�i�a�n�c�e� �p�a�r�a�m�e�t�e�r� �t�h�a�t� �i�s� �b�o�u�n�d�e�d� �b�y� �t�h�e� �v�a�l�u�e�s� �o�f� �o�t�h�e�r� �p�a�r�a�m�e�t�e�r�s�.� �T�h�i�s� �w�a�s� �t�h�e� �m�o�t�i�v�a�t�i�o�n� �f�o�r� �c�r�e�a�t�i�n�g� �t�h�i�s� �p�a�c�k�a�g�e�.� � � � �#�#� �I�n�s�t�a�l�l�a�t�i�o�n� � �T�h�e� �p�a�c�k�a�g�e� �c�a�n� �b�e� �i�n�s�t�a�l�l�e�d� �f�r�o�m� �g�i�t�h�u�b� �w�i�t�h� � �� �j�u�l�i�a� � �(�@�v�1�.�6�)� �p�k�g�>� �a�d�d� �h�t�t�p�s�:�/�/�g�i�t�h�u�b�.�c�o�m�/�e�s�p�e�n�m�e�i�/�G�R�E�M�L�M�o�d�e�l�s�.�j�l� � �� � � � �#�#� �E�x�a�m�p�l�e� �1� � �H�e�r�e� �i�s� �a�n� �e�x�a�m�p�l�e� �w�i�t�h� �s�i�m�u�l�a�t�e�d� �d�a�t�a� �s�h�o�w�i�n�g� �h�o�w� �t�h�e� �p�a�c�k�a�g�e� �m�a�y� �b�e� �u�s�e�d� �t�o� �f�i�t� �t�w�o� �v�a�r�i�a�n�c�e� �c�o�m�p�o�n�e�n�t�s� �a�n�d� �t�w�o� �c�o�v�a�r�i�a�t�e�s� �f�o�r� �t�h�e� �m�e�a�n�s�.� �T�h�e� �m�o�d�e�l� �f�o�r� �t�h�e� �m�e�a�n�s� �c�a�n� �b�e� �d�e�f�i�n�e�d� �u�s�i�n�g� �t�h�e� ��@�f�o�r�m�u�l�a�� �m�a�c�r�o� �t�o�g�e�t�h�e�r� �w�i�t�h� �a� ��D�a�t�a�F�r�a�m�e�� �s�t�o�r�i�n�g� �t�h�e� �c�o�v�a�r�i�a�t�e�s� �a�n�d� �t�h�e� �r�e�s�p�o�n�s�e� �v�a�r�i�a�b�l�e�.� �T�h�e� �m�o�d�e�l� �f�o�r� �t�h�e� �c�o�v�a�r�i�a�n�c�e� �c�a�n� �b�e� �d�e�f�i�n�e�d� �b�y� �p�r�o�v�i�d�i�n�g� �a� �v�e�c�t�o�r� �w�i�t�h� �r�e�l�a�t�i�o�n�s�h�i�p� �m�a�t�r�i�c�e�s�.� � ��j�u�l�i�a� � �u�s�i�n�g� �G�R�E�M�L�M�o�d�e�l�s�,� �L�i�n�e�a�r�A�l�g�e�b�r�a�,� �D�a�t�a�F�r�a�m�e�s�,� �S�t�a�t�s�M�o�d�e�l�s� � � � �n�,� �m� �=� �1�0�0�0�,� �2�0�0�0� � �v� �=� �r�a�n�d�n�(�n�,� �m�)� � �R�1� �=� �S�y�m�m�e�t�r�i�c�(�v� �*� �v�'� �/� �m�)� � �R�2� �=� �D�i�a�g�o�n�a�l�(�o�n�e�s�(�n�)�)� � �V� �=� �2�R�1� �+� �4�R�2� � �y� �=� �c�h�o�l�e�s�k�y�(�V�)�.�L� �*� �r�a�n�d�n�(�n�)� � �d�a�t� �=� �D�a�t�a�F�r�a�m�e�(�y� �=� �y�,� �x� �=� �r�a�n�d�(�n�)�,� �z� �=� �r�a�n�d�(�n�)�)� � �r� �=� �[�R�1�,� �R�2�]� � �m� �=� �f�i�t�(�G�R�E�M�L�M�o�d�e�l�,� �@�f�o�r�m�u�l�a�(�y� �~� �1� �+� �x� �+� �z�)�,� �d�a�t�,� �r�)� � �� � �T�h�i�s� �s�h�o�u�l�d� �g�i�v�e� �a�n� �o�u�t�p�u�t� �r�e�a�s�o�n�a�b�l�y� �s�i�m�i�l�a�r� �t�o� �t�h�i�s� � ��j�u�l�i�a� � �l�o�g�L�i�k� � � � � �-�2� �l�o�g�L�i�k� � �A�I�C� � � � � � � � �A�I�C�c� � � � � � � �B�I�C� � �-�2�3�1�4�.�5�5�9�0� �4�6�2�9�.�1�1�8�0� � �4�6�3�9�.�1�1�8�0� � �4�6�3�9�.�1�7�8�4� � �4�6�6�3�.�6�5�6�8� � � � � � � �V�a�r�i�a�n�c�e� �c�o�m�p�o�n�e�n�t� �p�a�r�a�m�e�t�e�r�s�:� � �C�o�m�p�.� � � �E�s�t�.� � � � �S�t�d�.� �E�r�r�o�r� � ���� � � � � � �2�.�0�9�4�3� � �-� � ���� � � � � � �4�.�0�6�7�4� � �-� � � � � �F�i�x�e�d�-�e�f�f�e�c�t�s� �p�a�r�a�m�e�t�e�r�s�:� � ��%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�% � � � � � � � � � � � � � � � � � � � �C�o�e�f�.� � �S�t�d�.� �E�r�r�o�r� � � � � � �z� � �P�r�(�>�|�z�|�)� � ��%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�% � �(�I�n�t�e�r�c�e�p�t�)� � �-�0�.�1�7�9�0�7�1� � � � � �0�.�2�0�3�2�1�3� � �-�0�.�8�8� � � � �0�.�3�7�8�2� � �x� � � � � � � � � � � � �-�0�.�0�0�2�9�3�2�7� � � � �0�.�2�6�2�8�4� � � �-�0�.�0�1� � � � �0�.�9�9�1�1� � �z� � � � � � � � � � � � � �0�.�4�3�5�4�4�5� � � � � �0�.�2�7�4�5�3�8� � � �1�.�5�9� � � � �0�.�1�1�2�7� � ��%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�% � �� � � � �N�o�t� �a�l�l� �a�r�g�u�m�e�n�t�s� �a�r�e� �a�v�a�i�l�a�b�l�e� �w�h�e�n� �f�i�t�t�i�n�g� �a� �m�o�d�e�l� �t�h�i�s� �w�a�y�.� �S�o�m�e�t�i�m�e�s� �i�t� �m�a�y� �b�e� �u�s�e�f�u�l� �t�o� �f�i�r�s�t� �d�e�f�i�n�e� �a� �m�o�d�e�l� �w�i�t�h� �t�h�e� �r�e�l�e�v�a�n�t� �a�r�g�u�m�e�n�t�s� �a�n�d� �i�n� �a� �n�e�x�t� �s�t�e�p� �t�o� �f�i�t� �t�h�a�t� �m�o�d�e�l�.� � � � ��j�u�l�i�a� � �y� �=� �d�a�t�[�!�,�:�y�]� � �X� �=� �M�a�t�r�i�x�(�h�c�a�t�(�o�n�e�s�(�n�)�,� �d�a�t�[�!�,�2�:�3�]�)�)� � �d�a�t� �=� �G�R�E�M�L�D�a�t�a�(�y�,� �X�,� �r�)� � �m�2� �=� �G�R�E�M�L�M�o�d�e�l�(�d�a�t�,� �[�1�.�0�,� �1�.�0�]�,� �[�0�.�0�,� �0�.�0�]�)� � �f�i�t�!�(�m�2�)� � �� � �T�h�i�s� �f�i�r�s�t� �c�r�e�a�t�e�s� �a�n� �o�b�j�e�c�t� �o�f� �t�y�p�e� ��G�R�E�M�L�D�a�t�a�� �t�h�a�t� �i�s� �u�s�e�d� �t�o� �h�o�l�d� �t�h�e� �d�a�t�a� �i�n�p�u�t� �t�o� �a� ��G�R�E�M�L�M�o�d�e�l��.� �T�h�e�n� �a� �m�o�d�e�l�,� �r�e�p�r�e�s�e�n�t�e�d� �b�y� �a�n� �o�b�j�e�c�t� �o�f� �t�y�p�e� ��G�R�E�M�L�M�o�d�e�l�� �i�s� �c�r�e�a�t�e�d� �f�r�o�m� �a� �c�o�n�s�t�r�u�c�t�o�r� �t�h�a�t� �a�c�c�e�p�t�s� �a� ��G�R�E�M�L�D�a�t�a��,� �a� �v�e�c�t�o�r� �o�f� �i�n�i�t�i�a�l� �v�a�l�u�e�s� �f�o�r� �t�h�e� �v�a�r�i�a�n�c�e� �c�o�m�p�o�n�e�n�t� �p�a�r�a�m�e�t�e�r�s� ��[�1�.�0�,� �1�.�0�]�� �a�n�d� �a� �v�e�c�t�o�r� �o�f� �l�o�w�e�r� �b�o�u�n�d�s� �f�o�r� �t�h�o�s�e� �p�a�r�a�m�e�t�e�r�s� ��[�0�.�0�,� �0�.�0�]��.� �A�n�d� �i�n� �t�h�e� �l�a�s�t� �l�i�n�e� �t�h�a�t� �m�o�d�e�l� �i�s� �f�i�t�t�e�d� ��f�i�t�!�(�m�2�)��.� �T�h�e� ��!�� �s�i�g�n�a�l�s� �t�h�a�t� �t�h�e� �f�u�n�c�t�i�o�n� �m�o�d�i�f�i�e�s� �t�h�e� �m�o�d�e�l�,� �w�h�i�c�h� �i�s� �a� �c�o�n�v�e�n�t�i�o�n� �i�n� �j�u�l�i�a�.� � � � �#�#� �E�x�a�m�p�l�e� �2� � �C�o�n�s�i�d�e�r� �t�h�e� �m�o�d�e�l� �f�o�r� �r�e�s�p�o�n�s�e�s� �o�b�t�a�i�n�e�d� �f�r�o�m� �o�f�f�s�p�r�i�n�g� � � � �$�$�\�b�o�l�d�s�y�m�b�o�l�{�y�}� �=� �\�b�f�{�X�}� �\�b�o�l�d�s�y�m�b�o�l�{�\�b�e�t�a�}� �+� �\�b�o�l�d�s�y�m�b�o�l�{�m�}� �+� �\�b�o�l�d�s�y�m�b�o�l�{�p�}� �+� �\�b�o�l�d�s�y�m�b�o�l�{�o�}� �+� �\�b�o�l�d�s�y�m�b�o�l�{�\�e�p�s�i�l�o�n�}�,�$�$� � � � �w�h�e�r�e� �$�\�b�o�l�d�s�y�m�b�o�l�{�o�}�$� �a�r�e� �d�i�r�e�c�t� �a�d�d�i�t�i�v�e� �e�f�f�e�c�t�s� �o�f� �o�f�f�s�p�r�i�n�g� �o�w�n� �g�e�n�e�s�,� �$�\�b�o�l�d�s�y�m�b�o�l�{�m�}�$� �a�n�d� �$�\�b�o�l�d�s�y�m�b�o�l�{�p�}�$� �a�r�e� �i�n�d�i�r�e�c�t� �e�f�f�e�c�t�s� �o�f� �m�a�t�e�r�n�a�l� �a�n�d� �p�a�t�e�r�n�a�l� �g�e�n�e�s� �a�n�d� �$�\�b�o�l�d�s�y�m�b�o�l�{�\�e�p�s�i�l�o�n�}�$� �e�r�r�o�r�s�.� �P�a�r�a�m�e�t�e�r�s� �o�f� �t�h�e� �m�o�d�e�l� �a�r�e� �t�h�e� �v�a�r�i�a�n�c�e�s� �a�n�d� �c�o�v�a�r�i�a�n�c�e�s� �a�m�o�n�g� �t�h�e� �g�e�n�e�t�i�c� �e�f�f�e�c�t�s� �a�n�d� �t�h�e� �e�r�r�o�r� �v�a�r�i�a�n�c�e�.� �A�m�o�n�g� �i�n�d�i�v�i�d�u�a�l�s� �t�h�e� �g�e�n�e�t�i�c� �e�f�f�e�c�t�s� �a�r�e� �c�o�r�r�e�l�a�t�e�d� �a�c�c�o�r�d�i�n�g� �t�o� �t�h�e�i�r� �g�e�n�e�t�i�c� �r�e�l�a�t�e�d�n�e�s�s�.� �C�o�n�d�i�t�i�o�n�a�l� �o�n� �t�h�e� �c�o�v�a�r�i�a�t�e�s�,� �t�h�i�s� �m�o�d�e�l� �h�a�s� �c�o�v�a�r�i�a�n�c�e� �m�a�t�r�i�x� � � � �$�$�\�m�a�t�h�b�f�{�V� �=� �\�s�i�g�m�a�_�o�^�2�{�A�}�_�{�m�m�}� �+� �\�s�i�g�m�a�_�m�^�2�A�_�{�p�p�}� �+� �\�s�i�g�m�a�_�p�^�2�A�_�{�o�o�}�\�\� � � �+� �\�s�i�g�m�a�_�{�m�p�}�(�A�_�{�m�p�}�+�A�_�{�p�m�}�)�+� �\�s�i�g�m�a�_�{�m�o�}�(�A�_�{�m�o�}� �+� �A�_�{�o�m�}�)� �+� �\�s�i�g�m�a�_�{�p�o�}�(�A�_�{�p�o�}�+�A�_�{�o�p�}�)� �+� �\�s�i�g�m�a�_�\�e�p�s�i�l�o�n�^�2�I�}�.�$�$� � � � �T�h�e� �g�e�n�e�t�i�c� �r�e�l�a�t�e�d�n�e�s�s� �m�a�t�r�i�c�e�s� �$�\�b�f�{�A�}�_�{�i�j�}�$� �a�r�e� �b�l�o�c�k�s� �o�f� �t�h�e� �g�e�n�e�t�i�c� �r�e�l�a�t�e�d�n�e�s�s� �m�a�t�r�i�x� �o�f� �a�l�l� �i�n�d�i�v�i�d�u�a�l�s� �w�h�e�n� �s�t�r�u�c�t�u�r�e�d� �a�c�c�o�r�d�i�n�g�l�y� � � � �$�$�\�b�f�{�A�}� �=� �\�b�e�g�i�n�{�b�m�a�t�r�i�x�}� �\�b�f�{�A�}�_�{�m�m�}� �&� �\�b�f�{�A�}�_�{�m�p�}� �&� �\�b�f�{�A�}�_�{�m�o�}�\�\� � �\�b�f�{�A�}�_�{�p�m�}� �&� �\�b�f�{�A�}�_�{�p�p�}� �&� �\�b�f�{�A�}�_�{�p�o�}�\�\� � �\�b�f�{�A�}�_�{�o�m�}� �&� �\�b�f�{�A�}�_�{�o�p�}� �&� �\�b�f�{�A�}�_�{�o�o�}� �\�e�n�d�{�b�m�a�t�r�i�x�}�.�$�$� � � � �I�n� �o�r�d�e�r� �t�o� �f�i�t� �t�h�e� �m�o�d�e�l�,� �t�h�e� �f�i�r�s�t� �s�t�e�p� �i�s� �t�o� �c�o�m�p�u�t�e� �t�h�e� �g�n�e�e�t�i�c� �r�e�l�a�t�i�o�n�s�h�i�p� �m�a�t�r�i�c�e�s�.� �T�h�e� �g�e�n�e�t�i�c� �r�e�l�a�t�i�o�n�s�h�i�p�s� �w�i�l�l� �b�e� �e�s�t�i�m�a�t�e�d� �f�r�o�m� �a� �s�i�m�u�l�a�t�e�d� �d�a�t�a�s�e�t� �o�f� �1�0�0�0�0� �S�N�P�s� �f�r�o�m� �6�0�0�0� �p�a�r�e�n�t�-�o�f�f�s�p�r�i�n�g� �t�r�i�o�s� �s�t�o�r�e�d� �i�n� ��p�l�i�n�k�� �f�o�r�m�a�t�.� �T�h�e� �j�u�l�i�a� �p�a�c�k�a�g�e� ��S�n�p�A�r�r�a�y�s�� �i�s� �u�s�e�d� �t�o� �r�e�a�d� �s�n�p�-�d�a�t�a� �a�n�d� �c�o�m�p�u�t�e� �g�e�n�e�t�i�c� �r�e�l�a�t�i�o�n�s�h�i�p�s�.� � � � ��j�u�l�i�a� � �u�s�i�n�g� �V�C�M�o�d�e�l�s�,� �L�i�n�e�a�r�A�l�g�e�b�r�a�,� �D�a�t�a�F�r�a�m�e�s�,� �S�t�a�t�s�M�o�d�e�l�s�,� �S�n�p�A�r�r�a�y�s� � � � �g�p�a�t�h� �=� �j�o�i�n�p�a�t�h�(�@�_�_�D�I�R�_�_�,� �"�d�a�t�a�"�,� �"�t�r�i�o�"�)� � �t�r�i�o� �=� �S�n�p�D�a�t�a�(�g�p�a�t�h�)� � �A� �=� �2� �*� �g�r�m�(�t�r�i�o�.�s�n�p�a�r�r�a�y�)� � �� � � � �B�e�c�a�u�s�e� �t�h�e� �s�n�p�-�d�a�t�a� �w�a�s� �a�r�r�a�n�g�e�d� �a�c�c�o�r�d�i�n�g� �t�o� �f�a�m�i�l�i�e�s� �w�i�t�h�i�n� �r�o�l�e�s�,� �t�h�e� �g�e�n�e�t�i�c� �r�e�l�a�t�i�o�n�s�h�i�p� �m�a�t�r�i�x� �c�a�n� �b�e� �p�a�r�t�i�t�i�o�n�e�d� �a�s� � � � ��j�u�l�i�a� � �K� �=� �d�i�v�(�t�r�i�o�.�p�e�o�p�l�e�,� �3�)� � �m�i�d� �=� �1�:�K� � �p�i�d� �=� �K�+�1�:�2�K� � �o�i�d� �=� �2�K�+�1�:�3�*�K� � � � �A�m�m� �=� �A�[�m�i�d�,� �m�i�d�]� � �A�p�p� �=� �A�[�p�i�d�,� �p�i�d�]� � �A�o�o� �=� �A�[�o�i�d�,� �o�i�d�]� � �D�p�m� �=� �A�[�p�i�d�,� �m�i�d�]� �+� �A�[�m�i�d�,� �p�i�d�]� � �D�o�m� �=� �A�[�o�i�d�,� �m�i�d�]� �+� �A�[�m�i�d�,� �o�i�d�]� � �D�o�p� �=� �A�[�o�i�d�,� �p�i�d�]� �+� �A�[�p�i�d�,� �o�i�d�]� � � �� � � � �T�h�e� �o�f�f�-�d�i�a�g�o�n�a�l� �b�l�o�c�k�s� �a�r�e� �a�d�d�e�d� �t�o� �c�o�m�p�u�t�e� �t�h�e� �e�x�p�r�e�s�s�i�o�n�s� �w�i�t�h�i�n� �p�a�r�e�n�t�h�e�s�i�s� �i�n� �t�h�e� �c�o�v�a�r�i�a�n�c�e� �e�x�p�r�e�s�s�i�o�n� �a�b�o�v�e�.� �T�h�e�s�e� �r�e�l�a�t�i�o�n�s�h�i�p� �m�a�t�r�i�c�e�s� �a�r�e� �c�o�l�l�e�c�t�e�d� �i�n� �a� �v�e�c�t�o�r� �a�s� �b�e�f�o�r�e� � � � ��j�u�l�i�a� � �R� �=� �D�i�a�g�o�n�a�l�(�o�n�e�s�(�K�)�)� � �r� �=� �[�A�m�m�,� �A�p�p�,� �A�o�o�,� �D�p�m�,� �D�o�m�,� �D�o�p�,� �R�]� � �� � � � �I�t� �i�s� �a�l�s�o� �n�e�c�e�s�s�a�r�y� �t�o� �d�e�f�i�n�e� �a� �m�a�t�r�i�x� �o�f� �c�o�v�a�r�i�a�t�e�s� �f�o�r� �t�h�e� �m�e�a�n�s�,� �w�h�i�c�h� �i�s� �j�u�s�t� �a� �c�o�n�s�t�a�n�t� �h�e�r�e� � � � ��j�u�l�i�a� � �X� �=� �o�n�e�s�(�K�)� � �� � �T�h�e� �d�a�t�a�s�e�t� �a�l�s�o� �c�o�n�t�a�i�n�s� �s�o�m�e� �s�i�m�u�l�a�t�e�d� �r�e�s�p�o�n�s�e�s� � � � ��j�u�l�i�a� � �y� �=� �p�a�r�s�e�.�(�F�l�o�a�t�6�4�,� �t�r�i�o�.�p�e�r�s�o�n�_�i�n�f�o�[�o�i�d�,� �:�p�h�e�n�o�t�y�p�e�]�)� � �� � � � �A� �c�h�a�l�l�e�n�g�e� �w�i�t�h� �t�h�i�s� �m�o�d�e�l� �i�s� �t�h�a�t� �s�o�m�e� �o�f� �t�h�e� �v�a�r�i�a�n�c�e� �t�e�r�m�s� �r�e�p�r�e�s�e�n�t� �c�o�v�a�r�i�a�n�c�e�s�,� �w�h�i�c�h� �p�u�t�s� �r�e�s�t�r�i�c�t�i�o�n�s� �o�n� �w�h�i�c�h� �v�a�l�u�e�s� �t�h�e� �p�a�r�a�m�e�t�e�r�s� �c�a�n� �t�a�k�e�.� �I�f� �w�e� �c�o�l�l�e�c�t� �t�h�e� �g�e�n�e�t�i�c� �p�a�r�a�m�e�t�e�r�s� �i�n� �a� �c�o�v�a�r�i�a�n�c�e� �m�a�t�r�i�x� � � � �$�$�\�b�f�{�G�}� �=� �\�b�e�g�i�n�{�b�m�a�t�r�i�x�}� �\�s�i�g�m�a�^�2�_�m� �&� �\�s�i�g�m�a�_�{�m�p�}� �&� �\�s�i�g�m�a�_�{�m�o�}�\�\� � �\�s�i�g�m�a�_�{�p�m�}� �&� �\�s�i�g�m�a�^�2�_�p� �&� �\�s�i�g�m�a�_�{�p�o�}�\�\� � �\�s�i�g�m�a�_�{�o�m�}� �&� �\�s�i�g�m�a�_�{�o�p�}� �&� �\�s�i�g�m�a�^�2�_�o� �\�e�n�d�{�b�m�a�t�r�i�x�}�,�$�$� � � � �t�h�i�s� �s�h�o�u�l�d� �b�e� �p�o�s�i�t�i�v�e� �s�e�m�i�d�e�f�i�n�i�t�e�.� �O�n�e� �w�a�y� �t�o� �e�n�f�o�r�c�e� �t�h�i�s� �i�s� �t�o� �l�e�t� �t�h�e� �p�a�r�a�m�e�t�e�r�s� �o�f� �t�h�e� �m�o�d�e�l� �b�e� �t�h�e� �e�l�e�m�e�n�t�s� �o�f� �t�h�e� �l�o�w�e�r�-�t�r�i�a�n�g�u�l�a�r� �C�h�o�l�e�s�k�y� �f�a�c�t�o�r�i�z�a�t�i�o�n� �$�\�b�f�{�L�}�$� �o�f� �t�h�e� �c�o�v�a�r�i�a�n�c�e� �m�a�t�r�i�x� �f�o�r� �g�e�n�e�t�i�c� �e�f�f�e�c�t�s� � � � �$�$�\�b�f�{�L�}� �\�b�f�{�L�}�'� �=� �\�b�f�{�G�}�.�$�$� � � � �I�f� �$�\�b�o�l�d�s�y�m�b�o�l�{�\�t�h�e�t�a�}�$� �i�s� �a� �v�e�c�t�o�r� �o�f� �u�n�i�q�u�e� �p�a�r�a�m�e�t�e�r�s� �o�f� �t�h�e� �m�o�d�e�l�,� �$�\�b�f�{�L�}�$� �c�a�n� �b�e� �d�e�s�c�r�i�b�e�d� �a�s� � � � �$�$�\�b�f�{�L�}� �=� �\�b�e�g�i�n�{�b�m�a�t�r�i�x�}� �\�t�h�e�t�a�_�1� �&� �0� �&� �0�\�\� � �\�t�h�e�t�a�_�2� �&� �\�t�h�e�t�a�_�4� �&� �0�\�\� � �\�t�h�e�t�a�_�3� �&� �\�t�h�e�t�a�_�5� �&� �\�t�h�e�t�a�_�6� �\�e�n�d�{�b�m�a�t�r�i�x�}�.�$�$� � � � �T�h�e�n� � � � �$�$�\�b�f�{�G�}� �=� �\�b�e�g�i�n�{�b�m�a�t�r�i�x�}� �\�s�i�g�m�a�^�2�_�m� �&� �\�s�i�g�m�a�_�{�m�p�}� �&� �\�s�i�g�m�a�_�{�m�o�}�\�\� � �\�s�i�g�m�a�_�{�p�m�}� �&� �\�s�i�g�m�a�^�2�_�p� �&� �\�s�i�g�m�a�_�{�p�o�}�\�\� � �\�s�i�g�m�a�_�{�o�m�}� �&� �\�s�i�g�m�a�_�{�o�p�}� �&� �\�s�i�g�m�a�^�2�_�o� �\�e�n�d�{�b�m�a�t�r�i�x�}� �=� � �\�b�e�g�i�n�{�b�m�a�t�r�i�x�}� �\�t�h�e�t�a�^�2�_�1� �&� �\�t�h�e�t�a�_�1�\�t�h�e�t�a�_�2� �&� �\�t�h�e�t�a�_�1�\�t�h�e�t�a�_�3�\�\� � �\�t�h�e�t�a�_�2�\�t�h�e�t�a�_�1� �&� �\�t�h�e�t�a�^�2�_�2� �+� �\�t�h�e�t�a�^�2�_�4� �&� �\�t�h�e�t�a�_�2�\�t�h�e�t�a�_�3� �+� �\�t�h�e�t�a�_�4�\�t�h�e�t�a�_�5�\�\� � �\�t�h�e�t�a�_�3�\�t�h�e�t�a�_�1� �&� �\�t�h�e�t�a�_�3�\�t�h�e�t�a�_�2� �+� �\�t�h�e�t�a�_�5�\�t�h�e�t�a�_�4� �&� �\�t�h�e�t�a�^�2�_�3� �+� �\�t�h�e�t�a�^�2�_�5� �+� �\�t�h�e�t�a�^�2�_�6�\�e�n�d�{�b�m�a�t�r�i�x�}�$�$� � � � �T�h�i�s� �c�a�n� �b�e� �i�m�p�l�e�m�e�n�t�e�d� �b�y� �o�v�e�r�r�i�d�i�n�g� �t�h�e� �d�e�f�a�u�l�t� ��V�C�M�o�d�e�l�s�.�t�r�a�n�s�f�o�r�m�!�� �m�e�t�h�o�d�.� �T�h�i�s� �i�s� �u�s�e�d� �t�o� �d�e�f�i�n�e� �t�h�e� �$�q�$� �d�i�m�e�n�s�i�o�n�a�l� �v�e�c�t�o�r� �o�f� �v�a�r�i�a�n�c�e� �c�o�m�p�o�n�e�n�t� �p�a�r�a�m�e�t�e�r�s� �$�\�d�e�l�t�a�$� �a�s� �a� �f�u�n�c�t�i�o�n� �o�f� �t�h�e� �$�q�$� �d�i�m�e�n�s�i�o�n�a�l� �v�e�c�t�o�r� �o�f� �p�a�r�a�m�e�t�e�r�s� �t�h�a�t� �a�r�e� �o�p�t�i�m�i�z�e�d� �$�\�t�h�e�t�a�$�.� �I�n� �t�h�i�s� �c�a�s�e�,� �t�h�i�s� �c�o�u�l�d� �b�e� �i�m�p�l�e�m�e�n�t�e�d� �a�s� � � � ��j�u�l�i�a� � �f�u�n�c�t�i�o�n� �G�R�E�M�L�M�o�d�e�l�s�.�t�r�a�n�s�f�o�r�m�!�(���:�:�V�e�c�t�o�r�,� ���:�:�V�e�c�t�o�r�)� � � � � � ���[�1�]� �=� ���[�1�]�^�2� � � � � � ���[�2�]� �=� ���[�2�]�^�2� �+� ���[�4�]�^�2� � � � � � ���[�3�]� �=� ���[�3�]�^�2� �+� ���[�5�]�^�2� �+� ���[�6�]�^�2� � � � � � ���[�4�]� �=� ���[�2�]� �*� ���[�1�]� � � � � � ���[�5�]� �=� ���[�3�]� �*� ���[�1�]� � � � � � ���[�6�]� �=� ���[�3�]� �*� ���[�2�]� �+� ���[�5�]� �*� ���[�4�]� � � � � � ���[�7�]� �=� ���[�7�]� � �e�n�d� � �� � � � �T�h�e� �l�a�s�t� �e�l�e�m�e�n�t� �i�s� �t�h�e� �r�e�s�i�d�u�a�l� �v�a�r�i�a�n�c�e�.� �W�e� �c�a�n� �t�r�y� �t�o� �f�i�t� �t�h�i�s� �m�o�d�e�l� �b�y� �i�n�i�t�i�a�l�i�z�i�n�g� �L� �t�o� �i�d�e�n�t�i�t�y� �a�n�d� �s�e�t�t�i�n�g� �t�h�e� �l�o�w�e�r� �b�o�u�n�d� �f�o�r� �t�h�e� �d�i�a�g�o�n�a�l� �e�l�e�m�e�n�t�s� �t�o� �z�e�r�o� � � � ��j�u�l�i�a� � �d�a�t� �=� �G�R�E�M�L�D�a�t�a�(�y�,� �X�,� �r�)� � �m�1� �=� �G�R�E�M�L�M�o�d�e�l�(�d�a�t�,� �[�1�.�0�,� �0�.�0�,� �0�.�0�,� �1�.�0�,� �0�.�0�,� �1�.�0�,� �1�.�0�]�,� �[�0�.�0�,� �-�I�n�f�,� �-�I�n�f�,� �0�.�0�,� �-�I�n�f�,� �0�.�0�,� �0�.�0�]�)� � �f�i�t�!�(�m�1�)� � �� � � � �T�h�i�s� �s�h�o�u�l�d� �g�i�v�e� �a�n� �o�u�t�p�u�t� �s�i�m�i�l�a�r� �t�o� � � � ��j�u�l�i�a� � �l�o�g�L�i�k� � � � � � �-�2� �l�o�g�L�i�k� � � �A�I�C� � � � � � � � � �A�I�C�c� � � � � � � � �B�I�C� � �-�1�7�7�7�2�.�7�0�3�2� �3�5�5�4�5�.�4�0�6�5� � �3�5�5�6�1�.�4�0�6�5� � �3�5�5�6�1�.�4�3�0�5� � �3�5�6�1�5�.�0�0�2�6� � � � � �V�a�r�i�a�n�c�e� �c�o�m�p�o�n�e�n�t� �p�a�r�a�m�e�t�e�r�s�:� � �C�o�m�p�.� � � �E�s�t�.� � � � �S�t�d�.� �E�r�r�o�r� � ���� � � � � � �1�.�9�2�1�4� � �-� � ���� � � � � � �0�.�1�8�4�8� � �-� � ���� � � � � � �1�.�1�3�4�6� � �-� � ���� � � � � � �-�1�.�6�8�1�4� �-� � ���� � � � � � �-�1�.�2�5�3�0� �-� � ���� � � � � � �0�.�4�9�4�7� � �-� � ���� � � � � � �9�.�3�3�0�9� � �-� � � � � �F�i�x�e�d�-�e�f�f�e�c�t�s� �p�a�r�a�m�e�t�e�r�s�:� � ��%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�% � � � � � � � �C�o�e�f�.� � �S�t�d�.� �E�r�r�o�r� � � � � � �z� � �P�r�(�>�|�z�|�)� � ��%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�% � �x�1� � �2�.�0�7�9�0�9� � � �0�.�0�4�4�8�4�6�2� � �4�6�.�3�6� � � � �<�1�e�-�9�9� � ��%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�%�% � �� � � � �N�o�t�e� �t�h�a�t� �t�h�e� �v�a�r�i�a�n�c�e� �c�o�m�p�o�n�e�n�t�s� �r�e�p�o�r�t�e�d� �i�n� �t�h�e� �o�u�t�p�u�t� �a�r�e� �n�o�t� �t�h�e� �v�a�r�i�a�n�c�e� �c�o�m�p�o�n�e�n�t�s� �i�n� �t�h�e� �m�o�d�e�l�,� �t�h�e�y� �a�r�e� �t�h�e� �e�l�e�m�e�n�t�s� �o�f� �$�\�b�f�{�L�}�$�.� �T�o� �o�b�t�a�i�n� �t�h�e� �v�a�r�i�a�n�c�e� �c�o�m�p�o�n�e�n�t�s� �w�e� �c�o�u�l�d� �d�o� � � � ��j�u�l�i�a� � �m�1�.��� � �7�-�e�l�e�m�e�n�t� �V�e�c�t�o�r�{�F�l�o�a�t�6�4�}�:� � � �3�.�6�9�1�9�2�7�3�3�7�0�8�7�6�9�9�7� � � �2�.�8�6�1�2�1�1�2�4�1�1�1�7�6�1�4� � � �3�.�1�0�1�8�6�0�3�6�1�1�7�8�3�7�5� � � �0�.�3�5�5�0�8�6�3�0�3�7�2�5�4�0�1�7�3� � � �2�.�1�7�9�9�8�2�0�9�9�4�3�4�4�6�3�6� � � �2�.�3�1�6�4�0�9�3�4�5�0�6�5�3�4�3�6� � � �9�.�3�3�0�8�5�4�7�1�4�7�3�1�7�4�3� � �� � �