Add more tests

javierbarbero · Jul 18, 2020 · 478f17a · 478f17a
1 parent 607a458
commit 478f17a
Show file tree

Hide file tree

Showing 12 changed files with 250 additions and 16 deletions.
diff --git a/src/deaadd.jl b/src/deaadd.jl
@@ -242,7 +242,7 @@ end
 
 function deaadd(X::Vector, Y::Matrix, model::Symbol = :Default; orient::Symbol = :Graph,
     rts::Symbol = :VRS,
-    rho::Vector = Array{Float64}(undef, 0), rhoY::Matrix = Array{Float64}(undef, 0, 0),
+    rhoX::Vector = Array{Float64}(undef, 0), rhoY::Matrix = Array{Float64}(undef, 0, 0),
     Xref::Vector = X, Yref::Matrix = Y,
     disposX::Symbol = :Strong, disposY::Symbol = :Strong,
     names::Vector{String} = Array{String}(undef, 0))::AdditiveDEAModel
@@ -291,7 +291,7 @@ function Base.show(io::IO, x::AdditiveDEAModel)
     disposX = x.disposX
     disposY = x.disposY
     dmunames = names(x)
-    
+
     eff = efficiency(x)
     slackX = slacks(x, :X)
     slackY = slacks(x, :Y)

