Python 3.x 使用数据集的Python索引和边界错误

我们的类使用Python作为模型的解决方案工具。然而，这是我自1997年VB以来第一次使用python或任何编程语言，所以我很挣扎。我们已经向我们提供了以下代码

from numpy import loadtxt, array, ones, column_stack
from numpy import dot, sqrt
from scipy.linalg import inv
from scipy.stats import norm, t


f = loadtxt('text data.raw')

y = f[:,4]
n = y.size
x = array([f[:,2],f[:,8],f[:,4]])
one = ones(n)
#xa = column_stack([one,f[:,3],f[:,4]])
xa = column_stack([one,x.T])
k = xa.shape[1] 

xx = dot(xa.T,xa)
invx = inv(xx)
xy = dot(xa.T,y)
b = dot(invx,xy)

# Compute cov(b)
e = y - dot(xa,b)
s2 = dot(e.T,e)/(n-k)
covb = invx*s2

# Compute t-stat
tstat = b[1]/sqrt(covb[1][1])

#compute p-value
p = 1 - norm.cdf(tstat,0,1)
pt = 1 - t.cdf(tstat,88)

我们的数据集是一个10x88矩阵。我们的目标是创建一个线性规划并找到一些答案。在我们的数据中，第1列已经设置为价格，这在我们的线性规划中是我们的期望产量，我需要使用第3、4和5列。作为我的x1、x2和x3。我不确定为了完成这项任务，第9行和第11行的值需要更改为什么或如何更改，目前我也不清楚这两行在程序中具体要求或执行什么。再说一遍，我不熟悉编程

我尝试的每件事通常都会产生类似的错误

IndexError: index 5 is out of bounds for axis 1 with size 5

---

有什么建议吗

请提供fullstack错误。请提供fullstack错误。