nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

Why does the Sum of Squares under "Effect Tests" not add up to the overall model's Sum of Squares?

Bashburz1

Hi all,

I fitted a simple model with one continuous and one nominal variable using Standard Least Squares. Now I am wondering why the sums of squares seem to mismatch under Analysis of Variance and Effect Tests. I am on JMP 18.0.1.

3.16 and 0.01 do not give the overall model's SS of 3.37. It's like some 0.2 are missing. Could anyone please explain why that is?

This particularly puzzles me because I am looking at analysis of the same data in minitab provided by colleagues, and there the factor "Month" indeed gets the increased SS of 0.209. All the other stuff like DF and SS of everything else matches. In the end, this has implications for the p-value to judge the significance of "month" for the fit, which is why this bothered me in the first place.

Regards,

Konstantin

Accepted answers

statman

It would really help to look over your results if you attach the actual data table instead of just the outputs. A little bit like "show your work" in school.

One observation I can instantly make in your comparison with Minitab is Minitab is reporting sequential SS (type 1) not partial SS (type 3). Sequential SS are order dependent, while partial SS are more robust to model order (and are preferred and usually the default). JMP will use type 3 as the default and also uses REML (I think Minitab is still on the EMS estimation platform).

I'm not sure what your factors are, and which are continuous and which are nominal. Neither batch nor month seem continuous? You have 9 batches and 2 months. but I can't tell where your DF for MSE is from?

If you look at the parameter estimates table and add the VIF column, what are those values?

All comments

statman

It would really help to look over your results if you attach the actual data table instead of just the outputs. A little bit like "show your work" in school.

If you look at the parameter estimates table and add the VIF column, what are those values?

Bashburz1

Thank you very much, it is much clearer already!

I have attached the data to this reply. The VIFs are fairly small.

Your observation is very enlightening to me. Yesterday I was actually wondering if the order of the model plays a role, and switched month and batch around, but JMP did not care. Perfectly makes sense with your explanation with Type 3. I now also understand that the effects in JMP are adjusted for the presence of other effects, while minitab does not do this, which explains the difference. After removing "batch" from the model, JMP gives an SS of 0.2 for "month", which is what minitab outputs from the model with both factors.

I believe this question is resolved.

However, I have also been wondering similarly as you about the DFs. Month is a continuous factor. So why does it end up with a DF of 1? Maybe you get an idea after looking at the data.

data.jmp

MRB3855

Hi @Bashburz1 : "However, I have also been wondering similarly as you about the DFs. Month is a continuous factor. So why does it end up with a DF of 1?"

That's because as a continuous factor, it only introduces one parameter to the model; the "slope". i.e., the term that is multiplied by Month.

Bashburz1

Got it, thank you!

statman

BTW, if you change the data type for month to nominal you will get 9 DF's for month. I'm not sure why you consider month a continuous variable?

Also if you click on the red triangle next to the response in the Fit Model output, select Estimates>Sequential Tests, you will get type 1 SS.

MRB3855

Hi @statman My educated guess is that “month” is indeed continuous; it is actually time since the sample was put on stability…I.e., in a chamber at a given temperature and relative humidity. But its units are months (could be seconds, minutes, days, weeks, etc). I’m guessing this is stability analysis (shelf-life) in the pharmaceutical industry. Am I correct @Bashburz1? if so, “month” is continuous.