Hi @StatlearnR,
I'm not entirely sure about the design you have used (a datatable would have been easier to figure out), but you may find useful answers about similar questions here :
How to analyze a 4x4 Latin square design?
Analyze Latin squares using JMP
Repeated randomized 4 by 4 latin square
It looks like you can analyze your data by specifying Treatment as a fixed effect in your model, and rows and columns of your LS design as random blocking effects, using the Fit Model platform.
I did a simple dataset and demo trying to fit to your infos and context :
Datatable is attached with the script "Model" to launch the correct modeling, and "Fit Least Squares Yield" with the result of the modeling. Random Effects are however not supported by ordinal logistic models, so you might have to transform your response, or use another one. I added as an example a binning formula column of Yield with ordinal data, you can try launching ordinal logistic regression but it won't work with random blocking effects.
Hope this may help you. If not, please bring more information and if possible, an (anonymized) datatable.
Victor GUILLER
"It is not unusual for a well-designed experiment to analyze itself" (Box, Hunter and Hunter)