From: Pedro de Barros <pbarros_at_ualg.pt>

Date: Wed 20 Jul 2005 - 22:32:36 EST

R-help@stat.math.ethz.ch mailing list

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Wed Jul 20 22:41:03 2005

Date: Wed 20 Jul 2005 - 22:32:36 EST

Dear All,

I am trying to build a randomization test for interaction

The problem is as follows: I have a set of stations where the ocurrence and biomass of each species being investigated was recorded.

For each group of stations, the relative importance of a set of species was calculated, as the Index of Relative Importance (IRI), given as

IRI(i)=(FO(i)*RW(i))/sum(FO(j)*RW(j))

FO(i) is the relative frequency of ocurrence of species i and RW(i) is the relative Biomass of species i (Total biomass of this species/Total biomass of all species). Stations are grouped according to 3 different criteria, say A, B and C, with 3 levels of A, 5 of B, and 3 of C.

The null hypothesis of interest is that the distribution of the IRI among the species is independent of these classification levels. I suggest using the average difference in IRI(i) of each grouping to the overall IRI(i) as the test statistic.

I have no difficulty building the randomization test for the main effects: Simply randomize the corresponding group membership factor. However, I am having difficulty building the test for the interaction effect: As far as I understand, the test should randomize the levels of say B, within the levels of A, so that the main effects should not be affected (all observations classified as level A1 would still be A1, and all observations classified as level B1 should stay as B1). My difficulty is how to calculate the "overall" distribution of the IRI in this case. This should consider all main effects, but no interaction.

I would really appreciate any pointer to a solution of this problem. I believe it is not complicated (and probably quite obvious) but the solution keeps out of reach, even though I have been searching for over a week.

Thanks,

Pedro

Pedro de Barros

Faculdade de Ciências do Mar e do Ambiente
Universidade do Algarve ,Campus de Gambelas
8000-117 FARO, Portugal

Tel: + 351 289 800 034

Fax: +351 289 818 353

e-mail: pbarros@ualg.pt

R-help@stat.math.ethz.ch mailing list

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Wed Jul 20 22:41:03 2005

*
This archive was generated by hypermail 2.1.8
: Fri 03 Mar 2006 - 03:33:51 EST
*