MFggplot: ggplot2 extension for a MF1_single or MF2_multiple object
In MF.beta4: Measuring Ecosystem Multi-Functionality and Its Decomposition
MFggplot R Documentation
ggplot2 extension for a MF1_single or MF2_multiple object
Description
MFggplot provides graphical BEF relationships based on the output from the function MF1_single or MF2_multiple.
Usage
MFggplot(output, model = "LMM.both", caption = "slope", by_group = NULL)
Arguments
output
the output obtained from MF1_single or MF2_multiple.
For output obtained from MF1_single, if BEF relationships are desired within each category specified by_group, the by_group column must be included in the input.
model
specifying the fitting model, model = "lm" for linear model; model = "LMM.intercept",
"LMM.slope" and "LMM.both" for linear mixed models with random effects for intercepts, slopes, and both, respectively.
Default is model = "LMM.both".
If output is obtained from MF2_multiple with by_pair = FALSE, model can only set as "lm".
caption
caption that will be shown in the BEF plots; caption = "slope" to show the estimated slopes in each plot,
or caption = "R.squared" to show the ordinary R-squared for linear models or estimated marginal and conditional R-squared for linear mixed models in each plot.
Default is caption = "slope".
by_group
the column name of the stratifying variable that is used to group data for model fitting. For example, if by_group = "country", then model will be fitted within each country. Default is NULL.
It is required if a linear mixed model is selected in the model.
If output is obtained from MF2_multiple with by_pair = TRUE, the by_group setting must be the same as that set in MF2_multiple.
If output is obtained from MF2_multiple with by_pair = FALSE, the by_group must be NULL, because alpha/beta/gamma data points are not sufficient to fit model in each group.
Value
For an MF1_single object, this function returns a figure that plots the BEF relationship between multifunctionality of
order q (= 0, 1 and 2) and species diversity of the same order q for two cases (i) correlations between functions are not corrected for, and (ii) correlations between
functions are corrected. The fitted lines for the chosen model are also shown in the figure.
For an MF2_multiple object of given individual function weights, this function returns a list of two objects ($corr_uncorrected and $corr_corrected)
that respectively for two cases: (i) correlations between functions are not corrected for, and (ii) correlations between functions are corrected for.
Each object consists of four figures: "$ALL" returns a figure that depicts the BEF relationship between alpha/beta/gamma multifunctionality of
order q (= 0, 1 and 2) and the corresponding species diversity of the same order q. The fitted lines for the chosen model are also shown in the figure.
"$Gamma" returns only the gamma part of "$ALL", "$Alpha" returns only the alpha part of "$ALL", and "$Beta" returns only the beta part of "$ALL".
Examples
library(dplyr)
### Use the entire of 209 plots data from six countries
## within-ecosystem MF computation
data("forest_function_data_normalized")
data("forest_biodiversity_data")
output1 <- MF1_single(func_data = forest_function_data_normalized[,6:31], weight = 1,
species_data = forest_biodiversity_data)
## BEF relationships: display fitted line of linear mixed model with random slopes and
## random intercepts
output1 <- data.frame(output1, country=rep(forest_function_data_normalized$country, each = 6))
MFggplot(output1, model = "LMM.both", by_group="country", caption = "slope")
### Use partial data to quickly obtain output
### (Take the first 18 plots in Germany and the last 18 plots in Italy)
## within-ecosystem MF computation for partial data
data("forest_function_data_raw")
data("forest_biodiversity_data")
GER_ITA_forest_function_raw <- filter(forest_function_data_raw,
country=="GER"|country=="ITA")[c(1:18,57:74),]
GER_ITA_forest_function_normalized <- function_normalization(data = GER_ITA_forest_function_raw,
fun_cols = 6:31,
negative = c("soil_cn_ff_10","wue"),
by_group = "country")
GER_ITA_forest_biodiversity <- forest_biodiversity_data[c(49:82,181:229),]
output2 <- MF1_single(func_data = GER_ITA_forest_function_normalized[,6:31], weight = 1,
species_data = GER_ITA_forest_biodiversity)
## BEF relationships: display fitted line of linear mixed model with random slopes and
## random intercepts
output2 <- data.frame(output2, country=rep(GER_ITA_forest_function_normalized$country, each = 6))
MFggplot(output2, model = "LMM.both", by_group="country", caption = "slope")
### Use data from plots in five countries (data in Finland are excluded)
## MF decomposition for all pairs of ecosystems
data("forest_function_data_normalized")
data("forest_biodiversity_data")
forest_function_data_normalized <- filter(forest_function_data_normalized, country != "FIN")
forest_biodiversity_data <- forest_biodiversity_data[-(1:48),]
output3 <- MF2_multiple(func_data = forest_function_data_normalized[,6:32],
species_data = forest_biodiversity_data,
weight = 1,
by_group = "country")
## BEF relationships: display fitted line of linear mixed model with random slopes and
## random intercepts
figure_LMM <- MFggplot(output3, model = "LMM.both", by_group = "country",
caption = "slope")
figure_LMM$corr_uncorrected$ALL
figure_LMM$corr_corrected$ALL
### Use partial data to quickly obtain output
### (Take the first 18 plots in Germany and the last 18 plots in Italy)
## MF decomposition for all pairs of ecosystems for partial data
data("forest_function_data_raw")
data("forest_biodiversity_data")
GER_ITA_forest_function_raw <- filter(forest_function_data_raw,
country=="GER"|country=="ITA")[c(1:18,57:74),]
GER_ITA_forest_function_normalized <- function_normalization(data = GER_ITA_forest_function_raw,
fun_cols = 6:31,
negative = c("soil_cn_ff_10","wue"),
by_group = "country")
GER_ITA_forest_biodiversity <- forest_biodiversity_data[c(49:82,181:229),]
output4 <- MF2_multiple(func_data = GER_ITA_forest_function_normalized[,6:32],
species_data = GER_ITA_forest_biodiversity,
weight = 1,
by_group = "country")
## BEF relationships: display fitted line of linear mixed model with random slopes and
## random intercepts for partial data
figure_LMM_GER_ITA <- MFggplot(output4, model = "LMM.both", by_group = "country",
caption = "slope")
figure_LMM_GER_ITA$corr_uncorrected$ALL
figure_LMM_GER_ITA$corr_corrected$ALL
### Use partial data to calculate multifunctionality based on 3 plots in each country, not by pairs
### (Take the first 3 plots in each country)
data("forest_function_data_raw")
data("forest_biodiversity_data")
forest_function_data_raw_3plots <- forest_function_data_raw[c(1:3,29:31,67:69,103:105,
146:148,174:176),]
forest_function_data_normalized_3plots <-
function_normalization(data = forest_function_data_raw_3plots,
fun_cols = 6:31,
negative = c("soil_cn_ff_10","wue"),
by_group = "country")
forest_biodiversity_data_3plots <-
forest_biodiversity_data[c(1:6,49:52,141:148,230:232,351:355,411:417),]
output5 <- MF2_multiple(func_data = forest_function_data_normalized_3plots[,6:32],
species_data = forest_biodiversity_data_3plots,
weight = 1,
by_group = "country", by_pair = FALSE)
## BEF relationships: display fitted line of linear model
figure_all_plots <- MFggplot(output5, model = "lm", caption = "slope")
figure_all_plots$corr_uncorrected$ALL
figure_all_plots$corr_corrected$ALL
MF.beta4 documentation built on Nov. 3, 2024, 5:06 p.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
| MFggplot | R Documentation |
ggplot2 extension for a MF1_single or MF2_multiple object
Description
MFggplot provides graphical BEF relationships based on the output from the function MF1_single or MF2_multiple.
Usage
MFggplot(output, model = "LMM.both", caption = "slope", by_group = NULL)
Arguments
output |
the output obtained from |
model |
specifying the fitting model, |
caption |
caption that will be shown in the BEF plots; |
by_group |
the column name of the stratifying variable that is used to group data for model fitting. For example, if |
Value
For an MF1_single object, this function returns a figure that plots the BEF relationship between multifunctionality of
order q (= 0, 1 and 2) and species diversity of the same order q for two cases (i) correlations between functions are not corrected for, and (ii) correlations between
functions are corrected. The fitted lines for the chosen model are also shown in the figure.
For an MF2_multiple object of given individual function weights, this function returns a list of two objects ($corr_uncorrected and $corr_corrected)
that respectively for two cases: (i) correlations between functions are not corrected for, and (ii) correlations between functions are corrected for.
Each object consists of four figures: "$ALL" returns a figure that depicts the BEF relationship between alpha/beta/gamma multifunctionality of
order q (= 0, 1 and 2) and the corresponding species diversity of the same order q. The fitted lines for the chosen model are also shown in the figure.
"$Gamma" returns only the gamma part of "$ALL", "$Alpha" returns only the alpha part of "$ALL", and "$Beta" returns only the beta part of "$ALL".
Examples
library(dplyr)
### Use the entire of 209 plots data from six countries
## within-ecosystem MF computation
data("forest_function_data_normalized")
data("forest_biodiversity_data")
output1 <- MF1_single(func_data = forest_function_data_normalized[,6:31], weight = 1,
species_data = forest_biodiversity_data)
## BEF relationships: display fitted line of linear mixed model with random slopes and
## random intercepts
output1 <- data.frame(output1, country=rep(forest_function_data_normalized$country, each = 6))
MFggplot(output1, model = "LMM.both", by_group="country", caption = "slope")
### Use partial data to quickly obtain output
### (Take the first 18 plots in Germany and the last 18 plots in Italy)
## within-ecosystem MF computation for partial data
data("forest_function_data_raw")
data("forest_biodiversity_data")
GER_ITA_forest_function_raw <- filter(forest_function_data_raw,
country=="GER"|country=="ITA")[c(1:18,57:74),]
GER_ITA_forest_function_normalized <- function_normalization(data = GER_ITA_forest_function_raw,
fun_cols = 6:31,
negative = c("soil_cn_ff_10","wue"),
by_group = "country")
GER_ITA_forest_biodiversity <- forest_biodiversity_data[c(49:82,181:229),]
output2 <- MF1_single(func_data = GER_ITA_forest_function_normalized[,6:31], weight = 1,
species_data = GER_ITA_forest_biodiversity)
## BEF relationships: display fitted line of linear mixed model with random slopes and
## random intercepts
output2 <- data.frame(output2, country=rep(GER_ITA_forest_function_normalized$country, each = 6))
MFggplot(output2, model = "LMM.both", by_group="country", caption = "slope")
### Use data from plots in five countries (data in Finland are excluded)
## MF decomposition for all pairs of ecosystems
data("forest_function_data_normalized")
data("forest_biodiversity_data")
forest_function_data_normalized <- filter(forest_function_data_normalized, country != "FIN")
forest_biodiversity_data <- forest_biodiversity_data[-(1:48),]
output3 <- MF2_multiple(func_data = forest_function_data_normalized[,6:32],
species_data = forest_biodiversity_data,
weight = 1,
by_group = "country")
## BEF relationships: display fitted line of linear mixed model with random slopes and
## random intercepts
figure_LMM <- MFggplot(output3, model = "LMM.both", by_group = "country",
caption = "slope")
figure_LMM$corr_uncorrected$ALL
figure_LMM$corr_corrected$ALL
### Use partial data to quickly obtain output
### (Take the first 18 plots in Germany and the last 18 plots in Italy)
## MF decomposition for all pairs of ecosystems for partial data
data("forest_function_data_raw")
data("forest_biodiversity_data")
GER_ITA_forest_function_raw <- filter(forest_function_data_raw,
country=="GER"|country=="ITA")[c(1:18,57:74),]
GER_ITA_forest_function_normalized <- function_normalization(data = GER_ITA_forest_function_raw,
fun_cols = 6:31,
negative = c("soil_cn_ff_10","wue"),
by_group = "country")
GER_ITA_forest_biodiversity <- forest_biodiversity_data[c(49:82,181:229),]
output4 <- MF2_multiple(func_data = GER_ITA_forest_function_normalized[,6:32],
species_data = GER_ITA_forest_biodiversity,
weight = 1,
by_group = "country")
## BEF relationships: display fitted line of linear mixed model with random slopes and
## random intercepts for partial data
figure_LMM_GER_ITA <- MFggplot(output4, model = "LMM.both", by_group = "country",
caption = "slope")
figure_LMM_GER_ITA$corr_uncorrected$ALL
figure_LMM_GER_ITA$corr_corrected$ALL
### Use partial data to calculate multifunctionality based on 3 plots in each country, not by pairs
### (Take the first 3 plots in each country)
data("forest_function_data_raw")
data("forest_biodiversity_data")
forest_function_data_raw_3plots <- forest_function_data_raw[c(1:3,29:31,67:69,103:105,
146:148,174:176),]
forest_function_data_normalized_3plots <-
function_normalization(data = forest_function_data_raw_3plots,
fun_cols = 6:31,
negative = c("soil_cn_ff_10","wue"),
by_group = "country")
forest_biodiversity_data_3plots <-
forest_biodiversity_data[c(1:6,49:52,141:148,230:232,351:355,411:417),]
output5 <- MF2_multiple(func_data = forest_function_data_normalized_3plots[,6:32],
species_data = forest_biodiversity_data_3plots,
weight = 1,
by_group = "country", by_pair = FALSE)
## BEF relationships: display fitted line of linear model
figure_all_plots <- MFggplot(output5, model = "lm", caption = "slope")
figure_all_plots$corr_uncorrected$ALL
figure_all_plots$corr_corrected$ALL
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
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.