Current preprocessing_model nests systems [j] within stories [i].
As a result, when running the analysis for multiple systems:
out <- frbca(eal_multiple_systems, cost_multiple_systems, params)
the output out is nested by stories rather than systems. In other words, out is a list with the length length(out) = number of unique story heights represented in the data.
It should be the other way around: nest stories within systems, so that length(out) = number of unique structural systems.
Current
preprocessing_modelnests systems[j]within stories[i].As a result, when running the analysis for multiple systems:
the output
outis nested by stories rather than systems. In other words,outis a list with the lengthlength(out) =number of unique story heights represented in the data.It should be the other way around: nest stories within systems, so that
length(out) =number of unique structural systems.