底辺私立医大を卒業した医者って頭悪いよね？ Part17

底辺私立医大を卒業した医者って頭悪いよね？ Part17at HOSP

底辺私立医大を卒業した医者って頭悪いよね？ Part17 - 暇つぶし2ch■コピペモード
□スレを通常表示
□オプションモード
□このｽﾚｯﾄﾞのURL
■項目テキスト

300:記録して精算して持ち金を０にリセットする。閏日はギャンブルなしとして  これを毎年繰り返したときの黒字の日数の分布は年数を増やすと正規分布になるか？

301:卵の名無しさん
 20/07/25 06:32:51.87 Wl4jphdd.net
flip <- function(N=1000,n=10){ # n+ heads in N coin flips 
library(gmp) 
library(expm) 
'%.%' <- function(x,y) sum(x*y) 
M=cbind(rbind(rep(1,n),diag(1,n,n)),c(rep(0,n),2)) 
v=as.vector.bigq(c(1,rep(0,n))) 
P=as.bigq((M/2)%^%N %*% v)[n+1] 
list(P=P,p=as.numeric(P)) 
} 
flip(1000,10)$p 
flip(1000,10)$p - flip(1000,11)$p 
flip.max <- function(N=1000,n=10){ 
P=flip(N,n)$P - flip(N,n+1)$P 
list(P=P,p=as.numeric(P)) 
} 
sapply(11:13, function(x) flip.max(1e4,x)$p) 
sapply(14:16, function(x) flip.max(1e5,x)$p) 
sapply(17:19, function(x) flip.max(1e6,x)$p) 
sapply(21:23, function(x) flip.max(1e7,x)$p) 
sapply(24:26, function(x) flip.max(1e8,x)$p)

302:卵の名無しさん
 20/07/25 06:42:53.28 Wl4jphdd.net
分数表示版がようやく完成。 
Flip <- function(N=100,n=5){ # n+ heads in N coin flips 
library(gmp) 
library(expm) 
'%.%' <- function(x,y) sum(x*y) 
M=as.bigq(cbind(rbind(rep(1,n),diag(1,n,n)),c(rep(0,n),2))) 
v=as.vector.bigq(c(1,rep(0,n))) 
d=diag(1,n+1,n+1) 
for(i in 1:N) d=d%*%M/2 
P=(d%*%v)[[n+1]] 
list(P=P,p=as.numeric(P)) 
} 
Flip(1000,10) 
Flip.max <- function(N=1000,n=10){ 
P=Flip(N,n)$P - Flip(N,n+1)$P 
list(P=P,p=as.numeric(P)) 
} 
Flip.max(1000,10)

次ページ