Skip to content
Prev 306199 / 398506 Next

Script to count unique values from two linked matricies

Hi,

Thanks for helping me with this one.

To save you time, the following is the code for the tables I uploaded as jpegs that you may not have received:

dat1<-data.frame(Species=paste("Species",1:3),Site1=c(5,NA,4),Site2=c(3,2,4),Site3=c(NA,5,NA))
dat2<-data.frame(Species=paste("Species",1:3),TraitType1=c(1,2,3),TraitType2=c(2,4,2),TraitType3=c(5,1,1))
dat3<-data.frame(Site=paste("Site",1:3),Trait_Richness=c(5,7,3))

So,

I have two matricies, dat1 and dat2.

Dat1 is a species abundance matrix. Dat2 is a species trait matrix.

I want to create dat3 through use of a script.

Dat 3 is a count of unique traits observed at each site. i.e. at site 1, species 1 and 3 are present (ass seen in dat1). Species 1 has traits: 1, 2 and 5 for trait types 1, 2 and 3 respectively. Species 3 has traits: 3, 2 and 1 for trait types 1, 2 and 3 respectively.

So, at site 1:

For trait type 1, 2 unique traits were observed. For trait type 2, 1 unique trait was observed (both species 1 and 3 were classed as "1") and for trait type 3, 2 unique traits (trait richness) were observed; thus, 2+1+2=5.

and so on... so at site 2, all three species were observed...

For trait type 1, 3 unique traits were observed (1, 2, 3), for trait type 2, 2 unique traits were observed (2, 4, 2) and for trait type 3, 2 unique traits were observed (5, 1, 1). So, for site 2, trait richness is 7 (3+2+2) traits.

I hope this helps to explain, please let me know if you need any further information,

Ben Gillespie
Research Postgraduate

School of Geography
University of Leeds
Leeds
LS2 9JT

Tel: +44(0)113 34 33345
Mob: +44(0)770 868 7641
http://www.geog.leeds.ac.uk/