【数理考古学】とある実数列(Real Sequance)の規定例③オイラーの等式(Eulerian identity)e^πi=-1が意味するもの？

(2021年2月27日)Rによる投稿
(2021年4月08日)texとpythonによる追記

とりあえず、以下までに実数族に下属する(部分集合を構成する)整数族とネイピア族の振る舞いをまとめてきました。それぞれが実数全体を規定する戦略を備えている辺りが不穏ですね。
[【数理考古学】とある実数列の規定例①等差数列から加法整数群へ](https://qiita.com/ochimusha01/items/07ae45b915bb1aa6c9f5)
[【数理考古学】とある実数列の規定例②等比数列から乗法群へ](https://qiita.com/ochimusha01/items/58afced6ea8bebf391c7)
今回は趣を変え、それらが置かれている座標系全体の成立過程について考察を深めて行きたいと思います。

#０次元(<font color=blue>Zero Dimension</font>)
(それ自体は観測対象となり得ない)**観測原点**(Observation Origin)$0$と**観測限界およびその向こう側**(Observation Limit)$\infty$の間に「<font color=red>何もない</font>」状態から出発します。**群論**(Group Theory)でいう**空環概念**(Empty/Nullary Ring Concept)に対応し、その**元**(Element)は空となります。
[【初心者向け】群論概念②空環から整数環へ。](https://qiita.com/ochimusha01/items/da766832c87a604bcb8a)

>**空環概念**(Empty/Nullary Ring Concept)は「**空和概念**(Empty/Nullary Sum Concept)」の規約と「**空積概念**(Empty/Nullary Product Concept)」の規約の統合結果である。

①ここでいう「**空和概念**(Empty/Nullary Sum Concept)**の規約**」とは、**群論**(Group Theory)上における、**空集合**(Empty/Nullary Set)$\emptyset_n(n=0)=()$について「**0個の因子を足し合わせた演算結果集合**(Operation Result Set)**で、加法単位元**(Additive Identity)0**に等しい**(():=0)」とする規約をいう。

* **二項演算**(Binary Operation)のうち、**加減算**(Addition and Subtraction)すなわち**加法単位元**(Additive Identity)$0+Α+Β(Α+(-Β))$においては任意の元(Any Element)$Α_n(n=0)()=Β_n(n=0)()=Γ_n(n=0)()=0$なので**結合則**(Associative Law)0+(Α+Β)+C+…=0+Α+(Β+C)+…も**交換則**(Commutative Law)0+Α+Β+…=0+Β+Α+…も成立し、その演算結果は閉じている(演算の都度0=()に戻る)。

②ここでいう「**空積概念**(Empty/Nullary Product Concept)**の規約**」とは、**群論**(Group Theory)における、**空集合概念**(Empty/Nullary Set)$\emptyset _n(n=0)=()$について「**0個の因子を掛け合わせた演算結果集合**(Operation Result Set)**で、乗法単位元**(Multiplicative Identity)1**に等しい**(():=1)」とする規約をいう。

* **二項演算**(Binary Operation)のうち、**乗除算**(Multiplying and Dividing)すなわち**乗法単位元**(Multiplicative Identity)$1×Α×Β(Α×(1/Β))$においては任意の元(Any Element)$Α_n(n=0)()=Β_n(n=0)()=Γ_n(n=0)()=1$となるので**結合則**(Associative Law)1×(Α×Β)×C+…=1×Α×(Β×C)×…も交換則(Commutative Law)1×Α×Β×…=1×Β×Α×…も成立し、その演算結果は閉じている(演算の都度1=()に戻る)。
* **二項演算**(Binary Operation)のうち、**冪算**(Exponentiation)それ自体こそ分配法則や交換法則を満たさないものの、その計算範囲自体は上掲の演算に内包され、自明の理としてその演算結果も同じ範囲で閉じている($…A^2=A×A,A^2=A,A^0=\frac{A}{A}=1,A^-1=\frac{1}{A},A^-2=\frac{1}{A×A},…・$)。

③この時点ではあくまで**空集合**$\emptyset _n(n=0)()=(加減算における単位元)0=(乗除算における単位元)1$と規約されるだけで$\emptyset _n(n=2)=(0,1)$が成立する訳ではない点に注意されたい。実際、**単元集合**(Singleton)/**単位集合**(Unit Set)$Α_n(n=1)=(0)$及び$Β_n(n=1)=(1)$、さらには二つを合算した集合$Γ_n(n=2)=(0,1)$全てが**空環**/**空群**$\emptyset(n=0)=()$と異なりそれ自体は群も環も構成しない。

* exp(-1)(0.3678794)を求める**ベルヌーイの式**(1-1/n)^nにおけるn=1(すなわち(1-1/1)^1=0)の場合に該当する。未開封のくじ引きセットの様なもので当たりなら必ず入っており「**外れる確率**」を計算しても0にしかならない。いわば観測自体成立していない、すなわち「<font color=red>閉じた目蓋の向こう側に全てが待っている</font>」状態。実数概念など成立しようがないから、この時点では全てが虚数領域に留まっているともいえる。

Rによる実装例
![image.png](https://camo.qiitausercontent.com/50b27a422f44a58437118c2e374d33e434d131a7/68747470733a2f2f71696974612d757365722d636f6e74656e74732e696d6769782e6e65742f687474707325334125324625324671696974612d696d6167652d73746f72652e73332e61702d6e6f727468656173742d312e616d617a6f6e6177732e636f6d2532463025324635393737373225324632343161396436302d333830632d643731652d303865332d3038356435393665633866652e6769663f69786c69623d72622d312e322e32266175746f3d666f726d6174266769662d713d363026713d373526773d31343030266669743d6d617826733d3038666230623165366336333930323230343830633262373564656134326539)

```R
#FN01=Find the number of Napier 01

FN01<-function(NoS){
pie.data<-rep(1, NoS)
names0<-NULL

for(i in 1:NoS){
names0<-c(names0,paste0(i,"/",NoS))
}
names(pie.data)<-names0
pie(pie.data, col = heat.colors(NoS),main = "Probability of losing  the Antipodal point")
}

#アニメーション
library("animation")
Time_Code=c(1,1,1,2,2,2,3,3,3,4,4,4,5,5,5,6,6,6,7,7,8,8,9,9,10,10,11,12,13,14,15,16)
saveGIF({
for (i in Time_Code){
FN01(i)
}
}, interval = 0.1, movie.name = "FN01.gif")
```
しかしもちろん、ただ「<font color=red>瞼を開いた</font>」だけで**観測原点**(Observation Origin)0と**観測限界**(Observation Limit)Inf(inity)の狭間に何かを見出せる訳ではありません。上掲の様にこの段階における**観測結果集合**はあくまで空集合、すなわち各観測結果の関係が全て**互いに素**なので次元抽出も相対設定も行えず、ある意味**観測原点**を中心に円弧状、もしくは球面状に均等に広がっているのみ。

* この状態はある意味、統計学における「<font color=red>統計対象に値しない</font>」**名義尺度**(Nominal scale, Stanley Smith Stevens,1946年)、あるいは仮説検定において棄却される**帰無仮説**、すなわち「**ここで証明したい事は証明不可能である**」とする論拠に対応する。まさしく「(<font color=blue>数理処理可能な有意味水準を超えるデータは</font>)<font color=red>存在しない</font>」状態の数理化という訳である。<br>[確率変数とは](https://bellcurve.jp/statistics/blog/14006.html)<br>[確率変数の尺度](https://bellcurve.jp/statistics/course/1562.html)<br>[帰無仮説(null hypothesis)](https://bellcurve.jp/statistics/glossary/899.html)

[【数理考古学】とある実数列の規定例①等差数列から加法整数群へ](https://qiita.com/ochimusha01/items/07ae45b915bb1aa6c9f5)
![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fqiita-user-contents.imgix.net%2Fhttps%253A%252F%252Fcdn-ak.f.st-hatena.com%252Fimages%252Ffotolife%252Fo%252Fochimusha01%252F20201111%252F20201111061314.gif%3Fixlib%3Drb-1.2.2%26auto%3Dformat%26gif-q%3D60%26q%3D75%26w%3D1400%26fit%3Dmax%26s%3D0ec5ad367a9238d99fae4325c3a95b29?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=61fff111f1cc073c26829b0b7b27df91)
![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fqiita-user-contents.imgix.net%2Fhttps%253A%252F%252Fcdn-ak.f.st-hatena.com%252Fimages%252Ffotolife%252Fo%252Fochimusha01%252F20201130%252F20201130003431.gif%3Fixlib%3Drb-1.2.2%26auto%3Dformat%26gif-q%3D60%26q%3D75%26w%3D1400%26fit%3Dmax%26s%3Dfa2431533fe04a79971e43e74315f0c5?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=8771e19e25b95ea04d07b6d2c5afa583)

まるで禅問答の様ですが、この概念を幾何学的イメージに置き換えたのが以下。**有向グラフ**(Directed graph)で表すと「**観測者側からの観測対象に対する一方的観察**」と図示される点が重要。 

**Rによる表現例**
<font color=blue>これを実装した当時はまだ双方向で考えてました。</font>

![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/2fdc8660-4121-caac-1053-d76ee279495c.png)

```R
library(igraph)
par(mfrow=c(1,2),mar=c(3,3,3,3))
gd<-graph(c(
"Zero(0)","Inf(Inity)",
"Inf(Inity)","Zero(0)"
))
plot(gd,layout=layout.fruchterman.reingold)
```
**pythonによる表現例**
[【Python演算処理】NetworkXライブラリを導入する。](https://qiita.com/ochimusha01/items/802453df0b97dfcb0984)

###空積概念(<font color=blue>Empty/Nullary Product Concept</font>)
![tmpja7jjltg.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/59b30270-549c-bc34-c297-46af8a8826e5.png)

* **観測限界**を点と置いた場合、(0~$\infty$の範囲に渡る)**観測原点**との距離$x$は、何倍にしても何分の1にしても値は変わらない。
* ただし例外的に$\frac{x}{x}=1$となり、これが**乗法単位元**(Multiplicative 
Identity)の概念の出発点になる。

###空和概念(<font color=blue>Empty/Nullary Sum Concept</font>)
![tmpog2lnpsb.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/13e592de-0d0e-6217-525d-d59f5c6f65c7.png)

* **観測限界**を(0~$\infty$の範囲に渡る)**観測原点**との距離$x$を半径とする**円弧**/**球面**と置いた場合、それをどんなに足し合わせても値は変わらない。
* ただし例外的に$x-x=0$となり、これが**加法単位元**(Additive Identity Identity)の概念の出発点になる。

まるである種の哲学論、というより般若心境における「**色即是空色即是空**」の境地? 科学的にはむしろ、実験しながら「**有意差出ろ!! 有意差出ろ!!**」と念じ続ける心境と重なります。つまりどちらかというと(乳幼児が生きているとされる、無意識と自我と超自我が未分離状態の)**呪術的世界観**に近いとも。そういえば機械学習の根幹部たる(人間の神経を模して活性化関数によって0と1の境界線をハンドリングする)**パーセプトロン**（Perceptron）概念もまた、この辺りの境界認識問題に関わってくるのです。
[パーセプトロン-Wikipedia](https://ja.wikipedia.org/wiki/%E3%83%91%E3%83%BC%E3%82%BB%E3%83%97%E3%83%88%E3%83%AD%E3%83%B3)

>心理学者・計算機科学者のフランク・ローゼンブラットが1957年に考案し、1958年に論文を発表した。モデルは同じく1958年に発表されたロジスティック回帰と等価である。
>
>視覚と脳の機能をモデル化したものであり、パターン認識を行う。ただし学習については自明ではなく、特に多層パーセプトロンの機械学習に関する歴史は、それがパーセプトロンの歴史だと言っても過言ではない。
>
>1960年代に爆発的なニューラルネットブームを巻き起こしたが、60年代末のミンスキーらによる、単層パーセプトロンは線形分離可能なものしか学習できないという指摘は、多層パーセプトロンの学習が当時まだよくわからなかったことから、一時研究を停滞させた。
>
>影響を受けた変種といえるニューラルネットワークも多数提案されているが、それらについてはここでは略す。パーセプトロン自体はその後、ボルツマンマシンや、多層パーセプトロンの機械学習をある程度実用化したバックプロパゲーションなどによって1980年代に再興隆した後、1990年代にはまた停滞した。
>
>21世紀に入った後の研究の進展などにより、より層を増やしたパーセプトロンでの機械学習を指すいわゆる「ディープラーニング」の実用化と、それによる画像認識などの成果により、2020年現在はブームにある。

#0.5次元(<font color=blue>0.5 Dimension</font>)

次元抽出や中心割り出しが可能そうな有望な観測結果集合が得られて初めて「**統計科学的**/**確率論的アプローチ**」が可能となります。

**Rによる実装例**
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/4abc44b6-6429-6ce5-f87b-72a4833bade9.png)

```R
library(igraph)
par(mfrow=c(1,2),mar=c(3,3,3,3))
gd<-graph(c(
"Zero(0)","Observation Result Set(x)",
"Observation Result Set(x)","Zero(0)",
"Observation Result Set(x)","Inf(Inity)",
"Inf(Inity)","Observation Result Set(x)"
))
plot(gd,layout=layout.fruchterman.reingold)
```
**pythonによる表現例**
[【Python演算処理】NetworkXライブラリを導入する。](https://qiita.com/ochimusha01/items/802453df0b97dfcb0984)

* 次元抽出や中心割り出しが可能そうな有望な**観測結果集合**(Observation Result Set)が得られて初めて「**統計科学的**/**確率論的アプローチ**」が可能となる。最初の評価軸であるから、当然それは**X軸**や**時間軸**(t軸)の形で与えられる。

![tmpsppviz17.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/599d8059-5cbb-5037-aee3-460378ee1963.png)

```python
import networkx as nx

G = nx.DiGraph()
G.add_edge("Observation Origin(0)","Elements(0<x<∞)")
G.add_edge("Elements(0<x<∞)","Observation Limit(∞)")
G.add_edge("Elements(0<x<∞)","Observation Origin(0)",)
G.add_edge("Observation Limit(∞)","Elements(0<x<∞)")
nx.nx_agraph.view_pygraphviz(G, prog='fdp') 
```
![tmp93bjyv9h.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/8260337b-5a14-f349-26bd-9fbee0fe8a4f.png)

```python
import networkx as nx

G = nx.DiGraph()
G.add_edge("Start(0)","Observation Result set(0<=t<=T)")
G.add_edge("Observation Result set(0<=t<=T)","End(T)")
G.add_edge("Observation Result set(0<=t<=T)","Start(0)")
G.add_edge("End(T)","Observation Result set(0<=t<=T)")
nx.nx_agraph.view_pygraphviz(G, prog='fdp') 
```
* 統計の対象となる実データには当然(観測不可能な)**両極限**は現れず、**最大値**(Max)と**最小値**(Min)がそれが推察可能な形で得られるとも限らない。

![tmpjjlbvf7n.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/b218b558-f4a7-9c2d-b507-fa43f670c812.png)

```python
import networkx as nx

G = nx.DiGraph()
G.add_edge("Min","Elements(Min<=x<=Max)")
G.add_edge("Elements(Min<=x<=Max)","Max")
G.add_edge("Elements(Min<=x<=Max)","Min",)
G.add_edge("Max","Elements(Min<=x<=Max)")
nx.nx_agraph.view_pygraphviz(G, prog='fdp') 
```
* この辺りのニュアンスも含め、**確率変数の尺度**(Stanley Smith Stevens,1946年)では「<font color=red>データ間の順位にのみ意味がある</font>」**順位尺度**(Ordinal scale)から「<font color=red>データ間の差にのみ意味がある</font>」**間隔尺度**(Interval scale)にかけてが該当する。<br>[確率変数の尺度](https://bellcurve.jp/statistics/course/1562.html)

この段階での中心的関心は「**最大値と最小値の狭間に**(平均や中央値や最頻値から)**中心を見出す**」事にあります。原則としてサンプル数nが多いほど**信頼度**(Reliability)が増えていきます。

Rによる実装例 **N=10**
![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20210117%2F20210117003848.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=37ccbadc6c100a2d89a877ae0590a7fd)

```R
ND<-function(x){
cx<-rnorm(10)
cy<-rnorm(10)
plot(cx,cy,xlim=c(-4,4),asp=1,ylim=c(-4,4),main="Normal Distribution",xlab="",ylab="")
c0<-seq(-pi,pi,length=61)
cx0<-cos(c0)
cy0<-sin(c0)
par(new=T)#上書き指定
plot(cx0+mean(cx),cy0+mean(cy),type="l",asp=1,xlim=c(-4,4),ylim=c(-4,4),main="",xlab="",ylab="",col=rgb(0,1,0))
par(new=T)#上書き指定
plot(cx0*2+mean(cx),cy0*2+mean(cy),type="l",asp=1,xlim=c(-4,4),ylim=c(-4,4),main="",xlab="",ylab="",col=rgb(0,1,0))
par(new=T)#上書き指定
plot(cx0*3+mean(cx),cy0*3+mean(cy),type="l",asp=1,xlim=c(-4,4),ylim=c(-4,4),main="",xlab="",ylab="",col=rgb(0,1,0))
text(mean(cx),mean(cy),"X",col=rgb(0,1,0))
}
#アニメーション
library("animation")
Time_Code=seq(1,16,length=16)
saveGIF({
for (i in Time_Code){
　ND(i)
}
}, interval = 0.1, movie.name = "ND04.gif")
```
Rによる実装例 **N=100**
![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20210117%2F20210117003225.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=38fcd483a692861fb49da6e992c5761a)

```R
ND<-function(x){
cx<-rnorm(100)
cy<-rnorm(100)
plot(cx,cy,xlim=c(-4,4),asp=1,ylim=c(-4,4),main="Normal Distribution",xlab="",ylab="")
c0<-seq(-pi,pi,length=61)
cx0<-cos(c0)
cy0<-sin(c0)
par(new=T)#上書き指定
plot(cx0+mean(cx),cy0+mean(cy),type="l",asp=1,xlim=c(-4,4),ylim=c(-4,4),main="",xlab="",ylab="",col=rgb(0,1,0))
par(new=T)#上書き指定
plot(cx0*2+mean(cx),cy0*2+mean(cy),type="l",asp=1,xlim=c(-4,4),ylim=c(-4,4),main="",xlab="",ylab="",col=rgb(0,1,0))
par(new=T)#上書き指定
plot(cx0*3+mean(cx),cy0*3+mean(cy),type="l",asp=1,xlim=c(-4,4),ylim=c(-4,4),main="",xlab="",ylab="",col=rgb(0,1,0))
text(mean(cx),mean(cy),"X",col=rgb(0,1,0))
}
#アニメーション
library("animation")
Time_Code=seq(1,16,length=16)
saveGIF({
for (i in Time_Code){
　ND(i)
}
}, interval = 0.1, movie.name = "ND03.gif")
```
Rによる実装例 **N=1000**
![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20210117%2F20210117002845.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=3fdd3f6ce5d5ffdccc5284ed07799871)

```R
ND<-function(x){
cx<-rnorm(1000)
cy<-rnorm(1000)
plot(cx,cy,xlim=c(-4,4),asp=1,ylim=c(-4,4),main="Normal Distribution",xlab="",ylab="")
c0<-seq(-pi,pi,length=61)
cx0<-cos(c0)
cy0<-sin(c0)
par(new=T)#上書き指定
plot(cx0+mean(cx),cy0+mean(cy),type="l",asp=1,xlim=c(-4,4),ylim=c(-4,4),main="",xlab="",ylab="",col=rgb(0,1,0))
par(new=T)#上書き指定
plot(cx0*2+mean(cx),cy0*2+mean(cy),type="l",asp=1,xlim=c(-4,4),ylim=c(-4,4),main="",xlab="",ylab="",col=rgb(0,1,0))
par(new=T)#上書き指定
plot(cx0*3+mean(cx),cy0*3+mean(cy),type="l",asp=1,xlim=c(-4,4),ylim=c(-4,4),main="",xlab="",ylab="",col=rgb(0,1,0))
text(mean(cx),mean(cy),"X",col=rgb(0,1,0))
}
#アニメーション
library("animation")
Time_Code=seq(1,16,length=16)
saveGIF({
for (i in Time_Code){
　ND(i)
}
}, interval = 0.1, movie.name = "ND02.gif")
```
Rによる実装例 **三次元表示**
![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20201110%2F20201110010727.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=90db2db66ddc81e814913c6046ffcf04)

```R
library(rgl)
Rtime<-seq(0,2*pi,length=60)
cx<-rnorm(1000)
cy<-rnorm(1000)
cz<-rnorm(1000)
plot3d(cx,cy,cz,type="p",xlim=c(-4,4),ylim=c(-4,4),zlim=c(-4,4))
movie3d(spin3d(axis=c(0,0,1),rpm=5),duration=10,fps=25,movie="~/Desktop/r06")
```
①**平均値**(Mean)に立脚する**中心検知技法**(Center Detecting Technique)は、**外れ値**(Outliers)の影響を受けやすいので、それを除去する技術が共に発展してきた。
[【初心者向け】誤差関数(ERF)と相補誤差関数 (ERFC)。](https://qiita.com/ochimusha01/items/ee8b32e4b2cc7ce4acb0)

②**中央値**(Median)は**平均値**(Mean)より外れ値に対して**ロバスト**(頑強)であり、かつ富が富裕層に集中している状況の様な中央値と平均値が解離している状況での中央検知に威力を発揮する。
[【初心者向け】平均と標準偏差、中央値と平均偏差](https://qiita.com/ochimusha01/items/971d62ebc5e4cac6a61b)

③正直、**最頻値**(Mode)に立脚する中央検知の意義はあまり研究されてないが、上掲二つのアプローチと併記して語られる事が多い。**ピアソンの経験則**によれば、ヒストグラムの分布が左右対称の場合は中央で重なる平均値・中央値・最頻値が、分布が左に偏る(右の裾野が広がる)場合には「**最頻値・中央値・平均値**」、分布が右に偏る(左の裾野が広がる)場合には「**平均値・中央値・最頻値**」の順に並ぶ事が多いという。
[【初心者向け】度数分布/度数分布図と最頻値](https://qiita.com/ochimusha01/items/ea7ca95bf05510004bc6)
[【プログラマーのための統計学】平均値・中央値・最頻値](https://qiita.com/Esfahan/items/5ca15a334d22007e10ef)
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/1f1275b1-8c3d-4532-bcf8-82f810e172bf.png)
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/e073c58e-b399-64b5-8d74-1a6894ce37b7.png)
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/7c6c3462-4bf0-9c68-e877-7ddd7b83c3cf.png)

④さらにややこしい事に、現代における「**統計科学的**/**確率論的アプローチ**」は、かかる中央検知問題ばかりに拘泥している訳でもない。**MCMC法**の広まりや**機械学習概念**の台頭はそういう観点に関連してくる。
[【初心者向け】「モンテカルロ法」概説](https://qiita.com/ochimusha01/items/979f43806ef00910ac81)
![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fqiita-image-store.s3.ap-northeast-1.amazonaws.com%2F0%2F597772%2F300dc4cc-97c6-dc12-9003-56d360995a03.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&s=8a5d33c311c3d5e8c77e8cb3f136f48c)
[【初心者向け】ロジスティック方程式とその関連範囲](https://qiita.com/ochimusha01/items/2f0118ecda2e17d366ba)

⑤**確率楕円**(Probability Ellipse)概念は、こうした考え方にさらに線形代数における固有値や固有ベクトル概念を加えた内容で、因子分析(Factor Analysis)の基礎。**主成分分析**(PCA=Principal Component Analysis)に至っては、同一平面上に展開する**観測結果集合**を複数の観測楕円に分割したりする。
**Rによる楕円関数の実装例①**
![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20201218%2F20201218180039.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=8b4eb55139403aefd9499ad4b2756f91)

```R
Cos_Waltz01<-function(Radius){
#格子描画
grx=seq(-3,3,length=4)
gry01=rep(c(-1,1),2)
gry02=rep(c(1,-1),2)
plot(grx,gry01,type="l",asp=1,xlim=c(-1,1),ylim=c(-1,1),main="Cosine Waltz",xlab="x=-Cos(θ)→Sin(θ)→Cos(θ)→Sin(θ)→-Cos(θ)",ylab="Cos(θ)",col=rgb(0,0,0))
par(new=T)
plot(grx,gry02,type="l",asp=1,xlim=c(-1,1),ylim=c(-1,1),main="",xlab="",ylab="",col=rgb(0,0,0))
#円描画
c0<-seq(0,pi*2,length=61)
cy<-cos(c0)
cx00<-cos(c0-pi*Radius)
cx01<-cos(c0-pi*Radius)-2
cx02<-cos(pi*Radius-c0)+2
par(new=T)
plot(cx00,cy,type="l",asp=1,xlim=c(-1,1),ylim=c(-1,1),main="",xlab="",ylab="",col=rgb(0,1,0))
for(i in 1:61)
{
segments(0,0,cx00[i],cy[i],col=rgb(0,1,0))
}
#補助線
abline(h=0,col=rgb(1,0,0))
abline(v=0,col=rgb(1,0,0))
}
#アニメーション
library("animation")
Time_Code=seq(0, 4,length=61)
saveGIF({
for (i in Time_Code[1:60]){
　Cos_Waltz01(i)
}
}, interval = 0.1, movie.name = "Cos_Waltz04.gif")
```
**Rによる楕円関数の実装例②**
![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20201218%2F20201218063808.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=8cb56725775f434df555872b0d640157)

```R
Cos_Waltz01<-function(Radius){
#格子描画
grx=seq(-3,3,length=4)
gry01=rep(c(-1,1),2)
gry02=rep(c(1,-1),2)
plot(grx,gry01,type="l",xlim=c(-3,3),ylim=c(-1,1),main="Cosine Waltz",xlab="x=-Cos(θ)→Sin(θ)→Cos(θ)→Sin(θ)→-Cos(θ)",ylab="Cos(θ)",col=rgb(0,0,0))
par(new=T)
plot(grx,gry02,type="l",xlim=c(-3,3),ylim=c(-1,1),main="",xlab="",ylab="",col=rgb(0,0,0))
#円描画
c0<-seq(0,pi*2,length=61)
cy<-cos(c0)
cx00<-cos(c0-pi*Radius)
cx01<-cos(c0-pi*Radius)-2
cx02<-cos(pi*Radius-c0)+2
par(new=T)
plot(cx00,cy,type="l",xlim=c(-3,3),ylim=c(-1,1),main="",xlab="",ylab="",col=rgb(0,1,0))
par(new=T)
plot(cx01,cy,type="l",xlim=c(-3,3),ylim=c(-1,1),main="",xlab="",ylab="",col=rgb(1,0,0))
par(new=T)
plot(cx02,cy,type="l",xlim=c(-3,3),ylim=c(-1,1),main="",xlab="",ylab="",col=rgb(0,0,1))
for(i in 1:61)
{
segments(0,0,cx00[i],cy[i],col=rgb(0,1,0))
}
for(i in 1:61)
{
segments(-2,0,cx01[i],cy[i],col=rgb(1,0,0))
}
for(i in 1:61)
{
segments(2,0,cx02[i],cy[i],col=rgb(0,0,1))
}
#補助線
abline(h=0,col=rgb(1,0,0))
abline(v=0,col=rgb(1,0,0))
}
#アニメーション
library("animation")
Time_Code=seq(0, 4,length=61)
saveGIF({
for (i in Time_Code[1:60]){
　Cos_Waltz01(i)
}
}, interval = 0.1, movie.name = "Cos_Waltz03.gif")
```
[Pythonで綴る多変量解析　4-1. 因子分析(scikit-learn)](https://qiita.com/y_itoh/items/227cb33317ceb09199c2)
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/76621615-8fcf-0b39-48fa-2d3d53368b37.png)
[主成分分析を Python で理解する](https://qiita.com/maskot1977/items/082557fcda78c4cdb41f)
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/72db9731-90e7-e8c0-074c-0c1e41321c08.png)

⑥適用する**距離関数**(Matric Function)の見直しは、かかる評価空間全体に影響を与える可能性がある。
[【初心者向け】方形描画関数②距離空間との関係。](https://qiita.com/ochimusha01/items/06e49a35dffbd12b35aa)
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/a630af03-e3b2-7801-7c58-6cc510bfd2fb.png)
[符号なし距離関数で作る魚眼レンズ効果](https://qiita.com/7CIT/items/6c931012d94184888ef7)
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/67710cd8-4e31-9d1a-e391-3948f467f8ea.png)
この投稿の目的は(こうした諸概念が現れる以前の)近世から近代にかけての数学者の足跡を追う事なので、とりあえずこのまま先を続けます。いずれにせよ興味深いのは、**統計学の世界**は**実軸**(Real Axis)を追加するごとに起こる**視界狭窄**(Vision Tunnel)と無縁な事で、これは統計学の座標系が原則として全て**虚軸**(Imaginal Axis)に分類される事を意味している可能性を示唆しているのです。

<font color=purple>0次元</font>を（乳幼児が生きているとされる、無意識と意識と超自我の境界線が存在しない)**呪術的世界**に擬えるなら、<font color=purple>0.5次元</font>は(ハイハイを始め、つかまり立ちを覚え、自分の足で歩き始める)**経験主義的世界**に該当します。哲学史の世界においては、人間が**白紙状態**(Tabula rasa)から経験の蓄積を始めると考える**経験主義的立場**と、経験知に対して**先天的/生得的/先験的に**(a priori)実存する認識が存在すると主張する**観念論的立場**が長年絶対的対立を続けてきましたが、**虚数**(Imaginal)概念が存在する数理の世界では、両者をこの様に連続的に統合的に扱う事が可能となるのですね。実際発達心理学によれば、人間の認知能力の成長は概ねこういう具合に進むとされています。
[タブラ・ラサ | 現代美術用語辞典ver.2.0 - Artscape](https://artscape.jp/artword/index.php/%E3%82%BF%E3%83%96%E3%83%A9%E3%83%BB%E3%83%A9%E3%82%B5)
[アプリオリ - Wikipedia](https://ja.wikipedia.org/wiki/%E3%82%A2%E3%83%97%E3%83%AA%E3%82%AA%E3%83%AA)

#1次元(<font color=blue>1.0 Dimension</font>)
**中心検出**に成功したら、いよいよこれを**単位元**(Identity Element)として**反元**(Inverse Elements)を設定し、**群**(Group)を定義する事が可能となります。**確率変数の尺度**(Stanley Smith Stevens,1946年)では**比例尺度**(Proportional Scale)に該当。これ以降、**関数列**(Function Sequance)を含む**演算結果集合**(Operation Result Set)が、**観測結果集合**(Ovsevation Result Set)と矛盾しない限り、その**族**(Famiry)を代表する形となります。
[確率変数の尺度](https://bellcurve.jp/statistics/course/1562.html)

**Rによる表現例**
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/0947d8fd-e780-0323-d7ab-9765240e9f71.png)

```R
library(igraph)
gd<-graph(edges=c("Zero(0)","Element(g)","Element(g)","Zero(0)","Element(g)","Inf(inity)","Inf(inity)","Element(g)","Zero(0)","Inverse(-g)","Inverse(-g)","Zero(0)","Inverse(-g)","-Inf(inity)","-Inf(inity)","Inverse(-g)"))
plot(gd,layout=layout.fruchterman.reingold) 
```
**pythonによる表現例**
[【Python演算処理】NetworkXライブラリを導入する。](https://qiita.com/ochimusha01/items/802453df0b97dfcb0984)

###加法実数群(Additive Real Group)の場合
[【数理考古学】とある実数列の規定例①等差数列から加法整数群へ](https://qiita.com/ochimusha01/items/07ae45b915bb1aa6c9f5)
**10進法**(Decimal System)、すなわち10を法とする**加法整数群**(Additive Integer Group)の**剰余類**(Residue Class)の多層添字化によって理論上、小数点下を含む実数全体を扱える様に拡張された**両側無限数列**(Both-sides Infinite Sequance)。
![tmpi_zdfl10.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/88b6634a-2cd4-36c3-9ede-b9445925c5e8.png)

##乗法実数群(Multiplicative Real Group)の場合
[【数理考古学】とある実数列の規定例②等比数列から乗法群へ](https://qiita.com/ochimusha01/items/58afced6ea8bebf391c7)
上掲の**加法実数群**(均等尺)を$-\infty→0,0→1,+\infty→\infty$に射影した**片側無限数列**(One-sides Infinite Sequance)。均等尺上の乗除算が加減算、冪乗算が乗算に置き換えられる。
[【数理考古学】常用対数表を使った計算](https://qiita.com/ochimusha01/items/b7ce3c7d86ab39263c68)
![tmpapj3itug.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/13fd7425-4d3c-932f-b0f3-e3491ec0bfaf.png)

おやおや？　正負の概念が導入され、**無限遠点**$\infty$が**無限大**$+\infty$と**無限小**$-\infty$に引き裂かれると同時に、**観測原点**(Observation Origin)が座標上から消えてしまいました。このタイミングで**偶奇**(Parity)の概念が導入されるのは決して偶然でなく、どうやら「**観測原点と観測対象の関係**」もまた**偶数系**(Even System)と**奇数系**(Odd System)に引き裂かれてしまう様なのです。哲学でいう「(独我論や無我論と対比される)**他我論**」の登場ですね。それを数理の世界は「<font color=purple>観測原点0が半径1の円弧を計測単位とする観測系を樹立すると、同時に観測対象たる1を観測原点とし、観測原点0を観測対象とする観測系を得る。前者を偶数系、後者を奇数系という</font>」といった具合に説明し、この説明を用いると独我論と他我論と無我論の境界線が相対化によって消失してしまう訳です。
[他我 - Wikipedia](https://ja.wikipedia.org/wiki/%E4%BB%96%E6%88%91)

#トーラス座標系(<font color=blue>Torus Coordinate System</font>)と円筒座標系(<font color=blue>Cylinder Coordinate System</font>)の狭間

以前の投稿ではこの「**独我**」状態から…
![image.png](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20201111%2F20201111061314.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=0ec5ad367a9238d99fae4325c3a95b29)
![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20201130%2F20201130003431.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=fa2431533fe04a79971e43e74315f0c5)
この「**我他相対**」状態に直接飛躍すると説明してきました。
[【数理考古学】とある実数列の規定例①等差数列から加法整数群へ](https://qiita.com/ochimusha01/items/07ae45b915bb1aa6c9f5)
![image.gif](https://camo.qiitausercontent.com/9e756e40dceca25c93c9517bc61e134658c8ae2e/68747470733a2f2f63646e2d616b2e662e73742d686174656e612e636f6d2f696d616765732f666f746f6c6966652f6f2f6f6368696d7573686130312f32303230313132362f32303230313132363032303830302e676966)
![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fqiita-user-contents.imgix.net%2Fhttps%253A%252F%252Fcdn-ak.f.st-hatena.com%252Fimages%252Ffotolife%252Fo%252Fochimusha01%252F20210105%252F20210105003122.gif%3Fixlib%3Drb-1.2.2%26auto%3Dformat%26gif-q%3D60%26q%3D75%26w%3D1400%26fit%3Dmax%26s%3D1c5e7641bb15c183f387fb9c6c92759a?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&s=8b3c1ac01fe5a4de036898e40dd0d907)
![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fqiita-user-contents.imgix.net%2Fhttps%253A%252F%252Fcdn-ak.f.st-hatena.com%252Fimages%252Ffotolife%252Fo%252Fochimusha01%252F20210105%252F20210105071334.gif%3Fixlib%3Drb-1.2.2%26auto%3Dformat%26gif-q%3D60%26q%3D75%26w%3D1400%26fit%3Dmax%26s%3D6cc2d838c056f87959b9d35678b0d60a?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&w=1400&fit=max&s=302dce6d5da2947a562c26895875c993)

しかし実は「**観測原点0が半径1の円弧を計測単位とする観測系を樹立した段階**」からさらにこれに直交する実軸を考える場合、そこに現れるべきは**円筒**(Unit Cylinder)でなくドーナツ型の**トーラス**(単数形:torus, 複数形:tori)なのです。

###Rによる実装例
![image.gif](https://cdn-ak.f.st-hatena.com/images/fotolife/o/ochimusha01/20210220/20210220215504.gif)
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/dec9e8fd-fa03-b7ad-dfa1-a1cf9f927839.png)

```R
library(rgl)
radians<-seq(0,60*pi,length=3000)
f0<-function(x){exp(x*(0+1i))}
f1<-function(x){sqrt(1-x^2)}
r0<-f0(radians)
c0<-seq(-1,1,length=1500)
#c1<-c(f1(c0)/2,-1*f1(c0)/2)
c1<-c(f1(c0),-1*f1(c0))
r1<-r0*(1+c1)
Real<-Re(r1)
Imaginal<-Im(r1)
Rtime<-c(seq(0,1,length=1500),seq(1,0,length=1500))
plot3d(Real,Imaginal,Rtime,type="l",xlim=c(-2,2),ylim=c(-2,2),zlim=c(0,2),xlab="",ylab="",zlab="")

#大円描写
Real01<-Re(f0(radians))
Imaginal01<-Im(f0(radians))
Rtime01<-seq(1/2,1/2,length=3000)
points3d(Real01,Imaginal01,Rtime01,col=rgb(0,0,1))

#小円描写
Real01<-Re(f0(radians))+1
Imaginal01<-Im(f0(radians))/2+1/2
Rtime01<-seq(0,0,length=3000)
points3d(Rtime01,Real01,Imaginal01,col=rgb(1,0,0))

movie3d(spin3d(axis=c(0,0,1),rpm=5),duration=10,fps=25,movie="~/Desktop/unit16")
```
###pythonによる実装例
[【Python演算処理】単位球面を巡る数理②そして単位トーラスへ](https://qiita.com/ochimusha01/items/c0f8c21e19d15f3e3569)
![image.gif](https://cdn-ak.f.st-hatena.com/images/fotolife/o/ochimusha01/20210412/20210412193831.gif)
![image.gif](https://cdn-ak.f.st-hatena.com/images/fotolife/o/ochimusha01/20210412/20210412193633.gif)
![image.gif](https://cdn-ak.f.st-hatena.com/images/fotolife/o/ochimusha01/20210412/20210412193720.gif)
トーラスは座標系(Coordinate System)としては媒介変数t,p(0≦t≦2π,0≦p≦2π)を用いて以下の様に表します。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/597772/6c7c7add-cf52-77b4-4ccb-569619feac1a.png)

* x=大半径R×cos(t)+小半径r×cos(p)×cos(t)
* y=大半径R×cos(t)+小半径r×cos(p)×sin(t)
* z=小半径r×sin(p)

トーラスではなくなる極限状態は以下。

* 大半径R=0,小半径R=1の時、**単位球面**(Unit Sphere)を二重に描く。従って**経緯度法**(経度-180度~180度に対して緯度-90度~90度)や**3次元極座標系**(水平角φ=-π~πに対して垂直角0~π)は適用範囲を半分に減らしてなお球面上の全座標を示せるのはこれによる。
* 逆に大半径R=1,小半径R=0の時、Z=0となって**ただの単位円**(Unit Circle)になる。

その中間過程に現れる大半径R=1,小半径R=1のトーラスを**単位トーラス**(Unit Torus)と呼ぶ事にしましょう。この単位トーラスをアイロンを掛けた様に平らに伸ばして円筒形座標系に移行するにはどうしたらいいでしょう？

* まずはトーラスを円筒状に展開する。この時高さは2π。中心線で分割するとπ×2。
* これを湾曲させて対蹠間の距離を2(片側の対蹠の座標が1なら-1)に到達させる。

おや、この操作には見覚えが…

```math
i^2=e^{πi}=-1
```
Rによる実装例
![image.gif](https://qiita-user-contents.imgix.net/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2Fo%2Fochimusha01%2F20210129%2F20210129205903.gif?ixlib=rb-1.2.2&auto=format&gif-q=60&q=75&s=60686f670bcfec0aa8598a240a67536d)

```R
Eulers_identity<-function(n){
rim01<-pi*complex(real=0,imaginary=1)
Tarm=c(complex(real=1,imaginary=0),seq(rim01/n,rim01/n,length=n))
Reduce(function(a,b) {a+a*b}, Tarm, accumulate = T)
}

#グラフを描画してみる。
Eulers_identity01_plot<-function(n){
ei01<-Eulers_identity(n)
ei02<-Conj(ei01)

#背景描画
c00<-seq(-pi,pi,length=61)
c01<-exp(c00*(0+1i))
plot(c01+(0+1i),type="l",asp=1,xlim=c(-pi,pi),ylim=c(0,pi),main="Even Scale to Logarithmic Scale",xlab="Real",ylab="Imaginal")
par(new=T) 
plot(c01+1,type="l",asp=1,xlim=c(-pi,pi),ylim=c(0,pi),main="",xlab="",ylab="")
par(new=T)
plot(c01-1,type="l",asp=1,xlim=c(-pi,pi),ylim=c(0,pi),main="",xlab="",ylab="")
par(new=T) # 上書き指定
plot(-Re(ei01)+1,Im(ei01),asp=1,xlim=c(-pi,pi),ylim=c(0,pi),type="l",main="", xlab="", ylab="",col=rgb(0,0,1))
par(new=T) # 上書き指定
plot(-Re(-ei02)-1,-Im(ei02),asp=1,xlim=c(-pi,pi),ylim=c(0,pi),type="l",main="", xlab="", ylab="",col=rgb(1,0,0))
segments(-2,0,2,0)
abline(h=(1+3/n)^n,col=c(200,200,200))
abline(h=(1+2.5/n)^n,col=c(200,200,200))
abline(h=(1+2/n)^n,col=c(200,200,200))
abline(h=(1+1.5/n)^n,col=c(200,200,200))
abline(h=(1+1/n)^n,col=rgb(0,1,0))
abline(h=(1+0.5/n)^n,,col=c(200,200,200))
abline(h=1,col=c(200,200,200))
abline(h=(1-0.5/n)^n,,col=c(200,200,200))
abline(h=(1-1/n)^n,col=rgb(0,1,0))
abline(h=(1-1.5/n)^n,col=c(200,200,200))
abline(h=(1-2/n)^n,col=c(200,200,200))
abline(h=(1-2.5/n)^n,col=c(200,200,200))
abline(h=(1-3/n)^n,col=c(200,200,200))
Text<-paste("(1±1/",n,")^",n,"=",round((1+1/n)^n,digits=4),"/",round((1-1/n)^n,digits=4))
legend("topleft", legend=c("π×(0+1i)","π×(0-1i)",Text), lty=c(1,1,1), col=c(rgb(0,0,1),rgb(1,0,0),rgb(0,1,0)))
}

#アニメーションさせてみる。

library("animation")
Time_Code=c(1,2,4,8,16,32,64,128)
saveGIF({
for (i in Time_Code){
　　Eulers_identity01_plot(i)
}
}, interval = 0.1, movie.name = "TEST.gif")
```
まだ関連数理を完全にはまとめ切れていませんが、基本方針としてはそう間違ってない気もしています。あくまで「補助線」の話に過ぎないので他の考え方も一杯出来そうだし。そんな感じで以下続報…