기본 콘텐츠로 건너뛰기

[Linear Algebra] 직교적 투영(Orthogonal Projection)

[data analysis] 데이터 라벨링(Data Labeling)

데이터 라벨링(Data Labeling)

클래스 표시 행렬

한 변수의 값들은 각 값에 적합된 그룹(클래스)에 1을, 나머지 그룹(클래스)에 0을 할당하는 방식으로 0과 1로 나타낼 수 있습니다. 표 1은 0, 1, 2, 3, 5의 값을 갖는 데이터를 5개의 클래스로 구분하여 각 값에 대응하는 클래스에만 1을 할당하는 과정을 나타낸 것입니다. 벡터 구조인 데이터의 각 값은 왼쪽을 시작으로 0 부터 부여되는 위치값을 가집니다. 이 위치값을 인덱스(index)라고 하며 클래스의 수와 인덱스의 갯수를 동일하게 하여 같은 위치에 1을 할당한 것입니다.

표 1 1차원 데이터의 클래스 표시행렬
data = [0, 1, 2, 3, 5]
valueindexclass 0 class1 class2 class 3 class 4
00 1 0 0 0 0
11 0 1 0 0 0
22 0 0 1 0 0
33 0 0 0 1 0
54 0 0 00 1

표 1의 결과는 식 1과 같이 벡터 형태인 데이터를 2차원 행렬 형태로 전환됩니다. 이 결과는 클래스표시행렬이라 합니다.

$$\begin{bmatrix}1&2&3&4 \end{bmatrix}\Rightarrow \begin{bmatrix}1&0&0&0\\ 0&1&0&0 \\0&0&1&0\\0&0&0&1\end{bmatrix}$$ (식 1)

식 1의 데이터는 고유한 값들로 구성되어 있는 특수한 형태입니다. 다음 코드의 객체 x와 같이 중접된 값들을 가지는 경우 클래스를 생성하기 위해 고유한 값들만을 추출할 필요가 있습니다. 이러한 과정을 포함하여 클래스 표시행렬을 생성하기 위해 sklearn.preprocessing.LabelBinarizer() 클래스를 사용합니다. 고유값은 이 클래스 속성 .class_으로 확인할 수 있습니다. 이 결과는 오름차순으로 정렬된 상태로 반환됩니다.

sklearn.preprocessing.LabelBinarizer()
  • 라벨표시 행렬을 생성하는 클래스
  • fit(): 메소드로 데이터를 클래스에 적합시킴
  • transform(): 메소드로 이진화된 라벨표시 행렬을 반환
  • 속성 .class_: 오름차순으로 정렬된 목록의 라벨을 확인
  • inverse_transform(변환된 객체): 변환된 객체의 원시형태를 반환
x=[5, 0, 2, 3, 0, 1, 3, 5, 3, 3]
lb = preprocessing.LabelBinarizer()
lb.fit(x)
print(lb.classes_)
[0 1 2 3 5]

.transform() 메서드로 변환된 클래스 표시행렬을 확인합니다.

print(lb.transform(x))
[[0 0 0 0 1]
 [1 0 0 0 0]
 [0 0 1 0 0]
 [0 0 0 1 0]
 [1 0 0 0 0]
 [0 1 0 0 0]
 [0 0 0 1 0]
 [0 0 0 0 1]
 [0 0 0 1 0]
 [0 0 0 1 0]]

이 행렬은 데이터의 고유값인 0, 1, 2, 3, 5에 대응하는 각 클래스에 대해 표시한 것으로 이 고유값외의 다른 값에 대해 클래스를 적합시킬 수 없습니다. 이 경우는 모든 클래스에 0이 할당됩니다. 다음 코드의 인수 4가 이 경우에 해당하는 것으로 이 값에 대응하는 클래스는 없는 상태이므로 모두 0으로 표시됩니다.

print(lb.transform([2, 4]))
[[0 0 1 0 0]
 [0 0 0 0 0]]

다중 클래스 표시 행렬

클래스표시행렬의 경우 1차원 데이터에 대응하는 것이라면 다중 클래스 표시 행렬은 2차원으로 구성된 라벨 데이터를 이진화하는 행렬입니다. 변환은 행 단위로 이루어집니다. 예를 들어 데이터 전체의 고유값을 오름차순으로 정리하여 [1, 2, 6, 7]이라면 각 요소의 인덱스에 대응하는 값은 1 나머지는 0이 됩니다. 이 고유값에서 7의 경우 [0, 0, 0,1]이 됩니다. 그러므로 다중클래스 표시행렬의 열차원은 4가 됩니다. 식 2는 위의 고유값을 지닌 데이터인 경우에서의 다중클래스표시 행렬을 나타낸 것입니다. 이 결과와 같이 이진화는 행단위로 이루어집니다. 그러므로 다중클래스표시행렬의 모양은 데이터의 행수 × 고유값의 수가 됩니다.

$$\begin{bmatrix} 1&7\\2&6\end{bmatrix} \Rightarrow $\begin{bmatrix} 1&0&0&1\\ 0&1&1&0\end{bmatrix}$$ (식 2)

다중클래스 데이터는 다음 코드의 y와 같이 모든 리스트의 크기가 같을 필요는 없습니다. 이 변형은 MultiLabelBinarizer() 클래스를 적용합니다.

