crop.dataorg | R Documentation |
LDAcrop.pro
and LDAcrop.plus
.
The function transforms raw archaeobotanical data calculating the square root of the weed seeds percentage within each entered sample.
crop.dataorg(dataframe, codes, samples)
dataframe |
The dataframe with the archaeobotanical samples |
codes |
The column containing the grouping variables of the different species (BHH, BFH etc) |
samples |
The column in which the samples data starts |
The function conducts a square-root transformation of the weed seeds percentages as per Jones (1984).
Elizabeth Stroud
Jones, G. 1984. Interpretation of archaeological plant remains: ethnographic models from Greece. In (eds) W. van Zeist & W.A. Casparie, Plants and Ancient Man: Studies in Palaeoethnobotany. Rotterdam: Balkema, pp 43-61
Jones, G., 1987. A statistical approach to the archaeological identification of crop processing. Journal of Archaeological Science, 14(3), pp.311-323.
## Example data
species<-c("Chenopodium album" , "Anthemis cotula", "Brassica rapa ssp campestris",
"Raphanus raphanistrum", "Agrostemma githago" , "Poa annua" )
category<-c("SFH", "BHH","SFH","BHH","BFH", "SFL" )
s.1246<-sample(1:3, 6, replace=TRUE)
s.46178<-sample(1:5, 6, replace=TRUE)
s.1<-sample(0:8, 6, replace=TRUE)
s.23<-sample(0:3, 6, replace=TRUE)
s.987<-sample(3:9, 6, replace=TRUE)
dataset<-data.frame(species,category,s.1246,s.46178,s.1,s.23,s.987)
## Usage
data<-crop.dataorg(dataset, codes=2, samples=3)
data
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.