pply(cbind(
lapply(wtr_null$GLM, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
}),
lapply(wtr_null$GOTM, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE) }),
lapply(wtr_null$Simstrat, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)})),
1, function(x) mean(x, na.rm = TRUE))
apply(cbind(
lapply(wtr_null$GLM, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
}),
lapply(wtr_null$GOTM, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE) }),
lapply(wtr_null$Simstrat, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)})),
1, function(x) mean(x, na.rm = TRUE))
lapply(wtr_null$GLM, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})
wtr_null$GLM
lapply(wtr_null$GLM, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})
lapply(wtr_null, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})
str(lapply(wtr_null, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
}))
lapply(wtr_null, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability
apply(cbind(
lapply(wtr_null, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_null, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_null, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE))
df <- data.frame('datetime' = wtr_null$GLM$datetime,
'ssi_null' = apply(cbind(
lapply(wtr_null, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_null, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_null, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_01' = apply(cbind(
lapply(wtr_01, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_01, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_01, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_05' = apply(cbind(
lapply(wtr_05, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_05, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_05, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_1' =apply(cbind(
lapply(wtr_1, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_1, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_1, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_15' = apply(cbind(
lapply(wtr_15, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_15, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_15, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_2' = apply(cbind(
lapply(wtr_2, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_2, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_2, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_25' = apply(cbind(
lapply(wtr_25, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_25, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_25, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_3' = apply(cbind(
lapply(wtr_3, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_3, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_3, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_5' = apply(cbind(
lapply(wtr_5, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_5, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_5, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_10' = apply(cbind(
lapply(wtr_10, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_10, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_10, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE))
)
df_ssi <- data.frame('datetime' = wtr_null$GLM$datetime,
'ssi_null' = apply(cbind(
lapply(wtr_null, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_null, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_null, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_01' = apply(cbind(
lapply(wtr_01, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_01, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_01, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_05' = apply(cbind(
lapply(wtr_05, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_05, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_05, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_1' =apply(cbind(
lapply(wtr_1, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_1, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_1, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_15' = apply(cbind(
lapply(wtr_15, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_15, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_15, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_2' = apply(cbind(
lapply(wtr_2, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_2, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_2, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_25' = apply(cbind(
lapply(wtr_25, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_25, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_25, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_3' = apply(cbind(
lapply(wtr_3, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_3, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_3, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_5' = apply(cbind(
lapply(wtr_5, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_5, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_5, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_10' = apply(cbind(
lapply(wtr_10, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_10, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_10, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE))
)
df <- data.frame('datetime' = wtr_null$GLM$datetime,
'temp_null' = apply(cbind(wtr_null$GLM$wtr_22 - wtr_null$GLM$wtr_2,
wtr_null$GOTM$wtr_22 - wtr_null$GOTM$wtr_2,
wtr_null$Simstrat$wtr_22 - wtr_null$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'temp_01' = apply(cbind(wtr_01$GLM$wtr_22 - wtr_01$GLM$wtr_2,
wtr_01$GOTM$wtr_22 - wtr_01$GOTM$wtr_2,
wtr_01$Simstrat$wtr_22 - wtr_01$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'temp_05' = apply(cbind(wtr_05$GLM$wtr_22 - wtr_05$GLM$wtr_2,
wtr_05$GOTM$wtr_22 - wtr_05$GOTM$wtr_2,
wtr_05$Simstrat$wtr_22 - wtr_05$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'temp_1' = apply(cbind(wtr_1$GLM$wtr_22 - wtr_1$GLM$wtr_2,
wtr_1$GOTM$wtr_22 - wtr_1$GOTM$wtr_2,
wtr_1$Simstrat$wtr_22 - wtr_1$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'temp_15' = apply(cbind(wtr_15$GLM$wtr_22 - wtr_15$GLM$wtr_2,
wtr_15$GOTM$wtr_22 - wtr_15$GOTM$wtr_2,
wtr_15$Simstrat$wtr_22 - wtr_15$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'temp_2' = apply(cbind(wtr_2$GLM$wtr_22 - wtr_2$GLM$wtr_2,
wtr_2$GOTM$wtr_22 - wtr_2$GOTM$wtr_2,
wtr_2$Simstrat$wtr_22 - wtr_2$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'temp_25' = apply(cbind(wtr_25$GLM$wtr_22 - wtr_25$GLM$wtr_2,
wtr_25$GOTM$wtr_22 - wtr_25$GOTM$wtr_2,
wtr_25$Simstrat$wtr_22 - wtr_25$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'temp_3' = apply(cbind(wtr_3$GLM$wtr_22 - wtr_3$GLM$wtr_2,
wtr_3$GOTM$wtr_22 - wtr_3$GOTM$wtr_2,
wtr_3$Simstrat$wtr_22 - wtr_3$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'temp_5' = apply(cbind(wtr_5$GLM$wtr_22 - wtr_5$GLM$wtr_2,
wtr_5$GOTM$wtr_22 - wtr_5$GOTM$wtr_2,
wtr_5$Simstrat$wtr_22 - wtr_5$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'temp_10' = apply(cbind(wtr_10$GLM$wtr_22 - wtr_10$GLM$wtr_2,
wtr_10$GOTM$wtr_22 - wtr_10$GOTM$wtr_2,
wtr_10$Simstrat$wtr_22 - wtr_10$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE))
)
df <- data.frame('datetime' = wtr_null$GLM$datetime,
'dens_null' = apply(cbind(density_null$GLM$wtr_22 - density_null$GLM$wtr_2,
density_null$GOTM$wtr_22 - density_null$GOTM$wtr_2,
density_null$Simstrat$wtr_22 - density_null$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'dens_01' = apply(cbind(density_01$GLM$wtr_22 - density_01$GLM$wtr_2,
density_01$GOTM$wtr_22 - density_01$GOTM$wtr_2,
density_01$Simstrat$wtr_22 - density_01$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'dens_05' = apply(cbind(density_05$GLM$wtr_22 - density_05$GLM$wtr_2,
density_05$GOTM$wtr_22 - density_05$GOTM$wtr_2,
density_05$Simstrat$wtr_22 - density_05$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'dens_1' = apply(cbind(density_1$GLM$wtr_22 - density_1$GLM$wtr_2,
density_1$GOTM$wtr_22 - density_1$GOTM$wtr_2,
density_1$Simstrat$wtr_22 - density_1$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'dens_15' = apply(cbind(density_15$GLM$wtr_22 - density_15$GLM$wtr_2,
density_15$GOTM$wtr_22 - density_15$GOTM$wtr_2,
density_15$Simstrat$wtr_22 - density_15$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'dens_2' = apply(cbind(density_2$GLM$wtr_22 - density_2$GLM$wtr_2,
density_2$GOTM$wtr_22 - density_2$GOTM$wtr_2,
density_2$Simstrat$wtr_22 - density_2$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'dens_25' = apply(cbind(density_25$GLM$wtr_22 - density_25$GLM$wtr_2,
density_25$GOTM$wtr_22 - density_25$GOTM$wtr_2,
density_25$Simstrat$wtr_22 - density_25$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'dens_3' = apply(cbind(density_3$GLM$wtr_22 - density_3$GLM$wtr_2,
density_3$GOTM$wtr_22 - density_3$GOTM$wtr_2,
density_3$Simstrat$wtr_22 - density_3$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'dens_5' = apply(cbind(density_5$GLM$wtr_22 - density_5$GLM$wtr_2,
density_5$GOTM$wtr_22 - density_5$GOTM$wtr_2,
density_5$Simstrat$wtr_22 - density_5$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'dens_10' = apply(cbind(density_10$GLM$wtr_22 - density_10$GLM$wtr_2,
density_10$GOTM$wtr_22 - density_10$GOTM$wtr_2,
density_10$Simstrat$wtr_22 - density_10$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE))
)
df_wtr <- data.frame('datetime' = wtr_null$GLM$datetime,
'temp_null' = apply(cbind(wtr_null$GLM$wtr_22 - wtr_null$GLM$wtr_2,
wtr_null$GOTM$wtr_22 - wtr_null$GOTM$wtr_2,
wtr_null$Simstrat$wtr_22 - wtr_null$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'temp_01' = apply(cbind(wtr_01$GLM$wtr_22 - wtr_01$GLM$wtr_2,
wtr_01$GOTM$wtr_22 - wtr_01$GOTM$wtr_2,
wtr_01$Simstrat$wtr_22 - wtr_01$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'temp_05' = apply(cbind(wtr_05$GLM$wtr_22 - wtr_05$GLM$wtr_2,
wtr_05$GOTM$wtr_22 - wtr_05$GOTM$wtr_2,
wtr_05$Simstrat$wtr_22 - wtr_05$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'temp_1' = apply(cbind(wtr_1$GLM$wtr_22 - wtr_1$GLM$wtr_2,
wtr_1$GOTM$wtr_22 - wtr_1$GOTM$wtr_2,
wtr_1$Simstrat$wtr_22 - wtr_1$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'temp_15' = apply(cbind(wtr_15$GLM$wtr_22 - wtr_15$GLM$wtr_2,
wtr_15$GOTM$wtr_22 - wtr_15$GOTM$wtr_2,
wtr_15$Simstrat$wtr_22 - wtr_15$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'temp_2' = apply(cbind(wtr_2$GLM$wtr_22 - wtr_2$GLM$wtr_2,
wtr_2$GOTM$wtr_22 - wtr_2$GOTM$wtr_2,
wtr_2$Simstrat$wtr_22 - wtr_2$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'temp_25' = apply(cbind(wtr_25$GLM$wtr_22 - wtr_25$GLM$wtr_2,
wtr_25$GOTM$wtr_22 - wtr_25$GOTM$wtr_2,
wtr_25$Simstrat$wtr_22 - wtr_25$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'temp_3' = apply(cbind(wtr_3$GLM$wtr_22 - wtr_3$GLM$wtr_2,
wtr_3$GOTM$wtr_22 - wtr_3$GOTM$wtr_2,
wtr_3$Simstrat$wtr_22 - wtr_3$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'temp_5' = apply(cbind(wtr_5$GLM$wtr_22 - wtr_5$GLM$wtr_2,
wtr_5$GOTM$wtr_22 - wtr_5$GOTM$wtr_2,
wtr_5$Simstrat$wtr_22 - wtr_5$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE)),
'temp_10' = apply(cbind(wtr_10$GLM$wtr_22 - wtr_10$GLM$wtr_2,
wtr_10$GOTM$wtr_22 - wtr_10$GOTM$wtr_2,
wtr_10$Simstrat$wtr_22 - wtr_10$Simstrat$wtr_2), 1, function(x) mean(x, na.rm = TRUE))
)
df_ssi <- data.frame('datetime' = wtr_null$GLM$datetime,
'ssi_null' = apply(cbind(
lapply(wtr_null, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_null, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_null, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_01' = apply(cbind(
lapply(wtr_01, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_01, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_01, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_05' = apply(cbind(
lapply(wtr_05, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_05, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_05, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_1' =apply(cbind(
lapply(wtr_1, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_1, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_1, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_15' = apply(cbind(
lapply(wtr_15, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_15, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_15, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_2' = apply(cbind(
lapply(wtr_2, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_2, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_2, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_25' = apply(cbind(
lapply(wtr_25, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_25, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_25, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_3' = apply(cbind(
lapply(wtr_3, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_3, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_3, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_5' = apply(cbind(
lapply(wtr_5, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_5, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_5, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE)),
'ssi_10' = apply(cbind(
lapply(wtr_10, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GLM$schmidt.stability,
lapply(wtr_10, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$GOTM$schmidt.stability,
lapply(wtr_10, function(x) {
ts.schmidt.stability(x, bathy = bathy, na.rm = TRUE)
})$Simstrat$schmidt.stability),
1, function(x) mean(x, na.rm = TRUE))
)
m.df <- reshape2::melt(df, id = 'datetime')
g <- ggplot(m.df) +
geom_line(aes(datetime, (value), col = variable)) +
ylab('Density differences [g/kg]') +
xlab('') +
geom_hline(yintercept = 1.1, linetype = 'dashed') +
theme_minimal(); g
m.df_ssi <- reshape2::melt(df_ssi, id = 'datetime')
g <- ggplot(m.df_ssi) +
geom_line(aes(datetime, (value), col = variable)) +
ylab('SSI [J/m2]') +
xlab('') +
theme_minimal(); g
df_ssi
df_ssi[,2:ncol(df_ssi)] <- apply(df_ssi, 2, function(x) x-df_ssi[,2])
df_ssi[,2]
apply(df_ssi, 2, function(x) x-df_ssi[,2])
apply(df_ssi, 2, function(x) x)
df_ssi[,2:ncol(df_ssi)]
str(df_ssi[,2:ncol(df_ssi)] )
as.matrix(df_ssi[,2:ncol(df_ssi)] )
df_ssi[,2:ncol(df_ssi)] <- apply(as.matrix(df_ssi[,2:ncol(df_ssi)] ), 2, function(x) x-as.matrix(df_ssi[,2:ncol(df_ssi)] )[,1])
df_ssi
m.df_ssi <- reshape2::melt(df_ssi, id = 'datetime')
g <- ggplot(m.df_ssi) +
geom_line(aes(datetime, (value), col = variable)) +
ylab('SSI [J/m2]') +
xlab('') +
theme_minimal(); g
g <- ggplot(m.df_ssi) +
geom_line(aes(datetime, (value), col = variable)) +
ylab('SSI normalized on null [J/m2]') +
xlab('') +
theme_minimal(); g
ggsave('output/schmidtensemble.png', g,  dpi = 300,width = 250,height = 100, units = 'mm')
