臨床統計もおもしろいですよ、その1at HOSP
臨床統計もおもしろいですよ、その1 - 暇つぶし2ch900:卵の名無しさん
18/10/30 00:28:02.10 Eer9XEMu.net
暴走せずに200万個を数えてくれる
>treasure 5 6 8
p1st = 1827737, q1st = 1825076, draw = 2200112

901:卵の名無しさん
18/10/30 01:00:07.00 Eer9XEMu.net
\hs>treasure 5 6 10
p1st = 7995426, q1st = 8023257, draw = 14026332

902:卵の名無しさん
18/10/30 07:08:46.76 Eer9XEMu.net
>treasure 5 6 9
p1st = 4130886, q1st = 4139080, draw = 6037184

903:卵の名無しさん
18/10/30 07:49:18.11 BVJsV4IZ.net
やはり事務員には難しすぎたようだね
さて 算数の問題です
URLリンク(m.youtube.com)
この動画は part1 が 2009年に投稿され
最終回まであと1話となっています
最終回の登校日はいつになるでしょう?
あらゆる手段を用いて求めましょう

904:卵の名無しさん
18/10/30 08:53:02.81 kXLKssbE.net
>>840
5×6マスで宝の数を10まで増やしていくと、
D:\bin>for %i in (1,2,3,4,5,6,7,8,9,10) do treasure 5 6 %i
D:\bin>treasure 5 6 1
p1st = 14, q1st = 14, draw = 2
D:\bin>treasure 5 6 2
p1st = 203, q1st = 197, draw = 35
D:\bin>treasure 5 6 3
p1st = 1801, q1st = 1727, draw = 532
D:\bin>treasure 5 6 4
p1st = 11418, q1st = 11008, draw = 4979
D:\bin>treasure 5 6 5
p1st = 55469, q1st = 54036, draw = 33001
D:\bin>treasure 5 6 6
p1st = 215265, q1st = 211894, draw = 166616
D:\bin>treasure 5 6 7
p1st = 685784, q1st = 680768, draw = 669248
D:\bin>treasure 5 6 8
p1st = 1827737, q1st = 1825076, draw = 2200112
D:\bin>treasure 5 6 9
p1st = 4130886, q1st = 4139080, draw = 6037184
D:\bin>treasure 5 6 10
p1st = 7995426, q1st = 8023257, draw = 14026332
1:同等
1~8:短軸探索有利
9、10:長軸探索有利
という結果になった。
久しぶりにバッチコマンドを書いた

905:卵の名無しさん
18/10/30 08:53:48.77 kXLKssbE.net
>>844
リンクは踏まないから問題文をかけよ。

906:卵の名無しさん
18/10/30 09:01:44.04 Eer9XEMu.net
>>844
>あらゆる手段を用いて求めましょう
答: お前が答を書け
解答終了!!!
>>840
数学板の賢者がcでプログラムするのを期待していたが、cだと最初にメモリを割り当てなくちゃいけないから難しいんじゃないかな?
動的メモリ確保にしてメモリ不足を回避する必要がでてくると思う。
Rで30個から15個取�


907:闖oす組み合わせをしようとしたら > nrow(combn(30,15)) Error: cannot allocate vector of size 8.7 Gb のメッセージがでて実行不能だった。



908:卵の名無しさん
18/10/30 10:25:52.99 Eer9XEMu.net
30個から15個取り出す組み合わせは155117520通りあるので
シミュレーションで頻度を出して確率を推測してみることにした。
m=5
n=6
mn=m*n
P=as.vector(matrix(1:mn,ncol=n,byrow=T))
Q=1:mn
sim <- function(k=15,all=FALSE){
tre=sample(mn,k)
p=q=numeric(k)
for(i in 1:k){
p[i]=which(tre[i]==P)
q[i]=which(tre[i]==Q)
}
p_order=ifelse(all,max(p),min(p))
q_order=ifelse(all,max(q),min(q))
p_order - q_order
}
choose(30,15)
re15=replicate(1e5,sim(15))
mean(re15<0) < mean(re15>0)
c(p1st_sim=mean(re15<0), q1st_sim=mean(re15>0) ,draw_sim=mean(re15==0))
> c(p1st_sim=mean(re15<0), q1st_sim=mean(re15>0) ,draw_sim=mean(re15==0))
p1st_sim q1st_sim draw_sim
0.16354 0.16515 0.67131

909:卵の名無しさん
18/10/30 12:49:00.84 BVJsV4IZ.net
おらっ 新しい算数の問題だ
自分で回答が出せなければ0点だぞW
URLリンク(www.nicovideo.jp)

910:卵の名無しさん
18/10/30 14:03:59.31 Eer9XEMu.net
treasures <- function(m,n,k){
mn=m*n
Q=1:mn
(mat=matrix(Q,ncol=n,byrow=T))
(P=as.vector(mat))
tre=combn(mn,k)
PQ <- function(x){
p=q=p_q=numeric(k)
for(i in 1:k){
p[i]=which(P==x[i])
q[i]=which(Q==x[i])
p_q[i]=p[i]-q[i]
}
return(p_q)
}
re=apply(tre,2,PQ)
result=matrix(rep(NA,k*3),ncol=3)
colnames(result)=c('p1st','q1st','even')
for(i in 1:k){
result[i,] = c(sum(re[i,]<0),sum(re[i,]>0),sum(re[i,]==0))
}
return(result)
}
treasures(4,5,5)
# URLリンク(tpcg.io)

