from sklearn.datasets import load_iris
from sklearn.preprocessing import PolynomialFeatures

# load the iris dataset
iris = load_iris()
X = iris.data

# create PolynomialFeatures object
poly = PolynomialFeatures(interaction_only=True)

# fit and transform the data
X_poly = poly.fit_transform(X)