@ckrapu Hi Chris, sorry to bother you again. The code I have below uses for loop to dot product of every two columns with beta vector, totally 8 columns in the matrix. I wonder whether this is better way to do this. This is inside the with part. Thanks a lot.
# Define utility for fixed effects
utility_fixed = tt.dot(X_fixed[:,[0,1]], beta)
utility_fixed = tt.reshape(utility_fixed, (69,1))
for i in range(3):
part = X_fixed[:, [2 + i, 3 + i]]
dot_prod = tt.dot(part, beta)
dot_prod = tt.reshape(dot_prod, (69,1))
utility_fixed = tt.concatenate([utility_fixed, dot_prod], axis = 1)