911:卵の名無しさん
18/10/30 14:06:10.19 Eer9XEMu.net
D:\bin>treasure 5 6 11
p1st = 13346984, q1st = 13395944, draw = 27884372
D:\bin>treasure 5 6 12
p1st = 19312228, q1st = 19372871, draw = 47808126

912:卵の名無しさん
18/10/30 14:08:53.77 Eer9XEMu.net
>>850
この結果は面白いな
全体としては互角なんだが、勝負のルールによって勝者がかわるのがこれで分かった。
> treasures(4,5,5)
p1st q1st even
[1,] 1948 9680 3876
[2,] 5488 10016 0
[3,] 7752 7752 0
[4,] 10016 5488 0
[5,] 9680 1948 3876

913:卵の名無しさん
18/10/30 16:01:45.84 Eer9XEMu.net
>>851
メモリ不足でクラッシュするかと思ったが案外、頑強だな。
エラー処理皆無でコードを書いたのだが。
1億を超える演算をこなしている。
Haskellのコンパイラーは優秀で感心。

D:\bin>treasure 5 6 13
p1st = 24301031, q1st = 24358063, draw = 71100756

914:卵の名無しさん
18/10/30 16:05:37.28 BVJsV4IZ.net
おい、婚活に失敗したド底辺
算数の問題の答えはまだぁ
まぁオマエの頭じゃ無理だろうな

915:卵の名無しさん
18/10/30 16:06:30.58 BVJsV4IZ.net
このスレは事務員が日がな一日妄想を垂れ流し
見物人たちがそれを見てフルボッコにするスレである
事務員 とは?
同一者の別名として、薄汚いジジイ、国試浪人、統計野郎、自称医科歯科、がある
自称医科歯科卒、実際は九州の駅弁国立出身で、卒業後は実家の東京に帰り国試浪人となり23年間経過、家庭教師バイトなどを経て現在は事務員、とされている
本人は、勤務医でたまに開業医の手伝いと内視鏡バイト、専門医なし、と主張しているが、連日連夜の異常な書き込み数の多さからは、勤務医とは考え難く、彼の職業がたとえ何であったとしても、人としてド底辺なことは間違いがない
自ら事務員であると告白したこともある
スレリンク(hosp板:108番)-109
自作自演も頻繁に行なっている
事務仕事の際に手にした他人の給与明細や同窓会費振込票を盗撮し、自分のものとしてアップロードしたこともある、犯罪なんじゃねーの?
統計絶対主義でマウントを取ろうとするが、大学卒業後の人生の大部分においてそれが何の役に立つのかはなはだ疑問である
統計に基づいた重病の治療成績を患者に説明しても、患者の希望する治療は患者の事情によって異なることは一定年数以上の臨床経験があるものなら誰でも知っていることだ
それが理解できない彼は、彼自身が素性を明らかにしているようなものである
コンプ薬屋同様に、私立医への執拗な誹謗中傷を数年余にわたって続けている
裏口批判をしているが、実際にそれが可能な資産家の子弟がいるであろう右派スレには近寄らない、一度書き込んだ形跡があったが指摘したら別人だと言って逃げた
弱いものが多そうなスレしか狙わないのもコンプ薬屋と共通する異常者の特徴である
劣等コンプレックスの塊なんだろうな
自分の主張に対して都合の悪い突込みが入ると、相手を私立医と決めつけて、さてはシリツだな、裏口バカ、というセリフでワンパターンに返すことが多い
国試問題を挙げて簡単だとほざいてますが、本番で通らなければお話になりません
いくらネットでわめいていても、医師国試の本番で通らなければお話になりません
大事なことなので2回言いました

916:卵の名無しさん
18/10/30 16:12:20.65 Eer9XEMu.net
>>854
問題文がないのだから、無理だよ。
>あらゆる手段を用いて求めましょう
答: お前が答を書け
解答終了!!!

917:卵の名無しさん
18/10/30 18:54:55.33 Eer9XEMu.net
treasure_order <- function(m,n,k){
mn=m*n
Q=1:mn
(mat=matrix(Q,ncol=n,byrow=T))
(P=as.vector(mat))
tre=combn(mn,k)
PQk <- function(x){
p=q=p_q=numeric(k)
for(i in 1:k){
p[i]=which(P==x[i])
q[i]=which(Q==x[i])
}
pq=sort(unique(c(p,q)))
np=as.numeric(pq[1:k] %in% p)
nq=as.numeric(pq[1:k] %in% q)
10*np + nq
}
(re=apply(tre,2,PQk))
(result=matrix(NA,nrow=k,ncol=3))
for(i in 1:k){
result[i,]=c(sum(re[i,]==10),sum(re[i,]==1 ),sum(re[i,]==11))
}
colnames(result)=c('P_gets','Q_gets','both_get')
result
}
treasure_order(4,5,5)

918:過去ログ ★
[過去ログ]
■ このスレッドは過去ログ倉庫に格納されています


最新レス表示
レスジャンプ
類似スレ一覧
スレッドの検索
話題のニュース
おまかせリスト
オプション
しおりを挟む
スレッドに書込
スレッドの一覧
暇つぶし2ch