Avoid using deprecated indexing

Author: BenChung

Project.toml

@@ -1,25 +1,28 @@
 name = "ModelTesting"
 uuid = "726cd793-0c59-4a7a-97d2-1baad301d093"
 authors = ["Ben Chung <benjamin.chung@juliahub.com> and contributors"]
-version = "1.0.6-DEV"
+version = "1.0.7-DEV"
 
 [deps]
 CSV = "336ed68f-0bac-5ca0-87d4-7b16caf5d00b"
 DataFrames = "a93c6f00-e57d-5684-b7b6-d8193f3e46c0"
+Distributions = "31c24e10-a181-5473-b8eb-7969acd0382f"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
+LowLevelParticleFilters = "d9d29d28-c116-5dba-9239-57a5fe23875b"
 ModelingToolkit = "961ee093-0014-501f-94e3-6117800e7a78"
 SciMLBase = "0bca4576-84f4-4d90-8ffe-ffa030f20462"
+SeeToDee = "1c904df7-48cd-41e7-921b-d889ed9a470c"
 StatsBase = "2913bbd2-ae8a-5f71-8c99-4fb6c76f3a91"
 SymbolicIndexingInterface = "2efcf032-c050-4f8e-a9bb-153293bab1f5"
 
 [compat]
-julia = "1.10"
 SymbolicIndexingInterface = "^0.3.15"
+julia = "1.10"
 
 [extras]
-Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
-ModelingToolkitStandardLibrary = "16a59e39-deab-5bd0-87e4-056b12336739"
 DifferentialEquations = "0c46a032-eb83-5123-abaf-570d42b7fbaa"
+ModelingToolkitStandardLibrary = "16a59e39-deab-5bd0-87e4-056b12336739"
+Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 
 [targets]
 test = ["Test", "ModelingToolkitStandardLibrary", "DifferentialEquations"]

src/test/compare.jl

@@ -35,7 +35,7 @@ function compare(
     @assert all(b ∈ names(reference) for b in reference_basis) "The reference basis must be a subset of the columns in the reference data"
     if new_sol.dense
         dat = new_sol(reference[:, :timestamp], idxs=basis)
-        obs = [[dat[j][i] for j=1:nrow(reference)] for i in eachindex(basis)]
+        obs = [[dat[i, j] for j=1:nrow(reference)] for i in eachindex(basis)]
         ref = collect.(eachcol(reference[:, Not(:timestamp)]))
         cmp(new_container, output_names, basis, reference[:, :timestamp], obs, ref)
     else