1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
|
"""Tests for module dr on Dimensionality Reduction """
# Author: Remi Flamary <remi.flamary@unice.fr>
# Minhui Huang <mhhuang@ucdavis.edu>
#
# License: MIT License
import numpy as np
import ot
import pytest
try: # test if autograd and pymanopt are installed
import ot.dr
nogo = False
except ImportError:
nogo = True
@pytest.mark.skipif(nogo, reason="Missing modules (autograd or pymanopt)")
def test_fda():
n_samples = 90 # nb samples in source and target datasets
np.random.seed(0)
# generate gaussian dataset
xs, ys = ot.datasets.make_data_classif('gaussrot', n_samples)
n_features_noise = 8
xs = np.hstack((xs, np.random.randn(n_samples, n_features_noise)))
p = 1
Pfda, projfda = ot.dr.fda(xs, ys, p)
projfda(xs)
np.testing.assert_allclose(np.sum(Pfda**2, 0), np.ones(p))
@pytest.mark.skipif(nogo, reason="Missing modules (autograd or pymanopt)")
def test_wda():
n_samples = 100 # nb samples in source and target datasets
np.random.seed(0)
# generate gaussian dataset
xs, ys = ot.datasets.make_data_classif('gaussrot', n_samples)
n_features_noise = 8
xs = np.hstack((xs, np.random.randn(n_samples, n_features_noise)))
p = 2
Pwda, projwda = ot.dr.wda(xs, ys, p, maxiter=10)
projwda(xs)
np.testing.assert_allclose(np.sum(Pwda**2, 0), np.ones(p))
@pytest.mark.skipif(nogo, reason="Missing modules (autograd or pymanopt)")
def test_wda_normalized():
n_samples = 100 # nb samples in source and target datasets
np.random.seed(0)
# generate gaussian dataset
xs, ys = ot.datasets.make_data_classif('gaussrot', n_samples)
n_features_noise = 8
xs = np.hstack((xs, np.random.randn(n_samples, n_features_noise)))
p = 2
P0 = np.random.randn(10, p)
P0 /= P0.sum(0, keepdims=True)
Pwda, projwda = ot.dr.wda(xs, ys, p, maxiter=10, P0=P0, normalize=True)
projwda(xs)
np.testing.assert_allclose(np.sum(Pwda**2, 0), np.ones(p))
@pytest.mark.skipif(nogo, reason="Missing modules (autograd or pymanopt)")
def test_prw():
d = 100 # Dimension
n = 100 # Number samples
k = 3 # Subspace dimension
dim = 3
def fragmented_hypercube(n, d, dim):
assert dim <= d
assert dim >= 1
assert dim == int(dim)
a = (1. / n) * np.ones(n)
b = (1. / n) * np.ones(n)
# First measure : uniform on the hypercube
X = np.random.uniform(-1, 1, size=(n, d))
# Second measure : fragmentation
tmp_y = np.random.uniform(-1, 1, size=(n, d))
Y = tmp_y + 2 * np.sign(tmp_y) * np.array(dim * [1] + (d - dim) * [0])
return a, b, X, Y
a, b, X, Y = fragmented_hypercube(n, d, dim)
tau = 0.002
reg = 0.2
pi, U = ot.dr.projection_robust_wasserstein(X, Y, a, b, tau, reg=reg, k=k, maxiter=1000, verbose=1)
U0 = np.random.randn(d, k)
U0, _ = np.linalg.qr(U0)
pi, U = ot.dr.projection_robust_wasserstein(X, Y, a, b, tau, U0=U0, reg=reg, k=k, maxiter=1000, verbose=1)
|
