Contents
Summary
- 분석 스크립트
- 4개 요지(임강요지, 정암리와요지, 관현리와요지, 왕진리와요지)를 modeling 한 후 그외 소비지를 한 그룹씩 입력하여 각 입력 도편의 사후확률 계산하였음
- 사후확률 : 모델링 된 데이터에, 새로운 데이터가 추가되었을 때, 새로운 데이터가 어느쪽으로 분류될 확률이 큰지를 알려주는 확률
- 또한 4개의 요지를 대상으로 모델링된 판별 모형에 도편 샘플 프로파일을 입력으로 하여 판별하고 그 결과를 LD1, LD2 좌표에 산개하고, 각 그룹의 중심점에서 가장 멀리 있는 샘플과의 거리를 반지름으로 하는 원(실선)과 그 반지름의 1.2배 되는 원(점선)을 그려, 사후 확률만으로 판단하기 어려운 생산지-소비지 연관성을 판단하는 기준을 제시함
모델
- 네 개의 요지를 training하여 얻은 모델은 크게 관현리와요지와 그 외 세 요지, 두 그룹으로 나뉘어질 수 있다. 즉 (관현리와요지)와 (임강요지, 정암리와요지, 왕진리와요지)는 서로 교차하지 않지만, (임강요지, 정암리와요지, 왕진리와요지)의 내부에서는 상호 교차함을 확인하였다. 좀 더 세분하 했을 경우 관현리와요지와, 임강요지, 그리고 나머지 두 요지로 나뉘어 질 수 있다.
- 이를 통해 관현리 외 세 개 요지는 상대적으로 유사한 특성을 갖는다고 판단할 수 있음
임강사지
- 10개 임강사지 도편 중 7개 도편이 임강요지로 분류될 확률이 높게 나왔지만, 1-3, 1-6은 정암리와요지로 분류될 확률(각각 9.993681e-01, 9.163543e-01)이 높게 나왔다. 1-7의 경우 판별식의 사후 확률(1.000000e+00) 상으로 관현리와요지로 분류될 수 있지만, 관현리와요지 모델 자료와의 거리가 멀어 실제 생산지가 관현리와요지로 판정하기 보다는, (임강요지, 정암리와요지, 왕진리와요지)에서 생산되지 않았다는 판단 근거가 됨
능산리사지
- 5개 능산리사지 도편 중 이상치(3-5)를 제거하였을 때, 3-2는 왕진리와요지로, 3-4는 임강요지로 분류될 확률이 높았으나 3-1, 3-3의 도편은 어느 그룹으로 분류된다 판단하기 어려움
군수리사지
- 5개 군수리사지의 도편 중 4개의 도편이 정암리와요지로 분류될 확률이 높게 나왔다. 하지만 4-5 도편의 경우 확률(1.000000e+00) 상으로 왕진리와요지로 분류될 확률이 높게 나왔지만 2차원 평면상 왕진리와요지로 보기 어렵다고 보여짐
용정리사지
- 5개 용정리사지의 도편 중 5-2, 5-3은 정암리와요지로 분류될 확률이 높게 나왔다. 하지만 그외 5-1, 5-5, 5-4의 경우 관현리와요지나 임강요지로 분류될 확률이 높게 나왔으나, 2차원평면상 어느 그룹으로도 분류된다고 판단하기 어려움
금강사지
- 5개 금강사지의 도편 중 6-1과 6-3은 왕진리와요지로, 6-4는 임강요지로 분류될 확률이 높다고 분류되었으나, 6-2, 6-5의 경우 어느 그룹으로 분류된다고 판단하기 어려움
왕흥사지
- 5개 왕흥사지의 도편 중 7-3의 경우 높은 확률로 정암리와요지로 분류되었으나, 그외 4개의 도편은 어느 그룹으로 분류된다고 판단하기 어려움
관북리
- 5개 관북리 도편 중 8-5의 경우 높은 확률로 임강요지로 분류되었으나, 그 외 4개 도편은 어느 그룹으로 분류된다고 판단하기 어려움
청마산성
- 5개 청마산성 도편 중 12-1, 12-2, 12-4 의 경우 임강요지로 분류되었으며, 12-5는 관현리와요지로 분류될 확률이 높았다. 하지만 12-3의 경우 어느 그룹으로 분류된다고 판단하기 어려움
부소산성
- 5개 부소산성 도편 중 13-3의 경우 임강요지로 분류될 확률이 높았으나 그외의 도편은 어느 그룹으로 분류된다고 판단하기 어려움
분류코드 범례(NAA백제시대)
분류 |
지역 |
Code |
요지 |
임강요지 |
IGYJ |
정암리와요지 |
JA |
|
관현리와요지 |
GH |
|
왕진리와요지 |
WJ |
|
소비지 |
임강사지 |
IGSJ |
능산리사지 |
NS |
|
군수리사지 |
GS |
|
용정리사지 |
YJ |
|
금강사지 |
GG |
|
왕흥사지 |
WH |
|
관북리 |
GB |
|
청마산성 |
CM |
|
부소산성 |
BS |
Model validation(LOOCV)
- 4개 요지(임강요지, 정암리와요지, 관현리와요지, 왕진리와요지) 바탕으로한 모델의 검증정보
. |
GH |
IGYJ |
JA |
WJ |
total |
GH |
6 |
0 |
0 |
0 |
1.0 |
IGYJ |
0 |
10 |
0 |
0 |
1.0 |
JA |
0 |
1 |
3 |
1 |
0.6 |
WJ |
0 |
1 |
2 |
2 |
0.4 |
입력:임강사지(I-IGSJ)
- 도편별 사후확률
IGSJ |
GH |
IGYJ |
JA |
WJ |
1-1 |
2.568072e-30 |
9.999999e-01 |
5.583141e-08 |
3.499069e-18 |
1-2 |
2.676254e-46 |
1.000000e+00 |
1.249613e-27 |
1.105761e-42 |
1-3 |
9.856512e-18 |
6.319320e-04 |
9.993681e-01 |
4.842642e-12 |
1-4 |
1.426383e-34 |
1.000000e+00 |
5.546896e-16 |
2.650947e-29 |
1-5 |
2.932354e-06 |
9.999971e-01 |
5.113388e-13 |
9.763064e-39 |
1-6 |
6.131797e-36 |
5.833263e-11 |
9.163543e-01 |
8.364567e-02 |
1-7 |
1.000000e+00 |
1.376397e-12 |
7.029946e-15 |
7.966918e-43 |
1-8 |
4.110082e-62 |
1.000000e+00 |
1.301682e-21 |
1.107095e-24 |
1-9 |
7.304652e-42 |
1.000000e+00 |
1.132820e-14 |
1.962671e-23 |
1-10 |
1.017196e-91 |
1.000000e+00 |
3.185674e-34 |
1.344886e-30 |
- 10개 임강사지 도편 중 7개 도편이 임강요지로 분류될 확률이 높게 나왔지만, 1-3, 1-6은 정암리와요지로 분류될 확률(각각 9.993681e-01, 9.163543e-01)이 높게 나왔다. 1-7의 경우 판별식의 사후 확률(1.000000e+00) 상으로 관현리와요지로 분류될 수 있지만, 관현리와요지 모델 자료와의 거리가 멀어 실제 생산지가 관현리와요지로 판정하기 보다는, (임강요지, 정암리와요지, 왕진리와요지)에서 생산되지 않았다는 판단 근거가
입력:능산리사지(I-NS)
- 도편별 사후확률
NS |
GH |
IGYJ |
JA |
WJ |
3-1 |
1.682103e-74 |
1.00000e+00 |
3.784822e-43 |
3.910748e-54 |
3-2 |
1.015043e-58 |
9.81102e-06 |
1.350360e-06 |
9.999888e-01 |
3-3 |
4.291460e-32 |
1.00000e+00 |
1.692435e-26 |
5.253458e-48 |
3-4 |
3.370387e-69 |
1.00000e+00 |
1.729070e-19 |
1.442709e-17 |
3-5 |
0.000000e+00 |
1.76492e-229 |
8.526613e-244 |
1.000000e+00 |
도편별 사후확률(이상치(3-5) 제거) <= 확률 변화 없음
|
GH |
IGYJ |
JA |
WJ |
3-1 |
1.682103e-74 |
1.00000e+00 |
3.784822e-43 |
3.910748e-54 |
3-2 |
1.015043e-58 |
9.81102e-06 |
1.350360e-06 |
9.999888e-01 |
3-3 |
4.291460e-32 |
1.00000e+00 |
1.692435e-26 |
5.253458e-48 |
3-4 |
3.370387e-69 |
1.00000e+00 |
1.729070e-19 |
1.442709e-17 |
- 5개 능산리사지 도편 중 이상치(3-5)를 제거하였을 때, 3-2는 왕진리와요지로, 3-4는 임강요지로 분류될 확률이 높았으나 3-1, 3-3의 도편은 어느 그룹으로 분류된다 판단하기 어려움
입력:군수리사지(I-GS)
- 도편별 사후확률
GS |
GH |
IGYJ |
JA |
WJ |
4-1 |
6.489118e-23 |
4.052153e-19 |
9.999996e-01 |
4.212926e-07 |
4-2 |
4.289309e-14 |
2.738210e-23 |
1.000000e+00 |
6.139267e-11 |
4-3 |
3.137342e-12 |
2.111425e-08 |
1.000000e+00 |
3.778152e-14 |
4-4 |
3.032998e-06 |
8.686740e-25 |
9.999970e-01 |
1.089192e-14 |
4-5 |
4.613587e-74 |
1.372079e-41 |
1.424402e-14 |
1.000000e+00 |
- 5개 군수리사지의 도편 중 4개의 도편이 정암리와요지로 분류될 확률이 높게 나왔다. 하지만 4-5 도편의 경우 확률(1.000000e+00) 상으로 왕진리와요지로 분류될 확률이 높게 나왔지만 2차원 평면상 왕진리와요지로 보기 어렵다고 보여짐
입력:용정리사지(I-YJ)
- 도편별 사후확률
YJ |
GH |
IGYJ |
JA |
WJ |
5-1 |
1.000000e+00 |
3.155456e-79 |
9.332851e-72 |
4.763343e-127 |
5-2 |
7.109075e-15 |
1.041676e-05 |
9.999896e-01 |
3.221504e-13 |
5-3 |
6.758722e-21 |
4.200643e-10 |
1.000000e+00 |
1.696601e-09 |
5-4 |
2.911624e-36 |
6.585683e-01 |
3.407343e-01 |
6.974859e-04 |
5-5 |
1.000000e+00 |
6.056954e-76 |
2.516887e-26 |
3.578774e-52 |
- 5개 용정리사지의 도편 중 5-2, 5-3은 정암리와요지로 분류될 확률이 높게 나왔다. 하지만 그외 5-1, 5-5, 5-4의 경우 관현리와요지나 임강요지로 분류될 확률이 높게 나왔으나, 2차원평면상 어느 그룹으로도 분류된다고 판단하기 어려움
입력:금강사지(I-GG)
- 도편별 사후확률
GG |
GH |
IGYJ |
JA |
WJ |
6-1 |
1.900559e-74 |
7.682224e-21 |
8.695597e-13 |
1.000000e+00 |
6-2 |
9.687179e-01 |
1.029406e-09 |
3.128209e-02 |
2.858163e-22 |
6-3 |
1.805594e-50 |
1.206323e-16 |
4.469955e-05 |
9.999553e-01 |
6-4 |
5.181481e-56 |
1.000000e+00 |
3.795474e-14 |
2.170800e-15 |
6-5 |
1.000000e+00 |
2.118989e-17 |
7.727864e-19 |
1.162716e-48 |
- 5개 금강사지의 도편 중 6-1과 6-3은 왕진리와요지로, 6-4는 임강요지로 분류될 확률이 높다고 분류되었으나, 6-2, 6-5의 경우 어느 그룹으로 분류된다고 판단하기 어려움
입력:왕흥사지(I-WH)
- 도편별 사후확률
WH |
GH |
IGYJ |
JA |
WJ |
7-1 |
1.000000e+00 |
3.026025e-119 |
1.171310e-88 |
5.336721e-149 |
7-2 |
1.000000e+00 |
7.315740e-193 |
6.018608e-169 |
5.271439e-271 |
7-3 |
3.853874e-20 |
9.675660e-14 |
1.000000e+00 |
2.511778e-09 |
7-4 |
9.120242e-02 |
9.087976e-01 |
8.299071e-16 |
1.230180e-45 |
7-5 |
1.000000e+00 |
5.873369e-121 |
5.148101e-87 |
5.173026e-146 |
- 5개 왕흥사지의 도편 중 7-3의 경우 높은 확률로 정암리와요지로 분류되었으나, 그외 4개의 도편은 어느 그룹으로 분류된다고 판단하기 어려움
입력:관북리(I-GB)
- 도편별 사후확률
GB |
GH |
IGYJ |
JA |
WJ |
8-1 |
1.000000e+00 |
1.455143e-84 |
7.817029e-86 |
1.363338e-149 |
8-2 |
1.000000e+00 |
1.583998e-155 |
7.206241e-125 |
1.481645e-203 |
8-3 |
1.376013e-176 |
2.913358e-48 |
1.644904e-46 |
1.000000e+00 |
8-4 |
7.074081e-30 |
1.000000e+00 |
4.010745e-51 |
5.496573e-90 |
8-5 |
1.316304e-52 |
1.000000e+00 |
1.727974e-18 |
2.553359e-24 |
- 5개 관북리 도편 중 8-5의 경우 높은 확률로 임강요지로 분류되었으나, 그 외 4개 도편은 어느 그룹으로 분류된다고 판단하기 어려움
입력:청마산성(I-CM)
- 도편별 사후확률
CM |
GH |
IGYJ |
JA |
WJ |
12-1 |
1.623808e-39 |
1.000000e+00 |
2.252062e-23 |
4.635575e-39 |
12-2 |
1.441774e-88 |
1.000000e+00 |
9.274613e-40 |
2.219970e-41 |
12-3 |
1.000000e+00 |
3.003541e-82 |
2.160137e-96 |
1.837852e-167 |
12-4 |
2.472651e-104 |
1.000000e+00 |
1.891780e-31 |
1.333612e-19 |
12-5 |
1.000000e+00 |
2.281754e-26 |
1.359897e-22 |
1.588973e-53 |
- 5개 청마산성 도편 중 12-1, 12-2, 12-4 의 경우 임강요지로 분류되었으며, 12-5는 관현리와요지로 분류될 확률이 높았다. 하지만 12-3의 경우 어느 그룹으로 분류된다고 판단하기 어려움
입력:부소산성(I-BS)
- 도편별 사후확률
BS |
GH |
IGYJ |
JA |
WJ |
13-1 |
4.938566e-156 |
8.856529e-81 |
3.261853e-43 |
1.000000e+00 |
13-2 |
1.007194e-37 |
2.754828e-11 |
6.133356e-01 |
3.866644e-01 |
13-3 |
1.223781e-32 |
1.000000e+00 |
2.079933e-15 |
2.486656e-29 |
13-4 |
4.558013e-89 |
2.037414e-09 |
3.313335e-16 |
1.000000e+00 |
13-5 |
3.400925e-103 |
1.000000e+00 |
1.712590e-37 |
3.322508e-30 |
- 5개 부소산성 도편 중 13-3의 경우 임강요지로 분류될 확률이 높았으나 그외의 도편은 어느 그룹으로 분류된다고 판단하기 어려움