diff --git a/src/deagdf.jl b/src/deagdf.jl
@@ -155,9 +155,8 @@ function deagdf(X::Matrix, Y::Matrix; alpha::Float64 = 0.5, rts::Symbol = :CRS,
                 # No contraint to add for constant returns to scale
             elseif rts == :VRS
                 @constraint(deamodel, sum(lambda) == 1)
-            else
-                error("Invalid returns to scale $rts. Returns to scale should be :CRS or :VRS")
             end
+            # No need for else statement as RTS parameter already checked in first stage
 
             # Optimize and return results
             JuMP.optimize!(deamodel)
@@ -169,7 +168,7 @@ function deagdf(X::Matrix, Y::Matrix; alpha::Float64 = 0.5, rts::Symbol = :CRS,
 
             # Check termination status
             if termination_status(deamodel) != MOI.LOCALLY_SOLVED
-                @warn ("DMU $i termination status: $(termination_status(deamodel)). Primal status: $(primal_status(deamodel)). Dual status: $(dual_status(deamodel))")
+                @warn ("DMU $i Slacks termination status: $(termination_status(deamodel)). Primal status: $(primal_status(deamodel)). Dual status: $(dual_status(deamodel))")
             end
 
         end

diff --git a/src/peers.jl b/src/peers.jl
@@ -80,10 +80,6 @@ struct DEAPeers <: AbstractDEAPeers
         dmunames = names(x)
         dmunamesref = String[]
 
-        if n != size(lambda,1)
-            error("Number of observation in the model and number of rows in peers matrix do not match")
-        end
-
         # Remove very small numbers
         droptol!(lambda, atol)
 

diff --git a/test/dea.jl b/test/dea.jl
@@ -2,7 +2,7 @@
 @testset "RadialDEAModel" begin
 
     ## Test Radial DEA Models with FLS Book data
-    # Test against results with R
+    # Test against results in R
     X = [5 13; 16 12; 16 26; 17 15; 18 14; 23 6; 25 10; 27 22; 37 14; 42 25; 5 17]
     Y = [12; 14; 25; 26; 8; 9; 27; 30; 31; 26; 12]
 
@@ -272,6 +272,16 @@
     @test slacks(deaooWeak, :X) ≈ [0; 0; 0; 0; 0] atol=1e-14
     @test slacks(deaooWeak, :Y) ≈ [0; 0; 0; 0; 0] atol=1e-14
 
+    # Test if weak disposability in inputs in the Input oriented model works
+    # In this example, same result as stron disposability
+    @test efficiency(dea(X, Y, orient = :Input, disposX = :Strong)) ==
+        efficiency(dea(X, Y, orient = :Input, disposX = :Weak))
+
+    # Test if weak disposability in outputs in the Output oriented model works
+    # In this example, same result as strong disposability
+    @test efficiency(dea(X, Y, orient = :Output, disposY = :Strong)) ==
+        efficiency(dea(X, Y, orient = :Output, disposY = :Weak))
+
     # ------------------
     # DMU names
     # ------------------
@@ -294,4 +304,27 @@
     @test occursin("Length of names greater than number of observations", string(logs))
     @test value == ["A", "B", "C", "D", "E"]
 
+    # ------------------
+    # Test Vector and Matrix inputs and outputs
+    # ------------------
+    # Tests against results in R
+
+    # Inputs is Matrix, Outputs is Vector
+    X = [2 2; 1 4; 4 1; 4 3; 5 5; 6 1; 2 5; 1.6	8]
+    Y = [1; 1; 1; 1; 1; 1; 1; 1]
+
+    @test efficiency(dea(X, Y, orient = :Input)) ≈ [1; 1; 1; 0.6; 0.4; 1; 0.6666666667; 0.625]
+
+    # Inputs is Vector, Output is Matrix
+    X = [1; 1; 1; 1; 1; 1; 1; 1]
+    Y = [7 7; 4 8; 8 4; 3 5; 3 3; 8 2; 6 4; 1.5 5]
+
+    @test efficiency(dea(X, Y, orient = :Output)) ≈ [1; 1; 1; 1.555555556; 2.333333333; 1; 1.272727273; 1.6]
+
+    # Inputs is Vector, Output is Vector
+    X = [2; 4; 8; 12; 6; 14; 14; 9.412]
+    Y = [1; 5; 8; 9; 3; 7; 9; 2.353]
+
+    @test efficiency(dea(X, Y, orient = :Input)) ≈ [0.4; 1; 0.8; 0.6; 0.4; 0.4; 0.5142857143; 0.2]
+
 end
diff --git a/test/deaadd.jl b/test/deaadd.jl
@@ -2,7 +2,7 @@
 @testset "AdditiveDEAModel" begin
 
     ## Test Additive DEA Models with FLS Book data
-    # Test against results with R
+    # Test against results in R
     X = [5 13; 16 12; 16 26; 17 15; 18 14; 23 6; 25 10; 27 22; 37 14; 42 25; 5 17]
     Y = [12; 14; 25; 26; 8; 9; 27; 30; 31; 26; 12]
 
@@ -464,14 +464,15 @@
     @test_throws ErrorException deaadd([1; 2; 3], [4; 5; 6], :Ones, rts = :Error) # Invalid returns to scale
     @test_throws ErrorException deaadd([1; 2; 3], [4; 5; 6], :Error) # Invalid model
     @test_throws ErrorException deaadd([1; 2; 3], [4; 5; 6], :Ones, orient = :Error) # Invalid orientation
+    @test_throws ErrorException deaadd([1; 2; 3], [4; 5; 6], :Custom, rhoX = [1; 1; 1], rhoY = [1; 1; 1], orient = :Error) # Invalid orientation with custom weights
 
     @test_throws ErrorException deaadd([1; 2 ; 3], [4 ; 5]) #  Different number of observations
     @test_throws ErrorException deaadd([1; 2], [4 ; 5], Xref = [1; 2; 3; 4]) # Different number of observations in reference sets
     @test_throws ErrorException deaadd([1 1; 2 2], [4 4; 5 5], Xref = [1 1 1; 2 2 2]) # Different number of inputs
     @test_throws ErrorException deaadd([1 1; 2 2], [4 4; 5 5], Yref = [4 4 4; 5 5 5]) # Different number of inputs
     @test_throws ErrorException deaadd([1; 2; 3], [4; 5; 6], rts = :Error) # Invalid returns to scale
-    @test_throws ErrorException deaadd([1; 2; 3], [4; 5; 6], :Custom, rhoX = [1; 2; 3; 4]) # Different size of weights
-    @test_throws ErrorException deaadd([1; 2; 3], [4; 5; 6], :Custom, rhoY = [4; 5; 6; 7]) # Different size of weights
+    @test_throws ErrorException deaadd([1; 2; 3], [4; 5; 6], :Custom, rhoX = [1; 2; 3; 4], rhoY = [1; 1; 1]) # Different size of weights
+    @test_throws ErrorException deaadd([1; 2; 3], [4; 5; 6], :Custom, rhoX = [1; 1; 1], rhoY = [4; 5; 6; 7]) # Different size of weights
 
     @test_throws ErrorException deaadd([1; 2; 3], [4; 5; 6], :Ones, rhoX = [1; 1; 1]) # Weights not allowed if model != :Custom
 
@@ -525,4 +526,42 @@
 
     @test_throws ErrorException deaadd([1; 2 ; 3], [4 ; 5; 6], orient = :Output, disposY = :Weak)  # Weak output disposability not possible in output oriented model
 
+    # ------------------
+    # Test Vector and Matrix inputs and outputs
+    # ------------------
+    # Tests against results in R
+
+    # Inputs is Matrix, Outputs is Vector
+    X = [2 2; 1 4; 4 1; 4 3; 5 5; 6 1; 2 5; 1.6	8]
+    Y = [1; 1; 1; 1; 1; 1; 1; 1]
+
+    @test efficiency(deaadd(X, Y, :Ones)) ≈ [0.0; 0.0; 0.0; 3.0; 6.0; 2.0; 3.0; 5.2]
+
+    # Inputs is Vector, Output is Matrix
+    X = [1; 1; 1; 1; 1; 1; 1; 1]
+    Y = [7 7; 4 8; 8 4; 3 5; 3 3; 8 2; 6 4; 1.5 5]
+
+    @test efficiency(deaadd(X, Y, :Ones)) ≈ [0.0; 0.0; 0.0; 6.0; 8.0; 2.0; 4.0; 7.5]
+
+    # Inputs is Vector, Output is Vector
+    X = [2; 4; 8; 12; 6; 14; 14; 9.412]
+    Y = [1; 5; 8; 9; 3; 7; 9; 2.353]
+
+    @test efficiency(deaadd(X, Y, :Ones)) ≈ [0; 0; 0; 0;4; 7.333333333; 2; 8.059]
+
+    # ------------------
+    # RAM and BAM with orientation
+    # ------------------
+
+    X = [1; 2; 3; 2; 4]
+    Y = [2; 3; 4; 1; 3]
+
+    @test efficiency(deaadd(X, Y, :RAM, orient = :Graph)) ≈ [0.0; 0.0; 0.0; 1/3; 1/3]
+    @test efficiency(deaadd(X, Y, :RAM, orient = :Input)) ≈ [0.0; 0.0; 0.0; 1/3; 2/3]
+    @test efficiency(deaadd(X, Y, :RAM, orient = :Output)) ≈ [0.0; 0.0; 0.0; 2/3; 1/3]
+
+    @test efficiency(deaadd(X, Y, :BAM, orient = :Graph)) ≈ [0.0; 0.0; 0.0; 2/3; 2/3]
+    @test efficiency(deaadd(X, Y, :BAM, orient = :Input)) ≈ [0.0; 0.0; 0.0; 1; 2/3]
+    @test efficiency(deaadd(X, Y, :BAM, orient = :Output)) ≈ [0.0; 0.0; 0.0; 2/3; 1]
+
 end
diff --git a/test/deacost.jl b/test/deacost.jl
@@ -17,7 +17,7 @@
 
 
     ## Test Cost DEA Model with Zofío and Prieto (2006) data.
-    # Test agains results with R
+    # Test agains results in R
     X = [5 3; 2 4; 4 2; 4 8; 7 9]
     Y = [7 4; 10 8; 8 10; 5 4; 3 6]
     W = [2 1; 2 1; 2 1; 2 1; 2 1.0]
@@ -70,6 +70,7 @@
     @test_throws ErrorException deacost([1; 2; 3], [4; 5; 6], [1; 2; 3], rts = :Error) # Invalid returns to scale
     @test_throws ErrorException deacost([1; 2; 3], [4; 5; 6], [1; 2; 3; 4]) # Different number of observation in prices
     @test_throws ErrorException deacost([1 1; 2 2; 3 3 ], [4; 5; 6], [1 1 1; 2 2 2; 3 3 3]) # Different number of input prices and inputs
+    @test_throws ErrorException deacost([1; 2; 3], [4; 5; 6], [1; 2; 3], dispos = :Error) # Invalid disposability
 
     # ------------------
     # Weak Disposability Tests
@@ -89,4 +90,30 @@
     @test efficiency(deacostWeak, :Technical) ≈ [1.0; 1.0; 1.0; 1.0; 0.5]
     @test efficiency(deacostWeak, :Allocative) ≈ [1.0; 1.0; 1.0; 1.0; 1.0]
 
+    # ------------------
+    # Test Vector and Matrix inputs and outputs
+    # ------------------
+    # Tests against results in R
+
+    # Inputs is Matrix, Outputs is Vector
+    X = [2 2; 1 4; 4 1; 4 3; 5 5; 6 1; 2 5; 1.6	8]
+    Y = [1; 1; 1; 1; 1; 1; 1; 1]
+    W = [1 1; 1 1; 1 1; 1 1; 1 1; 1 1; 1 1; 1 1]
+
+    @test efficiency(deacost(X, Y, W)) ≈ [1; 0.8; 0.8; 0.5714285714; 0.4; 0.5714285714; 0.5714285714; 0.4166666667]
+
+    # Inputs is Vector, Output is Matrix
+    X = [1; 1; 1; 1; 1; 1; 1; 1]
+    Y = [7 7; 4 8; 8 4; 3 5; 3 3; 8 2; 6 4; 1.5 5]
+    W = [1; 1; 1; 1; 1; 1; 1; 1]
+
+    @test efficiency(deacost(X, Y, W)) ≈ [1; 1; 1; 1; 1; 1; 1; 1]
+
+    # Inputs is Vector, Output is Vector
+    X = [2; 4; 8; 12; 6; 14; 14; 9.412]
+    Y = [1; 5; 8; 9; 3; 7; 9; 2.353]
+    W = [1; 1; 1; 1; 1; 1; 1; 1]
+
+    @test efficiency(deacost(X, Y, W)) ≈ [1; 1; 1; 1; 0.5; 0.4761904762; 0.8571428571; 0.2843710157]
+
 end
diff --git a/test/deaddf.jl b/test/deaddf.jl
@@ -2,7 +2,7 @@
 @testset "DirectionalDEAModel" begin
 
     ## Test Directional DF DEA Models with FLS Book data
-    # Test against results with R
+    # Test against results in R
     X = [5 13; 16 12; 16 26; 17 15; 18 14; 23 6; 25 10; 27 22; 37 14; 42 25; 5 17]
     Y = [12; 14; 25; 26; 8; 9; 27; 30; 31; 26; 12]
 
@@ -326,4 +326,27 @@
     @test_throws ErrorException deaddf([1 1; 2 2; 3 3], [4; 5; 6], [1 1 1; 2 2 2; 3 3 3], [4; 5; 6]) # Different size of inputs direction
     @test_throws ErrorException deaddf([1; 2; 3], [4 4; 5 5; 6 6], [1; 2; 3], [4 4 4; 5 5 5; 6 6 6]) # Different size of inputs direction
 
+    # ------------------
+    # Test Vector and Matrix inputs and outputs
+    # ------------------
+    # Tests against results in R
+
+    # Inputs is Matrix, Outputs is Vector
+    X = [2 2; 1 4; 4 1; 4 3; 5 5; 6 1; 2 5; 1.6	8]
+    Y = [1; 1; 1; 1; 1; 1; 1; 1]
+
+    @test efficiency(deaddf(X, Y, X, Y)) ≈ [0; 0; 0; 0.25; 0.4285714286; 0; 0.2; 0.2307692308]
+
+    # Inputs is Vector, Output is Matrix
+    X = [1; 1; 1; 1; 1; 1; 1; 1]
+    Y = [7 7; 4 8; 8 4; 3 5; 3 3; 8 2; 6 4; 1.5 5]
+
+    @test efficiency(deaddf(X, Y, X, Y)) ≈ [0; 0; 0; 0.2173913043; 0.4; 0; 0.12; 0.2307692308]
+
+    # Inputs is Vector, Output is Vector
+    X = [2; 4; 8; 12; 6; 14; 14; 9.412]
+    Y = [1; 5; 8; 9; 3; 7; 9; 2.353]
+
+    @test efficiency(deaddf(X, Y, X, Y)) ≈ [0.4285714286; 0; 0.1111111111; 0.25; 0.4285714286; 0.4285714286; 0.3207547170; 0.6666666667]
+
 end
diff --git a/test/deagdf.jl b/test/deagdf.jl
@@ -92,4 +92,26 @@
     @test_throws ErrorException deagdf([1 1; 2 2], [4 4; 5 5], alpha = 0.5, Yref = [4 4 4; 5 5 5]) # Different number of inputs
     @test_throws ErrorException deagdf([1; 2; 3], [4; 5; 6], alpha = 0.5, rts = :Error) # Invalid returns to scale
 
+    # ------------------
+    # Test Vector and Matrix inputs and outputs
+    # ------------------
+
+    # Inputs is Matrix, Outputs is Vector
+    X = [2 2; 1 4; 4 1; 4 3; 5 5; 6 1; 2 5; 1.6	8]
+    Y = [1; 1; 1; 1; 1; 1; 1; 1]
+
+    @test efficiency(deagdf(X, Y, rts = :VRS, slack = false)) ≈ [1; 1; 1; 1; 1; 1; 1; 1] atol = 1e-5
+
+    # Inputs is Vector, Output is Matrix
+    X = [1; 1; 1; 1; 1; 1; 1; 1]
+    Y = [7 7; 4 8; 8 4; 3 5; 3 3; 8 2; 6 4; 1.5 5]
+
+    @test efficiency(deagdf(X, Y, rts = :VRS, slack = false)) ≈ [1; 1; 1; 1; 1; 1; 1; 1] atol = 1e-5
+
+    # Inputs is Vector, Output is Vector
+    X = [2; 4; 8; 12; 6; 14; 14; 9.412]
+    Y = [1; 5; 8; 9; 3; 7; 9; 2.353]
+
+    @test efficiency(deagdf(X, Y, rts = :VRS)) ≈ [1; 1; 1; 1; 0.410097; 0.634489; 1; 0.20504] atol = 1e-5
+
 end
diff --git a/test/deaprofit.jl b/test/deaprofit.jl
@@ -31,4 +31,33 @@
     @test_throws ErrorException deaprofit([1 1; 2 2; 3 3], [4; 5; 6], [1 1; 2 2; 3 3], [4; 5; 6], [1 1 1; 2 2 2; 3 3 3], [4; 5; 6]) # Different size of inputs direction
     @test_throws ErrorException deaprofit([1; 2; 3], [4 4; 5 5; 6 6], [1; 2; 3], [4 4; 5 5; 6 6], [1; 2; 3], [4 4 4; 5 5 5; 6 6 6]) # Different size of inputs direction
 
+    # ------------------
+    # Test Vector and Matrix inputs and outputs
+    # ------------------
+    # Tests against results in R
+
+    # Inputs is Matrix, Outputs is Vector
+    X = [2 2; 1 4; 4 1; 4 3; 5 5; 6 1; 2 5; 1.6	8]
+    Y = [1; 1; 1; 1; 1; 1; 1; 1]
+    W = [1 1; 1 1; 1 1; 1 1; 1 1; 1 1; 1 1; 1 1]
+    P = [1; 1; 1; 1; 1; 1; 1; 1]
+
+    @test efficiency(deaprofit(X, Y, W, P, X, Y)) ≈ [0; 0.1666666667; 0.1666666667; 0.375; 0.5454545455; 0.375; 0.375; 0.5283018868]
+
+    # Inputs is Vector, Output is Matrix
+    X = [1; 1; 1; 1; 1; 1; 1; 1]
+    Y = [7 7; 4 8; 8 4; 3 5; 3 3; 8 2; 6 4; 1.5 5]
+    W = [1; 1; 1; 1; 1; 1; 1; 1]
+    P = [1 1; 1 1; 1 1; 1 1; 1 1; 1 1; 1 1; 1 1]
+
+    @test efficiency(deaprofit(X, Y, W, P, X, Y)) ≈ [0; 0.1538461538; 0.1538461538; 0.6666666667; 1.142857143; 0.3636363636; 0.3636363636; 1]
+
+    # Inputs is Vector, Output is Vector
+    X = [2; 4; 8; 12; 6; 14; 14; 9.412]
+    Y = [1; 5; 8; 9; 3; 7; 9; 2.353]
+    W = [1; 1; 1; 1; 1; 1; 1; 1]
+    P = [1; 1; 1; 1; 1; 1; 1; 1]
+
+    @test efficiency(deaprofit(X, Y, W, P, X, Y)) ≈ [0.6666666667; 0; 0.0625; 0.1904761905; 0.4444444444; 0.3809523810; 0.2608695652; 0.6849978751]
+
 end
diff --git a/test/deaprofitability.jl b/test/deaprofitability.jl
@@ -54,4 +54,32 @@
     @test_throws ErrorException deaprofitability([1; 2; 3], [4 4; 5 5; 6 6], [1; 2; 3], [4 4 4; 5 5 5; 6 6 6]) # Different number of oputput prices and outputs
     @test_throws ErrorException efficiency(deaprofbl, :Error)
 
+    # ------------------
+    # Test Vector and Matrix inputs and outputs
+    # ------------------
+
+    # Inputs is Matrix, Outputs is Vector
+    X = [2 2; 1 4; 4 1; 4 3; 5 5; 6 1; 2 5; 1.6	8]
+    Y = [1; 1; 1; 1; 1; 1; 1; 1]
+    W = [1 1; 1 1; 1 1; 1 1; 1 1; 1 1; 1 1; 1 1]
+    P = [1; 1; 1; 1; 1; 1; 1; 1]
+
+    @test efficiency(deaprofitability(X, Y, W, P)) ≈ ( sum(Y .* P, dims = 2) ./ sum(X .* W, dims = 2) )  / 0.25 atol = 1e-5
+
+    # Inputs is Vector, Output is Matrix
+    X = [1; 1; 1; 1; 1; 1; 1; 1]
+    Y = [7 7; 4 8; 8 4; 3 5; 3 3; 8 2; 6 4; 1.5 5]
+    W = [1; 1; 1; 1; 1; 1; 1; 1]
+    P = [1 1; 1 1; 1 1; 1 1; 1 1; 1 1; 1 1; 1 1]
+
+    @test efficiency(deaprofitability(X, Y, W, P)) ≈ ( sum(Y .* P, dims = 2) ./ sum(X .* W, dims = 2) )  / 14 atol = 1e-5
+
+    # Inputs is Vector, Output is Vector
+    X = [2; 4; 8; 12; 6; 14; 14; 9.412]
+    Y = [1; 5; 8; 9; 3; 7; 9; 2.353]
+    W = [1; 1; 1; 1; 1; 1; 1; 1]
+    P = [1; 1; 1; 1; 1; 1; 1; 1]
+
+    @test efficiency(deaprofitability(X, Y, W, P)) ≈ ( sum(Y .* P, dims = 2) ./ sum(X .* W, dims = 2) )  / 1.25 atol = 1e-5
+
 end
diff --git a/test/dearevenue.jl b/test/dearevenue.jl
@@ -17,7 +17,7 @@
 
 
     ## Test Revenue DEA Model with Zofío and Prieto (2006) data.
-    # Test agains results with R
+    # Test agains results in R
     X = [5 3; 2 4; 4 2; 4 8; 7 9]
     Y = [7 4; 10 8; 8 10; 5 4; 3 6]
     P = [3 2; 3 2; 3 2; 3 2; 3 2.0]
@@ -70,6 +70,7 @@
     @test_throws ErrorException dearevenue([1; 2; 3], [4; 5; 6], [1; 2; 3], rts = :Error) # Invalid returns to scale
     @test_throws ErrorException dearevenue([1; 2; 3], [4; 5; 6], [1; 2; 3; 4]) # Different number of observation in prices
     @test_throws ErrorException dearevenue([1; 2; 3], [4 4; 5 5; 6 6], [4 4 4; 5 5 5; 6 6 6]) # Different number of output prices and outputs
+    @test_throws ErrorException dearevenue([1; 2; 3], [4; 5; 6], [1; 2; 3], dispos = :Error) # Invalid disposability
 
     # ------------------
     # Weak Disposability Tests
@@ -89,4 +90,30 @@
     @test efficiency(dearevenueWeak, :Technical) ≈ [1.0; 1.0; 1.0; 0.3333333333333333; 1.0]
     @test efficiency(dearevenueWeak, :Allocative) ≈ [1.0; 1.0; 1.0; 1.0; 1.0]
 
+    # ------------------
+    # Test Vector and Matrix inputs and outputs
+    # ------------------
+    # Tests against results in R
+
+    # Inputs is Matrix, Outputs is Vector
+    X = [2 2; 1 4; 4 1; 4 3; 5 5; 6 1; 2 5; 1.6	8]
+    Y = [1; 1; 1; 1; 1; 1; 1; 1]
+    P = [1; 1; 1; 1; 1; 1; 1; 1]
+
+    @test efficiency(dearevenue(X, Y, P)) ≈ [1; 1; 1; 1; 1; 1; 1; 1]
+
+    # Inputs is Vector, Output is Matrix
+    X = [1; 1; 1; 1; 1; 1; 1; 1]
+    Y = [7 7; 4 8; 8 4; 3 5; 3 3; 8 2; 6 4; 1.5 5]
+    P = [1 1; 1 1; 1 1; 1 1; 1 1; 1 1; 1 1; 1 1]
+
+    @test efficiency(dearevenue(X, Y, P)) ≈ [1; 0.8571428571; 0.8571428571; 0.5714285714; 0.4285714286; 0.7142857143; 0.7142857143; 0.4642857143]
+
+    # Inputs is Vector, Output is Vector
+    X = [2; 4; 8; 12; 6; 14; 14; 9.412]
+    Y = [1; 5; 8; 9; 3; 7; 9; 2.353]
+    P = [1; 1; 1; 1; 1; 1; 1; 1]
+
+    @test efficiency(dearevenue(X, Y, P)) ≈ [1; 1; 1; 1; 0.4615384615; 0.7777777778; 1; 0.2816951993]
+
 end