Adding a New Column/Array to a Numpy Array

I am trying to add a column to a numpy array. Each row currently has four values, I would like each row to have five values. Below is the reproducible example which returns ValueError: all the input arrays must have same number of dimensions I do not understand why I get the error because Y has the same length as X, just as b has the same length as a in the documentation. Ultimately, I would like the most efficient way to add an array like Y to an existing array like X as a new column for each row.

import numpy as np
from sklearn import datasets

#Documentation
a = np.array([[1, 2], [3, 4]])
b = np.array([[5, 6]])
c = np.concatenate((a, b.T), axis=1)
print(c)

#My Case
iris = datasets.load_iris()
X = iris.data
Y = iris.target
Z = np.concatenate((X, Y.T), axis = 1) #Is transpose necessary for single dimension array? Throws error either way
print(Z)

Edit: I should add that in practice what I will be working with is predicted values from a model fit from sklearn. So what I am particularly interested in is what the most efficient way to add the predicted values to an existing array like X, i.e. the format used by sklearn. This solution below is from M4rtini's comment, which I believe is equivalent to one of Dietrich's solutions. Is this the fastest implementation?

#My Case
import numpy as np
from sklearn import datasets
from sklearn.linear_model import LinearRegression

iris = datasets.load_iris()
X = iris.data
Y = iris.target
model = LinearRegression()
model.fit(X,Y)
y_hat = model.predict(X).reshape(-1,1)
Z = np.concatenate((X, y_hat), axis = 1)