sklearn.preprocessing.MultiLabelBinarizer()
  • 다중클래스 표시 행렬을 생성
  • .fit() 메소드로 데이터에 적합
  • .transform() 메소드로 데이터로부터 다중클래스 표시행렬을 생성

다음은 5개의리스트를 가지고 6개의 고유값을 가지는 라벨 데이터입니다. 그러므로 클래스 표시 행렬의 차원은 5 × 6이 됩니다.

y = [[2, 3, 4, 6], [2], [0, 1, 3], [0, 1, 2, 3, 4], [0, 1, 2]]
mlb = preprocessing.MultiLabelBinarizer()
mlb.fit(y)
print(mlb.classes_)
[0 1 2 3 4 6]
print(mlb.transform(y))
[[0 0 1 1 1 1]
 [0 0 1 0 0 0]
 [1 1 0 1 0 0]
 [1 1 1 1 1 0]
 [1 1 1 0 0 0]]

이 블로그의 인기 게시물

유사변환과 대각화

내용 유사변환 유사행렬의 특성 대각화(Diagonalization) 유사변환(Similarity transformation) 유사변환 n×n 차원의 정방 행렬 A, B 그리고 가역 행렬 P 사이에 식 1의 관계가 성립하면 행렬 A와 B는 유사하다고 하며 이 변환을 유사 변환 (similarity transformation)이라고 합니다. $$\begin{equation}\tag{1} A = PBP^{-1} \Leftrightarrow P^{-1}AP = B \end{equation}$$ 식 1의 유사 변환은 다음과 같이 고유값을 적용하여 특성 방정식 형태로 정리할 수 있습니다. $$\begin{align} B - \lambda I &= P^{-1}AP – \lambda P^{-1}P\\ &= P^{-1}(AP – \lambda P)\\ &= P^{-1}(A - \lambda I)P \end{align}$$ 위 식의 행렬식은 다음과 같이 정리됩니다. $$\begin{align} &\begin{aligned}\textsf{det}(B - \lambda I ) & = \textsf{det}(P^{-1}(AP – \lambda P))\\ &= \textsf{det}(P^{-1}) \textsf{det}((A – \lambda I)) \textsf{det}(P)\\ &= \textsf{det}(P^{-1}) \textsf{det}(P) \textsf{det}((A – \lambda I))\\ &= \textsf{det}(A – \lambda I)\end{aligned}\\ &\begin{aligned}\because \; \textsf{det}(P^{-1}) \textsf{det}(P) &= \textsf{det}(P^{-1}P)\\ &= \t

[matplotlib] 히스토그램(Histogram)

히스토그램(Histogram) 히스토그램은 확률분포의 그래픽적인 표현이며 막대그래프의 종류입니다. 이 그래프가 확률분포와 관계가 있으므로 통계적 요소를 나타내기 위해 많이 사용됩니다. plt.hist(X, bins=10)함수를 사용합니다. x=np.random.randn(1000) plt.hist(x, 10) plt.show() 위 그래프의 y축은 각 구간에 해당하는 갯수이다. 빈도수 대신 확률밀도를 나타내기 위해서는 위 함수의 매개변수 normed=True로 조정하여 나타낼 수 있다. 또한 매개변수 bins의 인수를 숫자로 전달할 수 있지만 리스트 객체로 지정할 수 있다. 막대그래프의 경우와 마찬가지로 각 막대의 폭은 매개변수 width에 의해 조정된다. y=np.linspace(min(x)-1, max(x)+1, 10) y array([-4.48810153, -3.54351935, -2.59893717, -1.65435499, -0.70977282, 0.23480936, 1.17939154, 2.12397372, 3.0685559 , 4.01313807]) plt.hist(x, y, normed=True) plt.show()

R 미분과 적분

내용 expression 미분 2차 미분 mosaic를 사용한 미분 적분 미분과 적분 R에서의 미분과 적분 함수는 expression()함수에 의해 생성된 표현식을 대상으로 합니다. expression expression(문자, 또는 식) 이 표현식의 평가는 eval() 함수에 의해 실행됩니다. > ex1<-expression(1+0:9) > ex1 expression(1 + 0:9) > eval(ex1) [1] 1 2 3 4 5 6 7 8 9 10 > ex2<-expression(u, 2, u+0:9) > ex2 expression(u, 2, u + 0:9) > ex2[1] expression(u) > ex2[2] expression(2) > ex2[3] expression(u + 0:9) > u<-0.9 > eval(ex2[3]) [1] 0.9 1.9 2.9 3.9 4.9 5.9 6.9 7.9 8.9 9.9 미분 D(표현식, 미분 변수) 함수로 미분을 실행합니다. 이 함수의 표현식은 expression() 함수로 생성된 객체이며 미분 변수는 다음 식의 분모의 변수를 의미합니다. $$\frac{d}{d \text{변수}}\text{표현식}$$ 이 함수는 어떤 함수의 미분의 결과를 표현식으로 반환합니다. > D(expression(2*x^3), "x") 2 * (3 * x^2) > eq<-expression(log(x)) > eq expression(log(x)) > D(eq, "x") 1/x > eq2<-expression(a/(1+b*exp(-d*x))); eq2 expression(a/(1 + b * exp(-d * x))) > D(eq2, "x") a * (b * (exp(-d * x) * d))/(1 + b