def METHOD_NAME(self):
    demo_estimator = SEMEstimator(self.demo)
    for method in ["random", "std"]:
        B_init, zeta_init = demo_estimator.get_init_values(
            data=self.demo_data, method=method
        )
        demo_lisrel = self.demo.to_lisrel()
        m = len(demo_lisrel.eta)
        self.assertEqual(B_init.shape, (m, m))
        self.assertEqual(zeta_init.shape, (m, m))
        union_estimator = SEMEstimator(self.union)
        B_init, zeta_init = union_estimator.get_init_values(
            data=self.union_data, method=method
        )
        union_lisrel = self.union.to_lisrel()
        m = len(union_lisrel.eta)
        self.assertEqual(B_init.shape, (m, m))
        self.assertEqual(zeta_init.shape, (m, m))