diff --git a/src/plotting.jl b/src/plotting.jl
index 405c1485..417be4d9 100644
--- a/src/plotting.jl
+++ b/src/plotting.jl
@@ -380,8 +380,8 @@ autocorplot
     title --> "Input-Output cross correlation"
     seriestype := :corrplot
     @series begin
-        label --> ""
-        lags .* Ts, xc
+        label --> reshape(["u$ui -> y$yi" for ui in 1:size(u,2), yi in 1:size(y,2)], 1, :)
+        lags .* Ts, reshape(xc, size(xc, 1), :)
     end
 end
 
diff --git a/src/subspace2.jl b/src/subspace2.jl
index d1e1c890..bbf6251f 100644
--- a/src/subspace2.jl
+++ b/src/subspace2.jl
@@ -529,6 +529,7 @@ function find_PK(L1,L2,Or,n,p,m,r,s1,s2,A,C)
     
     K0 = vl*pinv(hl)
     W = (vl - K0*hl)*(vl-K0*hl)'
+    W ./= (length(vl) - length(K0))
     
     Q = W[1:n,1:n] |> Hermitian
     S = W[1:n,n+1:n+p]