Become an expert in R — Interactive courses, Cheat Sheets, certificates and more!
Get Started for Free

murdoch

Indicator Species Analysis by Murdoch Preference Function


Description

Calculates the indicator value of species in a single cluster or environment type using the Murdoch Preference Function

Usage

murdoch(comm,type,minval=0,minplt=10)
## S3 method for class 'murdoch'
summary(object,pval=0.05,digits=3,...)
## S3 method for class 'murdoch'
plot(x,axtype=1,pval=0.05,...)
## S3 method for class 'murdoch'
print(x,digits = 5, ...)

Arguments

comm

a matrix or data.frame of samples with species as columns and samples as rows

type

a logical vector with values of TRUE for samples in a specific cluster or type

minval

a threshold minimum abundance value to count as a presence

minplt

the minimum number of presences to include a species in the calculation

object

and object of class ‘murdoch’

pval

the maximum probability to include a species in the summary table

digits

the number of digits to report

...

ancillary arguments to maintain compatibility with the generic summary function

x

an object of class ‘murdoch’

axtype

a switch to control scaling of the x axis in the plot. 1=number of plots in the data set, other = number of presences in the type

Details

Calculates the indicator value of species for a specific type using the modified Murdoch statistic:

log((p/a) * (n-p_i)/n_i)

where: p = number of samples where species is present, a = number of samples where species is absent, n = total number of samples (p+a), p_i = number of samples in type i where species is present, n_i = number of samples in type i.

Probabilities are based on the hypergeometric distribution calculation of having as many or more presences in a type as observed.

Value

a list object of class ‘murdoch’ with components:

minplt

the minimum number of presences to be included

nplots

the number of plots a species occurs in

type

the plot membership vector for the type

pres

the number of presences for species in the type

abs

the number of absences of species in the type

murdoch

the Murdoch value for species in the type

pval

the probability of getting such a high murdoch value

Note

Indicator value analysis is a set of techniques designed to identify species of special interest in clusters or types. The most widely used indicator species analysis was proposed by Dufrene and Legendre (1997), and is included in package ‘labdsv’ as indval. murdoch differs significantly from indval in assumption and objective, seeking to identify species that have improbable occurrences in types, regardless of their relative frequency in the type

Author(s)

David W. Roberts with help from Ken Aho droberts@montana.edu http://ecology.msu.montana.edu/labdsv/R

See Also

Examples

data(shoshveg)                           # returns a vegetation dataframe
dis.bc <- dsvdis(shoshveg,'bray/curtis') # returns a dissimilarity 
                                         #   matrix
opt.5 <- optpart(5,dis.bc)
plot(murdoch(shoshveg,opt.5$clustering==1))

optpart

Optimal Partitioning of Similarity Relations

v3.0-3
GPL (>= 2)
Authors
David W. Roberts <droberts@montana.edu>
Initial release

We don't support your browser anymore

Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.