wdata_org | R Documentation |
wmodel.LDA
Changes archaeobotanical data or other cases the investigator wishes to classify, such as survey data from other farming regimes, into presence/absence data and extracts the functional data for the listed species, and calculates the average of each functional attribute for each sample.
wdata_org(dataframe, samples, codes,codename, model, fl_pr=NULL, sp_av=NULL, vg_pr=NULL)
dataframe |
The data set to be transformed. Data should be arranged with each sample as a column and each species as rows. One column needs to contain the species codes of the taxa (see LOCATION for list of taxa and corresponding codes). |
samples |
The column number in which the sample data starts |
codes |
The column number containing the species codes |
codename |
The name of the column containing the species codes |
model |
The model being used; 1, 2 or 3 - this will determine which functional attributes are used |
fl_pr |
The flowering period data for the species (needs to be set up as two columns: one containing species codes and one containing flowering data labelled FLOWPER). This is only required for model 1 and model 3 |
sp_av |
The averaged species data of composite taxa - this can be the ouput of |
vg_pr |
Modified vegetative proprogation values - all values VEGPROP values need to be re-entered if this argument is used |
A data frame containing the averaged functional attributes values for each sample. Which attributes are returned will be dependant on the model entered:
model 1: |
SLA, ARNODE, LOGCANH, LOGCAND, FLOWPER |
model 2: |
SLA, ARNODE, LOGCANH, LOGCAND |
model 3: |
FLOWPER, VEGPROP |
Elizabeth Stroud
PAPER REFERENCE
## Random data
species<-c("Chenopodium album" , "Anthemis cotula", "Brassica rapa ssp campestris",
"Raphanus raphanistrum", "Agrostemma githago" , "Poa annua" )
code<-c("chenalb", "anthcot", "brasrap","raphrap","agrogit", "poa_ann")
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,code,s.1246,s.46178,s.1,s.23,s.987)
code<-c("chenalb", "anthcot", "brasrap","raphrap","agrogit", "poa_ann")
FLOWPER<-sample(3:9, 6, replace=TRUE)
x<-data.frame(code,FLOWPER)
## Use
results<-wdata_org(dataset, samples=3, codes=2, codename="code", model=1, fl_pr=x)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.