Creates model latent variables

        Returns
        ----------
        None (changes model attributes)

def _create_latent_variables(self):
        """ Creates model latent variables

        Returns
        ----------
        None (changes model attributes)
        """

        self.latent_variables.add_z('Vol Constant', fam.Normal(0,3,transform=None), fam.Normal(0,3))

        for component in range(2):
            increment = 0.05
            for p_term in range(self.p):
                self.latent_variables.add_z("Component " + str(component+1) + ' p(' + str(p_term+1) + ')', fam.Normal(0,0.5,transform='logit'), fam.Normal(0,3))
                if p_term == 0:
                    self.latent_variables.z_list[1+p_term+component*(self.p+self.q)].start = 3.00
                else:
                    self.latent_variables.z_list[1+p_term+component*(self.p+self.q)].start = 2.00

            for q_term in range(self.q):
                self.latent_variables.add_z("Component " + str(component+1) + ' q(' + str(q_term+1) + ')', fam.Normal(0,0.5,transform='logit'), fam.Normal(0,3))
                if p_term == 0 and component == 0:
                    self.latent_variables.z_list[1+self.p+q_term+component*(self.p+self.q)].start = -4.00  
                elif p_term == 0 and component == 1:
                    self.latent_variables.z_list[1+self.p+q_term+component*(self.p+self.q)].start = -3.00  

        self.latent_variables.add_z('v', fam.Flat(transform='exp'), fam.Normal(0,3))
        self.latent_variables.add_z('Returns Constant', fam.Normal(0, 3, transform=None), fam.Normal(0,3))
        self.latent_variables.z_list[-2].start = 2.0