BookmarkSubscribeRSS Feed
Frank_PFC
Calcite | Level 5
I have 68 2km x 2km plots for which the proportions of 5 land cover classes (proportions add to 1.0) have been determined by 2 methods. How can I determine that the two methods generate a significantly different result?

I can use a pair-wise t-test to determine if the proportions generated by the two methods in the same class are significantly different. What I am looking for is an overall assessment of whether or not the two methods generate a significantly different result across all 5 land cover classes.

Thanks in advance for any input.
4 REPLIES 4
StatDave
SAS Super FREQ
The table to be analyzed here is a 5x5 table (method1 x method2) and each of the 68 plots falls into one of the 25 cells. As stated, it sounds like you want a test of marginal homogeneity. That is, whether the 5 row proportions (marginal proportions for one method) are the same as the 5 column proportions (marginals for the other method.

You can test for marginal homogeneity as shown in the example titled "Repeated Measures, 4 Response Levels, 1 Population" in the PROC CATMOD documentation. You will need to create a data set with one observation per plot and with METHOD1 and METHOD2 variables containing the plots' classes under each method. Then METHOD1 and METHOD2 are used like the RIGHT and LEFT variables in the example. See this usage note for additional interpretation:

http://support.sas.com/kb/39243
Frank_PFC
Calcite | Level 5
Thank you for your suggestion.

The approach you suggest works for the case where the 2 classifications are intersected spatially in a geographic information system and the area common to each of the 5 land cover classes in the two methods is evaluated in a contingency table to calculate a user's and producer's accuracy and an overall accuracy.

In my case, I would like to test the similarity of the two sets of proportions as if they are describing two independent populations. In this case, there are two sets of proportions that add up to 1.0 in each plot. I want to know if the two sets of proportions are significantly different from one another.
Dale
Pyrite | Level 9
Frank,

You can use the MIXED procedure to test multiple correlated land cover class differences. First, set up your data as within-plot differences for the first four cover classes. (The fifth cover class is an exact linear combination of the first four, so it can be excluded from use. If you had just two cover classes, you would only use the within-plot difference of the first land cover class as your response variable.) You will need each within-plot difference as a separate result with a label indicating which land cover class the difference belongs to. That is, you want your data to be constructed as shown below:

      plot       CoverClass           diff    
        1                 1             diff[1,1]
        1                 2             diff[1,2]
        1                 3             diff[1,3]
        1                 4             diff[1,4]
        2                 1             diff[2,1]
        2                 2             diff[2,2]
        2                 3             diff[2,3]
        2                 4             diff[2,4]
       ...                ...                ...

With your data in this form, the following PROC MIXED code would be appropriate:

proc mixed data=mydata;
  class plot CoverClass;
  model diff = CoverClass / s noint;
  repeated CoverClass / subject=plot type=un;
run;
Frank_PFC
Calcite | Level 5
Thanks very much for your suggestion, Dale. I haven't had a chance yet to try it out but my colleagues think it should do the trick.

sas-innovate-2024.png

Join us for SAS Innovate April 16-19 at the Aria in Las Vegas. Bring the team and save big with our group pricing for a limited time only.

Pre-conference courses and tutorials are filling up fast and are always a sellout. Register today to reserve your seat.

 

Register now!

What is ANOVA?

ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.

Find more tutorials on the SAS Users YouTube channel.

Discussion stats
  • 4 replies
  • 1255 views
  • 0 likes
  • 3 in conversation