Planning and analyzing a 2-level full factorial design in Python

🎓 Who I am:

I have a bachelors degree in coating science and a masters degree in material science. Currently I am doing my PhD in coatings technology. Also been passionate about teaching, I want to share my knowledge to people everywhere around the world.

[PierriDellaBruna]

Great video! Really valuable material. Congratulations :)

[fightchris]

Nice, explained very well. I like the approach to let ChatGPT write most of the code. Keep going!

[Imran_et_al]

Great Explanation. Thanks for the video and code...!

[sebastianoverri]

Love your work. Thanks, your teaching is very clear even for a newby inthe subject like me.

[christophersaid]

This is great, thanks for posting. I work at a startup that's extremely cheap, so having DOE capabilities in python is very cool. Would you be able to share the actual Jupyter Notebook? I couldn't find it on your site.

[victorrossing]

I have a problem when doing DOE in python. Lets say we have a 2^2 experiment single replicate with 5 centerpoints, then we have 4 degrees of freedom. One each of x1, x2, x1x2, whereby we only have one combined for x1^2 and x2^2, hence we need a pure quadratic term for these two combined, as they cant be separated. However OLS cant give me the pure quadratic term (just gives me an estimate of x1^2, which is equal to the estimate of x2^2), do you have a solution or another package you know that works? Code attached (chatGPT cant find a solution):

test=pd.DataFrame({})
test["x1"]=[-1, -1, 1, 1, 0, 0, 0, 0, 0]
test["x2"]=[-1, 1, -1, 1, 0, 0, 0, 0, 0]
test["response"]=[39.3, 40, 40.9, 41.5, 40.3, 40.5, 40.7, 40.2, 40.6]

model_poly = ols('Q("response") ~ x1 + x2 + I(x1 * x2) + I(x1**2 +x2**2)', data=test).fit()
print(model_poly.summary())

[rusty-yr]

pydoe2 is not not importing well in anaconda python.. it says imp is needed but i cannot install imp in newer version of python. Can someone share what python version works for pydoe2?

[BarisDuygu-qu]

Hey Marcel, Ive got the issue, that the results I have had some slip region included. Now of course when I remove these points, my cubic won't be complete. I tried to make a model based on the remaining runs, Used anova for this.

[mohammedbenali]

Thank you for this video.
Please, if possible how we can do doehlert design in python,
Video the scientific community please.