i need to convert the following code to python in this pmax.pcomp and pumpmark are columns in the dataframe STD

// do the following calculations on different fields to find needed values for chart -------------------

// find linear model b/w fileds pmax.pcomp & orthogonal polynomials of pumpmark with degree 2 --------

fit2aa <- lm(as.numeric(STD$pmax.pcomp) ~ poly(as.numeric(STD$pumpmark), 2, raw=TRUE))