diff --git a/Project.toml b/Project.toml
index fe41ff81..ef52eded 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
 name = "FinEtools"
 uuid = "91bb5406-6c9a-523d-811d-0644c4229550"
 authors = ["Petr Krysl <pkrysl@ucsd.edu>"]
-version = "7.0.0"
+version = "7.0.1"
 
 [deps]
 ChunkSplitters = "ae650224-84b6-46f8-82ea-d812ca08434e"
@@ -14,7 +14,6 @@ Printf = "de0858da-6303-5e67-8744-51eddeeeb8d7"
 Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
 SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
 Statistics = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"
-StatsBase = "2913bbd2-ae8a-5f71-8c99-4fb6c76f3a91"
 SymRCM = "286e6d88-80af-4590-acc9-0001b223b9bd"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 Unicode = "4ec0a83e-493e-50e2-b9ac-8f72acf5a8f5"
@@ -26,7 +25,6 @@ DataDrop = "^0.1"
 DelimitedFiles = "^1.7"
 HDF5 = "^0.16"
 LoopVectorization = "^0.12"
-StatsBase = "^0.34"
 SymRCM = "^0.2"
 WriteVTK = "^1"
 julia = "^1.9"
diff --git a/README.md b/README.md
index 7294a67f..4e1f9c74 100644
--- a/README.md
+++ b/README.md
@@ -30,8 +30,7 @@ The package supports application packages, for instance:
 - 12/27/2022: Added T4 and T10 cylinder meshing routines.
 - 11/03/2022: Added extraction of outer surface of solid.
 - 06/29/2022: Reconfigured documentation.
-- 06/15/2022: Added point partitioning.
-- 06/15/2022: Added dual connection matrix implementation.
+- 06/15/2022: Added point partitioning. Added dual connection matrix implementation.
 - 05/03/2022: Updated project to Julia 1.7.2: Julia 1.6 is not supported from version 5.4.0 on.
 - 01/26/2022: Added assembler for sparse diagonal matrices (such as the mass
   matrix for explicit finite element methods).
diff --git a/test/test_miscellaneous.jl b/test/test_miscellaneous.jl
index 6e30f8d3..60651f27 100644
--- a/test/test_miscellaneous.jl
+++ b/test/test_miscellaneous.jl
@@ -1408,8 +1408,7 @@ module mmqtrapm3
 using FinEtools
 using Test
 using FinEtools.AlgoBaseModule: qtrap, qcovariance, qvariance
-using StatsBase
-import Statistics: mean, cov
+using Statistics: mean, cov
 function test()
     ps = vec(
         [0.0201161 0.0567767 0.123692 0.141182 0.196189 0.311646 0.463708 0.798094 0.832338 0.875213 0.880719 0.947033 0.993938 0.99884],