Description Usage Arguments Details Value Author(s) References Examples

Calculates a set of values for a particular statistic or sets of observations, typically for observed values and multiple sets of randomized observations.

1 |

`rand` |
A list for which the elements are either sets of randomized variables or randomized index values. |

`marks` |
Should either be left empty, be a vector or a matrix. When |

`FUN` |
Any function used to calculate the statistic of interest (e.g. |

`...` |
Optional arguments to |

`cal.stat`

is designed to calculate statistics for spatial explicit data for which randomized data sets are generated with `fgrand`

.

`cal.stat`

returns a vector of statistics. If `rand`

is the output of `fgrand`

and `add.obs=TRUE`

, the first value is the statistic for the observed data and the rest for randomizations.

Reinder Radersma

Reinder Radersma & Ben C. Sheldon, 2014. A new permutation test for dealing with and exploring spatial autocorrelation submitted to MEE

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 | ```
#### Example for fgrand
## 200 random geographical locations
xy <- array(runif(400,0,2), dim=c(200,2))
## run fgrand to produce 99 randomizations for scale 1
test <- fgrand(xy, scale=1, iter=99, add.obs=TRUE)
## run fgrand to produce 99 bootstraps for scale 1
test <- fgrand(xy, scale=1, iter=99, FUN=function(x){x[sample.int(length(x),replace=TRUE)]}, add.obs=TRUE)
## 200 times 200 random distances (e.g. genetic relatedness between mated pairs)
trait <- array(rnorm(200*200,0.6,0.1), dim=c(200,200))
## make the observed pairs more alike
diag(trait) <- diag(trait)+0.02
## make two rows and two colums empty
trait[,3] <- NA
trait[,50] <- NA
trait[6,] <- NA
trait[12,] <- NA
## calculate means; will give NAs because there are missing values
calc <- cal.stat(test,trait,mean)
## calculate means
calc <- cal.stat(test,trait,mean, na.rm=TRUE)
## plot means
hist(calc)
abline(v=calc[1], col="red", lwd=2)
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.