我想拟合一个好几条线性关系一次性

我想用LsqFit.jl来拟合多组线性关系,但是现在一组线性关系都成问题,就是我有64个自变量,32个应变量;每个应变量都和64个自变量有线性关系。神经网络应该也是可以做得到,但是没时间调参数了我怕调的不好。其实我也不知道怎么会更快捷一些。

希望大家可以提一些意见,然后帮我看一下我哪里错了。

function multimodel(x, p)
    b=0
    n=5
    for i=1:n
        b=b .+ @.x[:,i]*p[i]
    end
    b=b.+p[n+1]
end
x= [1 2 3 4 5;1 2 3 4 5;1 2 3 4 5;1 2 3 4 5;]
p= [1 2 3 4 5 6]
y= [1 1 1 1]
multimodel(x,p)
fit = curve_fit(multimodel, x, y, p)