theme_minimal() +
geom_point(data=df_FEA_smilo, aes(shape = Category,
size = 3,
color = color_common[1]),
alpha=0.8) +
geom_line(data=df_FEA_smilo, size = 1.2, alpha = 0.8, color = color_common[1]) +
labs( x = "Mandible length",
y = "Mechanical efficiency") +
scale_x_continuous(limits = limit_mandible) +
geom_ribbon(data = df_FEA_smilo, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = color_common[1]), alpha=0.1)  +
theme(axis.line = element_line(color = "darkgrey", linewidth = 0.5),legend.position = "none")
# Add lion
plot_ME <- plot_ME  +
geom_point(data = df_FEA_lion, aes(x=MandibleL, y=ME,
shape = Category,
size = 3,
color = color_common[3]),
alpha=0.8) +
scale_shape_manual(values=shape_to_plot) +
geom_ribbon(data = df_FEA_lion, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = color_common[3]), alpha=0.1)  +
geom_line(data = df_FEA_lion,aes(x=MandibleL, y=ME),  size = 1.2, alpha = 0.8, color = color_common[3], linetype = "dotted")
plot_ME
plot_ME <- ggplot(data=df_FEA_smilo, aes(x=MandibleL, y=ME)) +
theme_minimal() +
geom_point(data=df_FEA_smilo, aes(shape = Category,
size = 3,
color = color_common[1]),
alpha=0.8) +
geom_line(data=df_FEA_smilo, size = 1.2, alpha = 0.8, color = color_common[1]) +
labs( x = "Mandible length",
y = "Mechanical efficiency") +
scale_x_continuous(limits = limit_mandible) +
geom_ribbon(data = df_FEA_smilo, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = color_common[1]), alpha=0.1)  +
theme(axis.line = element_line(color = "darkgrey", linewidth = 0.5),legend.position = "none")
# Add lion
plot_ME <- plot_ME  +
geom_point(data = df_FEA_lion, aes(x=MandibleL, y=ME,
shape = Category,
size = 3,
color = color_common[3]),
alpha=0.8)
plot_ME
plot_ME <- plot_ME  +
geom_point(data = df_FEA_lion, aes(x=MandibleL, y=ME,
shape = Category,
size = 3,
color = color_common[3]),
alpha=0.8) +
scale_shape_manual(values=shape_to_plot) +
geom_ribbon(data = df_FEA_lion, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = color_common[3]), alpha=0.1)  +
geom_line(data = df_FEA_lion,aes(x=MandibleL, y=ME),  size = 1.2, alpha = 0.8, color = color_common[3], linetype = "dotted")
# Adjusted internal energy
# Smilodon
plot_adjIE <- ggplot(data=df_FEA_smilo, aes(x=MandibleL, y=AdjustedIE)) +
theme_minimal() +
geom_point(aes(shape = Category, size = 3, color = color_common[1]), alpha=0.8) +
geom_line(size = 1.2, alpha = 0.8, color = color_common[1]) +
labs( x = "Mandible length",
y = "Adjusted internal energy") +
scale_x_continuous(limits = limit_mandible) +
geom_ribbon(data = df_FEA_smilo, aes(ymin=IEneg95_CI, ymax=IEpos_95_CI, fill = color_common[1]), alpha=0.1)  +
theme(axis.line = element_line(color = "darkgrey", linewidth = 0.5),legend.position = "none")
plot_ME <- plot_ME  +
geom_point(data = df_FEA_lion, aes(x=MandibleL, y=ME,
shape = Category,
size = 3,
color = color_common[3]),
alpha=0.8) +
scale_shape_manual(values=shape_to_plot) +
geom_ribbon(data = df_FEA_lion, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = color_common[3]), alpha=0.1)  +
geom_line(data = df_FEA_lion,aes(x=MandibleL, y=ME),  size = 1.2, alpha = 0.8, color = color_common[3], linetype = "dotted") +
theme(axis.line = element_line(color = "darkgrey", linewidth = 0.5),legend.position = "none")
# Adjusted internal energy
# Smilodon
plot_adjIE <- ggplot(data=df_FEA_smilo, aes(x=MandibleL, y=AdjustedIE)) +
theme_minimal() +
geom_point(aes(shape = Category, size = 3, color = color_common[1]), alpha=0.8) +
geom_line(size = 1.2, alpha = 0.8, color = color_common[1]) +
labs( x = "Mandible length",
y = "Adjusted internal energy") +
scale_x_continuous(limits = limit_mandible) +
geom_ribbon(data = df_FEA_smilo, aes(ymin=IEneg95_CI, ymax=IEpos_95_CI, fill = color_common[1]), alpha=0.1)
plot_ME <- ggplot(data=df_FEA_smilo, aes(x=MandibleL, y=ME)) +
theme_minimal() +
geom_point(data=df_FEA_smilo, aes(shape = Category,
size = 3,
color = color_common[1]),
alpha=0.8) +
geom_line(data=df_FEA_smilo, size = 1.2, alpha = 0.8, color = color_common[1]) +
labs( x = "Mandible length",
y = "Mechanical efficiency") +
scale_x_continuous(limits = limit_mandible) +
geom_ribbon(data = df_FEA_smilo, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = color_common[1]), alpha=0.1)  +
theme(axis.line = element_line(color = "darkgrey", linewidth = 0.5),legend.position = "none")
# Add lion
plot_ME <- plot_ME  +
geom_point(data = df_FEA_lion, aes(x=MandibleL, y=ME,
shape = Category,
size = 3,
color = color_common[3]),
alpha=0.8) +
scale_shape_manual(values=shape_to_plot) +
geom_ribbon(data = df_FEA_lion, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = color_common[3]), alpha=0.1)  +
geom_line(data = df_FEA_lion,aes(x=MandibleL, y=ME),  size = 1.2, alpha = 0.8, color = color_common[3], linetype = "dotted") +
plot_ME
plot_ME
plot_ME <- ggplot(data=df_FEA_smilo, aes(x=MandibleL, y=ME)) +
theme_minimal() +
geom_point(data=df_FEA_smilo, aes(shape = Category,
size = 3,
color = color_common[1]),
alpha=0.8) +
geom_line(data=df_FEA_smilo, size = 1.2, alpha = 0.8, color = color_common[1]) +
labs( x = "Mandible length",
y = "Mechanical efficiency") +
scale_x_continuous(limits = limit_mandible) +
geom_ribbon(data = df_FEA_smilo, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = color_common[1]), alpha=0.1)  +
theme(axis.line = element_line(color = "darkgrey", linewidth = 0.5),legend.position = "none")
# Add lion
plot_ME <- plot_ME  +
geom_point(data = df_FEA_lion, aes(x=MandibleL, y=ME,
shape = Category,
size = 3,
color = color_common[3]),
alpha=0.8) +
scale_shape_manual(values=shape_to_plot) +
geom_ribbon(data = df_FEA_lion, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = color_common[3]), alpha=0.1)  +
geom_line(data = df_FEA_lion,aes(x=MandibleL, y=ME),  size = 1.2, alpha = 0.8, color = color_common[3], linetype = "dotted")
plot_ME
df_FEA_lion <- cbind(FEA_lion[,1:5],
'MEpos_95_CI' = (FEA_lion$ME + FEA_lion$ME_95CI),
'MEneg95_CI' = (FEA_lion$ME - FEA_lion$ME_95CI),
'IEpos_95_CI' = (FEA_lion$AdjustedIE + FEA_lion$AdjustedIE_95CI),
'IEneg95_CI' = (FEA_lion$AdjustedIE - FEA_lion$AdjustedIE_95CI),
'species' = rep('Lion', legnth(FEA_lion$Specimen)))
df_FEA_lion <- cbind(FEA_lion[,1:5],
'MEpos_95_CI' = (FEA_lion$ME + FEA_lion$ME_95CI),
'MEneg95_CI' = (FEA_lion$ME - FEA_lion$ME_95CI),
'IEpos_95_CI' = (FEA_lion$AdjustedIE + FEA_lion$AdjustedIE_95CI),
'IEneg95_CI' = (FEA_lion$AdjustedIE - FEA_lion$AdjustedIE_95CI),
'species' = rep('Lion', length(FEA_lion$Specimen)))
df_FEA_lion
df_FEA_smilo <- cbind(FEA_smilo[,1:5],
'MEpos_95_CI' = (FEA_smilo$ME + FEA_smilo$ME_95CI),
'MEneg95_CI' = (FEA_smilo$ME - FEA_smilo$ME_95CI),
'IEpos_95_CI' = (FEA_smilo$AdjustedIE + FEA_smilo$AdjustedIE_95CI),
'IEneg95_CI' = (FEA_smilo$AdjustedIE - FEA_smilo$AdjustedIE_95CI),
'species' = rep('Smilo', length(FEA_lion$Specimen)))
df_FEA_smilo <- cbind(FEA_smilo[,1:5],
'MEpos_95_CI' = (FEA_smilo$ME + FEA_smilo$ME_95CI),
'MEneg95_CI' = (FEA_smilo$ME - FEA_smilo$ME_95CI),
'IEpos_95_CI' = (FEA_smilo$AdjustedIE + FEA_smilo$AdjustedIE_95CI),
'IEneg95_CI' = (FEA_smilo$AdjustedIE - FEA_smilo$AdjustedIE_95CI),
'species' = rep('Smilo', length(FEA_smilo$Specimen)))
df_FEA_lion <- cbind(FEA_lion[,1:5],
'MEpos_95_CI' = (FEA_lion$ME + FEA_lion$ME_95CI),
'MEneg95_CI' = (FEA_lion$ME - FEA_lion$ME_95CI),
'IEpos_95_CI' = (FEA_lion$AdjustedIE + FEA_lion$AdjustedIE_95CI),
'IEneg95_CI' = (FEA_lion$AdjustedIE - FEA_lion$AdjustedIE_95CI),
'species' = rep('Lion', length(FEA_lion$Specimen)))
color_common2 <- c("Adult_Smilo" = inferno(3)[2],
"adult_lion" = viridis(3)[2])
# Mechanial efficiency
# Smilodon
plot_ME <- ggplot(data=df_FEA_smilo, aes(x=MandibleL, y=ME)) +
theme_minimal() +
geom_point(data=df_FEA_smilo, aes(shape = Category,
size = 3,
color = species),
alpha=0.8) +
geom_line(data=df_FEA_smilo, size = 1.2, alpha = 0.8, color = species) +
labs( x = "Mandible length",
y = "Mechanical efficiency") +
scale_x_continuous(limits = limit_mandible) +
geom_ribbon(data = df_FEA_smilo, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = species), alpha=0.1)  +
theme(axis.line = element_line(color = "darkgrey", linewidth = 0.5),legend.position = "none")
# Add lion
plot_ME <- plot_ME  +
geom_point(data = df_FEA_lion, aes(x=MandibleL, y=ME,
shape = Category,
size = 3,
color = species,
alpha=0.8)) +
scale_shape_manual(values=shape_to_plot) +
geom_ribbon(data = df_FEA_lion, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = species), alpha=0.1)  +
geom_line(data = df_FEA_lion,aes(x=MandibleL, y=ME),  size = 1.2, alpha = 0.8, color = species, linetype = "dotted") +
scale_color_manual(values=color_common2) +
scale_fill_manual(values=color_common2)
df_FEA_smilo <- cbind(FEA_smilo[,1:5],
'MEpos_95_CI' = (FEA_smilo$ME + FEA_smilo$ME_95CI),
'MEneg95_CI' = (FEA_smilo$ME - FEA_smilo$ME_95CI),
'IEpos_95_CI' = (FEA_smilo$AdjustedIE + FEA_smilo$AdjustedIE_95CI),
'IEneg95_CI' = (FEA_smilo$AdjustedIE - FEA_smilo$AdjustedIE_95CI),
'species' = rep('Smilo', length(FEA_smilo$Specimen)))
df_FEA_lion <- cbind(FEA_lion[,1:5],
'MEpos_95_CI' = (FEA_lion$ME + FEA_lion$ME_95CI),
'MEneg95_CI' = (FEA_lion$ME - FEA_lion$ME_95CI),
'IEpos_95_CI' = (FEA_lion$AdjustedIE + FEA_lion$AdjustedIE_95CI),
'IEneg95_CI' = (FEA_lion$AdjustedIE - FEA_lion$AdjustedIE_95CI),
'species' = rep('Lion', length(FEA_lion$Specimen)))
# Mechanial efficiency
# Smilodon
plot_ME <- ggplot(data=df_FEA_smilo, aes(x=MandibleL, y=ME)) +
theme_minimal() +
geom_point(data=df_FEA_smilo, aes(shape = Category,
size = 3,
color = species),
alpha=0.8) +
geom_line(data=df_FEA_smilo, size = 1.2, alpha = 0.8, color = species) +
labs( x = "Mandible length",
y = "Mechanical efficiency") +
scale_x_continuous(limits = limit_mandible) +
geom_ribbon(data = df_FEA_smilo, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = species), alpha=0.1)  +
theme(axis.line = element_line(color = "darkgrey", linewidth = 0.5),legend.position = "none")
# Add lion
plot_ME <- plot_ME  +
geom_point(data = df_FEA_lion, aes(x=MandibleL, y=ME,
shape = Category,
size = 3,
color = species,
alpha=0.8)) +
scale_shape_manual(values=shape_to_plot) +
geom_ribbon(data = df_FEA_lion, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = species), alpha=0.1)  +
geom_line(data = df_FEA_lion,aes(x=MandibleL, y=ME),  size = 1.2, alpha = 0.8, color = species, linetype = "dotted") +
scale_color_manual(values=color_common2) +
scale_fill_manual(values=color_common2)
df_FEA_smilo
df_FEA_lion
plot_ME <- ggplot(data=df_FEA_smilo, aes(x=MandibleL, y=ME)) +
theme_minimal() +
geom_point(data=df_FEA_smilo, aes(shape = Category,
size = 3,
color = species),
alpha=0.8) +
geom_line(data=df_FEA_smilo, size = 1.2, alpha = 0.8, color = species) +
labs( x = "Mandible length",
y = "Mechanical efficiency") +
scale_x_continuous(limits = limit_mandible) +
geom_ribbon(data = df_FEA_smilo, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = species), alpha=0.1)  +
theme(axis.line = element_line(color = "darkgrey", linewidth = 0.5),legend.position = "none")
# Mechanial efficiency
# Smilodon
plot_ME <- ggplot(data=df_FEA_smilo, aes(x=MandibleL, y=ME)) +
theme_minimal() +
geom_point(data=df_FEA_smilo, aes(shape = Category,
size = 3,
color = species),
alpha=0.8) +
geom_line(data=df_FEA_smilo, size = 1.2, alpha = 0.8, color = species) +
labs( x = "Mandible length",
y = "Mechanical efficiency") +
scale_x_continuous(limits = limit_mandible) +
geom_ribbon(data = df_FEA_smilo, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = species), alpha=0.1)  +
theme(axis.line = element_line(color = "darkgrey", linewidth = 0.5),legend.position = "none")
# Mechanial efficiency
# Smilodon
plot_ME <- ggplot(data=df_FEA_smilo, aes(x=MandibleL, y=ME)) +
theme_minimal() +
geom_point(data=df_FEA_smilo, aes(shape = Category,
size = 3,
color = df_FEA_smilo$species),
alpha=0.8) +
geom_line(data=df_FEA_smilo, size = 1.2, alpha = 0.8, color = df_FEA_smilo$species) +
labs( x = "Mandible length",
y = "Mechanical efficiency") +
scale_x_continuous(limits = limit_mandible) +
geom_ribbon(data = df_FEA_smilo, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = df_FEA_smilo$species), alpha=0.1)  +
theme(axis.line = element_line(color = "darkgrey", linewidth = 0.5),legend.position = "none")
plot_ME
color_common2 <- c("Smilo" = inferno(3)[2],
"Lion" = viridis(3)[2])
# Mechanial efficiency
# Smilodon
plot_ME <- ggplot(data=df_FEA_smilo, aes(x=MandibleL, y=ME)) +
theme_minimal() +
geom_point(data=df_FEA_smilo, aes(shape = Category,
size = 3,
color = species),
alpha=0.8) +
geom_line(data=df_FEA_smilo, size = 1.2, alpha = 0.8, color = species) +
labs( x = "Mandible length",
y = "Mechanical efficiency") +
scale_x_continuous(limits = limit_mandible) +
geom_ribbon(data = df_FEA_smilo, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = species), alpha=0.1)  +
theme(axis.line = element_line(color = "darkgrey", linewidth = 0.5),legend.position = "none")
df_FEA_smilo <- cbind(FEA_smilo[,1:5],
'MEpos_95_CI' = (FEA_smilo$ME + FEA_smilo$ME_95CI),
'MEneg95_CI' = (FEA_smilo$ME - FEA_smilo$ME_95CI),
'IEpos_95_CI' = (FEA_smilo$AdjustedIE + FEA_smilo$AdjustedIE_95CI),
'IEneg95_CI' = (FEA_smilo$AdjustedIE - FEA_smilo$AdjustedIE_95CI),
'Species' = rep('Smilo', length(FEA_smilo$Specimen)))
df_FEA_lion <- cbind(FEA_lion[,1:5],
'MEpos_95_CI' = (FEA_lion$ME + FEA_lion$ME_95CI),
'MEneg95_CI' = (FEA_lion$ME - FEA_lion$ME_95CI),
'IEpos_95_CI' = (FEA_lion$AdjustedIE + FEA_lion$AdjustedIE_95CI),
'IEneg95_CI' = (FEA_lion$AdjustedIE - FEA_lion$AdjustedIE_95CI),
'Species' = rep('Lion', length(FEA_lion$Specimen)))
df_FEA_lion
View(df_FEA_lion)
# Mechanial efficiency
# Smilodon
plot_ME <- ggplot(data=df_FEA_smilo, aes(x=MandibleL, y=ME)) +
theme_minimal() +
geom_point(data=df_FEA_smilo, aes(shape = Category,
size = 3,
color = Species),
alpha=0.8) +
geom_line(data=df_FEA_smilo, size = 1.2, alpha = 0.8, color = Species) +
labs( x = "Mandible length",
y = "Mechanical efficiency") +
scale_x_continuous(limits = limit_mandible) +
geom_ribbon(data = df_FEA_smilo, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = Species), alpha=0.1)  +
theme(axis.line = element_line(color = "darkgrey", linewidth = 0.5),legend.position = "none")
df_FEA_smilo
View(df_FEA_smilo)
# Mechanial efficiency
# Smilodon
plot_ME <- ggplot(data=df_FEA_smilo, aes(x=MandibleL, y=ME)) +
theme_minimal() +
geom_point(data=df_FEA_smilo, aes(shape = Category,
size = 3,
color = Species),
alpha=0.8) +
geom_line(data=df_FEA_smilo, size = 1.2, alpha = 0.8, color = Species) +
labs( x = "Mandible length",
y = "Mechanical efficiency") +
scale_x_continuous(limits = limit_mandible) +
geom_ribbon(data = df_FEA_smilo, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = Species), alpha=0.1)  +
theme(axis.line = element_line(color = "darkgrey", linewidth = 0.5),legend.position = "none")
plot_ME <- ggplot(data=df_FEA_smilo, aes(x=MandibleL, y=ME)) +
theme_minimal() +
geom_point(data=df_FEA_smilo, aes(shape = Category,
size = 3,
color = color_common[3]),
alpha=0.8) +
geom_line(data=df_FEA_smilo, size = 1.2, alpha = 0.8, color = color_common[3]) +
labs( x = "Mandible length",
y = "Mechanical efficiency") +
scale_x_continuous(limits = limit_mandible) +
geom_ribbon(data = df_FEA_smilo, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = color_common[3]), alpha=0.1)  +
theme(axis.line = element_line(color = "darkgrey", linewidth = 0.5),legend.position = "none")
plot_ME
# Mechanial efficiency
# Smilodon
plot_ME <- ggplot(data=df_FEA_smilo, aes(x=MandibleL, y=ME)) +
theme_minimal() +
geom_point(data=df_FEA_smilo, aes(shape = Category,
size = 3,
color = color_common[3]),
alpha=0.8) +
geom_line(data=df_FEA_smilo, size = 1.2, alpha = 0.8, color = color_common[3]) +
labs( x = "Mandible length",
y = "Mechanical efficiency") +
scale_x_continuous(limits = limit_mandible) +
geom_ribbon(data = df_FEA_smilo, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, color = color_common[3]), alpha=0.1)  +
theme(axis.line = element_line(color = "darkgrey", linewidth = 0.5),legend.position = "none")
# Add lion
plot_ME <- plot_ME  +
geom_point(data = df_FEA_lion, aes(x=MandibleL, y=ME,
shape = Category,
size = 3,
color = species,
alpha=0.8)) +
scale_shape_manual(values=shape_to_plot) +
geom_ribbon(data = df_FEA_lion, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = species), alpha=0.1)  +
geom_line(data = df_FEA_lion,aes(x=MandibleL, y=ME),  size = 1.2, alpha = 0.8, color = species, linetype = "dotted") +
scale_color_manual(values=color_common2) +
scale_fill_manual(values=color_common2)
plot_ME
# Mechanial efficiency
# Smilodon
plot_ME <- ggplot(data=df_FEA_smilo, aes(x=MandibleL, y=ME)) +
theme_minimal() +
geom_point(aes(shape = Category,
size = 3,
color = color_common[3]),
alpha=0.8) +
geom_line(size = 1.2, alpha = 0.8, color = color_common[3]) +
labs( x = "Mandible length",
y = "Mechanical efficiency") +
scale_x_continuous(limits = limit_mandible) +
geom_ribbon(aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, color = color_common[3]), alpha=0.1)  +
theme(axis.line = element_line(color = "darkgrey", linewidth = 0.5),legend.position = "none")
plot_ME
# Mechanial efficiency
# Smilodon
plot_ME <- ggplot(data=df_FEA_smilo, aes(x=MandibleL, y=ME)) +
theme_minimal() +
geom_point(aes(shape = Category,
size = 3,
color = color_common[3]),
alpha=0.8) +
geom_line(size = 1.2, alpha = 0.8, color = color_common[3]) +
labs( x = "Mandible length",
y = "Mechanical efficiency") +
scale_x_continuous(limits = limit_mandible) +
geom_ribbon(aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = color_common[3]), alpha=0.1)  +
theme(axis.line = element_line(color = "darkgrey", linewidth = 0.5),legend.position = "none")
plot_ME
plot_ME <- ggplot(data=df_FEA_smilo, aes(x=MandibleL, y=ME)) +
theme_minimal() +
geom_point(aes(shape = Category,
size = 3,
color = color_common[3]),
alpha=0.8) +
geom_line(size = 1.2, alpha = 0.8, color = color_common[3]) +
labs( x = "Mandible length",
y = "Mechanical efficiency") +
scale_x_continuous(limits = limit_mandible) +
geom_ribbon(aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = color_common[3]), alpha=0.1)  +
theme(axis.line = element_line(color = "darkgrey", linewidth = 0.5),legend.position = "none")
# Add lion
plot_ME <- plot_ME  +
geom_point(data = df_FEA_lion, aes(x=MandibleL, y=ME,
shape = Category,
size = 3,
color = color_common[1],
alpha=0.8)) +
scale_shape_manual(values=shape_to_plot) +
geom_ribbon(data = df_FEA_lion, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = color_common[1]), alpha=0.1)  +
geom_line(data = df_FEA_lion,aes(x=MandibleL, y=ME),  size = 1.2, alpha = 0.8, color = color_common[1], linetype = "dotted")
plot_ME
plot_ME <- plot_ME  +
geom_point(data = df_FEA_lion, aes(x=MandibleL, y=ME,
shape = Category,
size = 3,
color = color_common[1],
alpha=0.8)) +
scale_shape_manual(values=shape_to_plot) +
geom_ribbon(data = df_FEA_lion, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = color_common[1]), alpha=0.1)  +
geom_line(data = df_FEA_lion,aes(x=MandibleL, y=ME),  size = 1.2, alpha = 0.8, color = color_common[1], linetype = "dotted")
# Adjusted internal energy
# Smilodon
plot_adjIE <- ggplot(data=df_FEA_smilo, aes(x=MandibleL, y=AdjustedIE)) +
theme_minimal() +
geom_point(aes(shape = Category, size = 3, color = color_common[1]), alpha=0.8) +
geom_line(size = 1.2, alpha = 0.8, color = color_common[1]) +
labs( x = "Mandible length",
y = "Adjusted internal energy") +
scale_x_continuous(limits = limit_mandible) +
geom_ribbon(data = df_FEA_smilo, aes(ymin=IEneg95_CI, ymax=IEpos_95_CI, fill = color_common[1]), alpha=0.1)
# Add lion
plot_adjIE <- plot_adjIE +
geom_point(data = df_FEA_lion, aes(x=MandibleL, y=AdjustedIE, shape = Category, size = 3, color = color_common[3]), alpha=0.8) +
scale_shape_manual(values=shape_to_plot) +
scal
plot_ME <- plot_ME  +
geom_point(data = df_FEA_lion, aes(x=MandibleL, y=ME,
shape = Category,
size = 3,
color = color_common[1],
alpha=0.8)) +
scale_shape_manual(values=shape_to_plot) +
geom_ribbon(data = df_FEA_lion, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = color_common[1]), alpha=0.1)  +
geom_line(data = df_FEA_lion,aes(x=MandibleL, y=ME),  size = 1.2, alpha = 0.8, color = color_common[1], linetype = "dotted")
# Adjusted internal energy
# Smilodon
plot_adjIE <- ggplot(data=df_FEA_smilo, aes(x=MandibleL, y=AdjustedIE)) +
theme_minimal() +
geom_point(aes(shape = Category, size = 3, color = color_common[1]), alpha=0.8) +
geom_line(size = 1.2, alpha = 0.8, color = color_common[1]) +
labs( x = "Mandible length",
y = "Adjusted internal energy") +
scale_x_continuous(limits = limit_mandible) +
geom_ribbon(data = df_FEA_smilo, aes(ymin=IEneg95_CI, ymax=IEpos_95_CI, fill = color_common[1]), alpha=0.1)
# Add lion
plot_adjIE <- plot_adjIE +
geom_point(data = df_FEA_lion, aes(x=MandibleL, y=AdjustedIE, shape = Category, size = 3, color = color_common[3]), alpha=0.8) +
scale_shape_manual(values=shape_to_plot) +
geom_ribbon(data = df_FEA_lion, aes(ymin=IEneg95_CI, ymax=IEpos_95_CI, fill = color_common[3]), alpha=0.1)  +
geom_line(data = df_FEA_lion,aes(x=MandibleL, y=AdjustedIE),  size = 1.2, alpha = 0.8, color = color_common[3], linetype = "dotted")
ggarrange(plot_ME, plot_adjIE, nrow = 2)
plot_ME <- ggplot(data=df_FEA_smilo, aes(x=MandibleL, y=ME)) +
theme_minimal() +
geom_point(aes(shape = Category,
size = 3,
color = color_common[3]),
alpha=0.8) +
geom_line(size = 1.2, alpha = 0.8, color = color_common[3]) +
labs( x = "Mandible length",
y = "Mechanical efficiency") +
scale_x_continuous(limits = limit_mandible) +
geom_ribbon(aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = color_common[3]), alpha=0.1)  +
theme(axis.line = element_line(color = "darkgrey", linewidth = 0.5),legend.position = "none")
# Add lion
plot_ME <- plot_ME  +
geom_point(data = df_FEA_lion, aes(x=MandibleL, y=ME,
shape = Category,
size = 3,
color = color_common[1],
alpha=0.8)) +
scale_shape_manual(values=shape_to_plot) +
geom_ribbon(data = df_FEA_lion, aes(ymin=MEneg95_CI, ymax=MEpos_95_CI, fill = color_common[1]), alpha=0.1)  +
geom_line(data = df_FEA_lion,aes(x=MandibleL, y=ME),  size = 1.2, alpha = 0.8, color = color_common[1], linetype = "dotted")
# Adjusted internal energy
# Smilodon
plot_adjIE <- ggplot(data=df_FEA_smilo, aes(x=MandibleL, y=AdjustedIE)) +
theme_minimal() +
geom_point(aes(shape = Category, size = 3, color = color_common[1]), alpha=0.8) +
geom_line(size = 1.2, alpha = 0.8, color = color_common[1]) +
labs( x = "Mandible length",
y = "Adjusted internal energy") +
scale_x_continuous(limits = limit_mandible) +
geom_ribbon(data = df_FEA_smilo, aes(ymin=IEneg95_CI, ymax=IEpos_95_CI, fill = color_common[1]), alpha=0.1)
# Add lion
plot_adjIE <- plot_adjIE +
geom_point(data = df_FEA_lion, aes(x=MandibleL, y=AdjustedIE, shape = Category, size = 3, color = color_common[3]), alpha=0.8) +
scale_shape_manual(values=shape_to_plot) +
geom_ribbon(data = df_FEA_lion, aes(ymin=IEneg95_CI, ymax=IEpos_95_CI, fill = color_common[3]), alpha=0.1)  +
geom_line(data = df_FEA_lion,aes(x=MandibleL, y=AdjustedIE),  size = 1.2, alpha = 0.8, color = color_common[3], linetype = "dotted") +
theme(axis.line = element_line(color = "darkgrey", linewidth = 0.5),legend.position = "none")
ggarrange(plot_ME, plot_adjIE, nrow = 2)
