- 底辺私立医大を卒業した医者って頭悪いよね？ Part17
763 ：卵の名無しさん[sage]：2020/10/18(日) 04:42:31.02 ID:PYky+T0u - options(digits=22)
par(bty='l')
par(mfrow=c(2,1))
f <- function(x) x^x
curve(f(x),-0.1,1)
D(expression(x^x),'x')
f1 <- function(x) x^(x - 1) * x + x^x * log(x)
curve(f1(x),-0.1,1) ; abline(h=0,lty=3)
uniroot(f1,c(0.01,1)) ; 1/exp(1)
f(1/exp(1)) ; exp(-1/exp(1))
| - 底辺私立医大を卒業した医者って頭悪いよね？ Part17
764 ：卵の名無しさん[sage]：2020/10/18(日) 04:44:21.95 ID:PYky+T0u - >>762
値に０や１があるとカイ二乗検定が使えないから、ベイズ統計で処理することになるね。
| - 底辺私立医大を卒業した医者って頭悪いよね？ Part17
766 ：卵の名無しさん[sage]：2020/10/18(日) 15:49:59.13 ID:PYky+T0u - "
時刻t=0に、黒板に数字0が書かれている。
各時刻t=1,2,3,...において、数字0,1,2を下図のように1つずつ横一列に書いていく。
どの数字が書かれるかは同様に確からしい。
（t=0のとき）0
（t=5のときの一例）011212
各時刻tにおいて、数字を1つ書き終えた後、以下の操作を一回行う。
【操作】
この時刻に書かれた数字が0であった場合、そこから最も近い位置の0を探す。
この2つの0で挟まれた数字列の中に1がある場合、それら全てを0に書き換える。
ただしこの数字列の中に2が1つでも含まれる場合は、数字列に対する何らの操作も行わない。
"
| - 底辺私立医大を卒業した医者って頭悪いよね？ Part17
767 ：卵の名無しさん[sage]：2020/10/18(日) 15:50:06.26 ID:PYky+T0u - rm(list=ls())
f<- function(t){ # 時刻tでの数列aを返す
a=numeric()
a[1]=0 # a[1]=0, t=0での数列
if(t==0) return(a)
for(i in 2:(t+1)){
x=sample(0:2,1) # 追加する数字0,1,2
a[i]=x
if(x==0){ # 0であれば
i_max=max(which(a[1:(i-1)]==0)) # i_max:追加前の最後の0のindex
b=a[(i_max+1):(i-1)] # b:追加前の最後の0からの追加した0を除く数列
if(all(b!=2)){ # bに2が含まれなければ
b[b==1]=0 # bの１を0に置き換えて
a[(i_max+1):(i-1)]=b # aの配列に置換する
}
}
}
return(a)
}
f(10)
sim <- function(n) sum(f(n)==0)
mean(replicate(1e4,sim(10)))
| - 底辺私立医大を卒業した医者って頭悪いよね？ Part17
768 ：卵の名無しさん[sage]：2020/10/18(日) 16:28:25.78 ID:PYky+T0u - # 書かれた数列から変換後の結果を返す
f<- function(y=c(1,1,0,1,2,0,1)){
t=length(y)
a=numeric(t+1)
a[1]=0 # a[1]=0, t=0での数列
if(t==0) return(a)
for(i in 2:(t+1)){
x=y[i-1]
a[i]=x
if(x==0){ # 0であれば
i_max=max(which(a[1:(i-1)]==0)) # i_max:追加前の最後の0のindex
b=a[(i_max+1):(i-1)] # b:追加前の最後の0からの追加した0を除く数列
if(all(b!=2)){ # bに2が含まれなければ
b[b==1]=0 # bの１を0に置き換えて
a[(i_max+1):(i-1)]=b # aの配列に置換する
}
}
}
return(a)
}
f(c(1,1,0))
f(c(1,1,0,1,2,0,1))
| - 底辺私立医大を卒業した医者って頭悪いよね？ Part17
769 ：卵の名無しさん[sage]：2020/10/18(日) 16:28:45.18 ID:PYky+T0u - > f(c(1,1,0))
[1] 0 0 0 0
> f(c(1,1,0,1,2,0,1))
[1] 0 0 0 0 1 2 0 1
|
|