不使用statsmodels.api完成多因变量的多元线性回归可行吗，请给出简单案例和代码，请给出多个答案

是的，可以使用其他 Python 库以及自己编写代码来完成多因变量的多元线性回归。以下是一些可行的方案：

1. 使用 numpy 和 scipy 进行矩阵运算：

import numpy as np
from scipy.linalg import inv

# 构造数据
X = np.array([[1, 2, 3], [4, 5, 6], [7, 8, 9]])
Y = np.array([[10, 20], [30, 40], [50, 60]])

# 增加常数列
X = np.concatenate([np.ones((X.shape[0], 1)), X], axis=1)

# 计算回归系数
beta = np.dot(np.dot(inv(np.dot(X.T, X)), X.T), Y)

# 输出结果
print('Coefficients:', beta)

在这个示例中，我们使用 numpy 构造了一个形状为 (3,3) 的自变量矩阵 X，和一个形状为 (3,2) 的因变量矩阵 Y。然后，我们在 X 中增加一个常数列，以便拟合截距。接下来，我们使用矩阵运算求解回归系数，其中 inv 函数用于计算矩阵的逆。最后，我们输出了回归系数。

2. 使用 tensorflow 进行计算图构建和优化：

import tensorflow as tf

# 构造数据
X = np.array([[1, 2, 3], [4, 5, 6], [7, 8, 9]])
Y = np.array([[10, 20], [30, 40], [50, 60]])

# 增加常数列
X = np.concatenate([np.ones((X.shape[0], 1)), X], axis=1)

# 创建计算图
graph = tf.Graph()
with graph.as_default():
    # 定义变量和占位符
    w = tf.Variable(tf.zeros([X.shape[1], Y.shape[1]]))
    x = tf.placeholder(tf.float32, [None, X.shape[1]])
    y = tf.placeholder(tf.float32, [None, Y.shape[1]])

    # 构建模型
    y_pred = tf.matmul(x, w)

    # 定义损失函数和优化器
    loss = tf.reduce_mean(tf.square(y_pred - y))
    optimizer = tf.train.GradientDescentOptimizer(0.01).minimize(loss)

    # 训练模型
    with tf.Session() as sess:
        sess.run(tf.global_variables_initializer())
        for i in range(1000):
            _, l = sess.run([optimizer, loss], feed_dict={x: X, y: Y})
            if i % 100 == 0:
                print('Iteration:', i, ' Loss:', l)

        # 输出结果
        beta = sess.run(w)
        print('Coefficients:', beta)

在这个示例中，我们使用 tensorflow 构建了一个计算图来执行多元线性回归。首先，我们使用 numpy 构造了一个形状为 (3,3) 的自变量矩阵 X，和一个形状为 (3,2) 的因变量矩阵 Y。然后，我们在 X 中增加一个常数列，以便拟合截距。接下来，我们在计算图中定义了变量、占位符、模型、损失函数和优化器。最后，我们使用 tf.Session 进行训练，并输出了回归系数。

需要注意的是，以上示例仅仅是展示了多元线性回归的实现方式，对于真实数据集，还需要进行数据预处理、特征工程、模型选择、交叉验证等步骤，以获得更准确和鲁棒的结果。