Before in general moving onto splines and other smoothing methods and applying more transforms, I wanted to continue looking at how a polynomial representation of the model would affect the fit.
Since the last time I was comparing fits just based off a polynomial for obesity, this will compare both, using just inactivity, and using inactivity and obesity together.

We see here on the basis of the P values that we actually only need up to a quadratic function.
Next I tried combining polynomials for both inactivity and obesity to see if it affected the fit in a meaningful way.

As seen here by the R square values we can tell that it was not a significant increase.
After this I tried taking log functions of the polynomials to see if that would alter my results significantly.

As we can see from the R squared values again this was not producing anything drastically different in terms of what the fit was.
To further modify things I have tried to add an interaction term and as well as taking a log of diabetes as well.

As we can see from our latest test that the addition of the interaction term changed things rather significantly as compared to the previous transforms and this also helped increase the R squared value.