Skip to main content
Commonmark migration
Source Link

###This was a part of the original question, but it can be misleading, see the answer below.

This was a part of the original question, but it can be misleading, see the answer below.

###This was a part of the original question, but it can be misleading, see the answer below.

This was a part of the original question, but it can be misleading, see the answer below.

deleted 4 characters in body
Source Link
James S.
  • 451
  • 1
  • 4
  • 9

First herehere is my data. I had to link to it because of its size. Note that the biomass column has been standardized (mean = 0, sd = 1), hence the negative values. This does not alter the output. Once downloaded and the working directory has been specified, the file can be read in as follows:

First here is my data. I had to link to it because of its size. Note that the biomass column has been standardized (mean = 0, sd = 1), hence the negative values. This does not alter the output. Once downloaded and the working directory has been specified, the file can be read in as follows:

First here is my data. I had to link to it because of its size. Note that the biomass column has been standardized (mean = 0, sd = 1), hence the negative values. This does not alter the output. Once downloaded and the working directory has been specified, the file can be read in as follows:

edited title
Source Link
amoeba
  • 107.3k
  • 36
  • 321
  • 347

The problem

These results clearly differ. The B:C interaction is no longer significant and the P-value for the A:B interaction is quite a bit higher. Both models should be computing the P-values in similar ways and so it's hard to imagine them being so different.

Why are they different?


###This was a part of the original question, but it can be misleading, see the answer below.

In fact, it seems that the anova(model, type = 3) function is actually using type 2 SS, which we can verify by running anova(model, type = 2).

 

The problem

These results clearly differ. The B:C interaction is no longer significant and the P-value for the A:B interaction is quite a bit higher. Both models should be computing the P-values in similar ways and so it's hard to imagine them being so different. In fact, it seems that the anova(model, type = 3) function is actually using type 2 SS, which we can verify by running anova(model, type = 2).

These results clearly differ. The B:C interaction is no longer significant and the P-value for the A:B interaction is quite a bit higher. Both models should be computing the P-values in similar ways and so it's hard to imagine them being so different.

Why are they different?


###This was a part of the original question, but it can be misleading, see the answer below.

In fact, it seems that the anova(model, type = 3) function is actually using type 2 SS, which we can verify by running anova(model, type = 2).

 
edited title
Link
amoeba
  • 107.3k
  • 36
  • 321
  • 347
Loading
edited tags; edited title
Link
amoeba
  • 107.3k
  • 36
  • 321
  • 347
Loading
replaced http://stats.stackexchange.com/ with https://stats.stackexchange.com/
Source Link
Loading
deleted 617 characters in body
Source Link
James S.
  • 451
  • 1
  • 4
  • 9
Loading
Added answer to problem
Source Link
James S.
  • 451
  • 1
  • 4
  • 9
Loading
Post Reopened by amoeba, gung - Reinstate Monica, Jake Westfall, Matt Krause, mdewey
added 1052 characters in body
Source Link
James S.
  • 451
  • 1
  • 4
  • 9
Loading
Added data (linked) for reproducibility and clarified the language in the last paragraph
Source Link
James S.
  • 451
  • 1
  • 4
  • 9
Loading
added 8 characters in body
Source Link
James S.
  • 451
  • 1
  • 4
  • 9
Loading
Source Link
James S.
  • 451
  • 1
  • 4
  • 9
Loading