ryamada22

---
title: "フィッシャーの正確確率検定とパーミュテーションテスト"
output: html_document
---

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
```

## フィッシャーの正確確率検定

ryamada22

---
title: "フォームを確かめる"
author: "ryamada"
date: "2021年9月10日"
output: html_document
---

```{r}
### packageのインストールがまだなら・・・
# install.packages(c("knitr","rgl","onion"))

ryamada22

---
title: "ラケット操作"
author: "ryamada"
date: "2021年9月4日"
output: html_document
---

#```{r setup, include=FALSE}
```{r}
library(knitr)

ryamada22

---
title: "Birkhoff polytopeと二重確率行列とその周辺"
author: "ryamada"
date: "2021年3月8日"
output: html_document
---

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
```

ryamada22

---
title: "Grassmann-Plucker Relation"
author: "ryamada"
date: "2021年3月3日"
output: html_document
---

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
```

ryamada22

library(randtests)

## ポアソン乱数を発生させる
n <- 1000
lambda <- 1.3
R <- rpois(n,lambda)
plot(R,type="h")
## ポアソン乱数列に対して設計したわけではない！runs.test()で評価してみる
out <- runs.test(R)
out # この結果が妥当だとは思えない

ryamada22

---
title: "ホロサイクルを使って双曲上半平面測地距離の有限値を計算する"
author: "ryamada"
date: "2020年10月19日"
output: html_document
---

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
```

ryamada22

---
title: "双曲幾何と射影幾何変換"
author: "ryamada"
date: "2020年10月18日"
output: html_document
---

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
```

ryamada22

---
title: "双曲幾何と射影幾何変換"
author: "ryamada"
date: "2020年10月18日"
output: html_document
---

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
```

ryamada22

my.B.mut <- function(B,k){
  new.B <- B
  n <- length(B[,1])
  new.B[k,] <- - B[k,]
  new.B[,k] <- - B[,k]
  for(i in 1:n){
    for(j in 1:n){
      if(i != k){
        if(j != k){
          new.B[i,j] <- B[i,j] + sign(B[i,k]) * max(B[i,k]*B[k,j],0)