View source: R/get_sampling_periods.R
get_sampling_periods | R Documentation |
Converts several different types of data storage into sampling periods for calculating or permuting networks
get_sampling_periods(association_data, association_times, sampling_period,
identities = NULL, location = NULL, within_locations = FALSE,
data_format = c("gbi", "groups", "individuals"), return="SP")
association_data |
Can be either a group by individual matrix, a list containing group members in each element, or a two-column data frame with individual ID in the first column and group ID in the second column |
association_times |
Because sampling periods are inferred over time, each group must contain some time data (can be in any format, such as seconds, days, etc.). One time must be provided for each row of the association data. |
sampling_period |
The number of time periods over which data are combined (for example 10 days, 3600 seconds) |
identities |
Optional identities for each individual in the dataset |
location |
If spatial disaggregation need to be maintained, samping periods can be calculated per time per location |
within_locations |
Flag whether to include location information |
data_format |
Format of the input data |
return |
By default ("SP") returns the sampling periods. Anything else will return the occurrence data (see get_network function) |
This function will calculate an association matrix for each sampling period. If locations are included, these will be treated independently.
Returns a K x N x N
stack of matrices, where each N x N
slice is an association matrix. Row names and Column names of these slices are given the identity where available. The K slice names are given either the time or time_location for each sampling period. Alternatively (return != "SP") the function returns the occurrence of each individual in each sampling period, with individuals as rows and sampling periods as columns.
Damien R. Farine
## define group memberships (these would be read from a file)
individuals <- data.frame(ID=c("C695905","H300253","H300253",
"H300283","H839876","F464557","H300296","H300253",
"F464557","H300296","C695905","H300283","H839876"),
GROUP=c(1,1,2,2,2,3,3,4,5,5,6,6,6))
## create a time column
individuals <- cbind(individuals,
DAY=c(1,1,1,1,1,2,2,2,3,3,3,3,3))
SPs <- get_sampling_periods(individuals[,c(1,2)],
individuals[,3],1,data_format="individuals")
occurs <- get_sampling_periods(individuals[,c(1,2)],
individuals[,3],1,data_format="individuals", return="occ")
## define group memberships (these would be read from a file)
groups <- list(G1=c("C695905","H300253"),
G2=c("H300253","H300283","H839876"),
G3=c("F464557","H300296"),
G4=c("H300253"),
G5=c("F464557","H300296"),
G6=c("C695905","H300283","H839876"))
## create a time variable
days <- c(1,1,2,2,3,3)
SPs <- get_sampling_periods(groups,
days,1,data_format="groups")
occurs <- get_sampling_periods(groups,
days,1,data_format="groups", return="occ")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.