臨床統計もおもしろいですよ、その１ [無断転載禁止]©2ch.net

**卵の名無しさん** · 2017/05/03(水) 20:04:54.62

　
　内科認定医受験の最低限の知識、
　製薬会社の示してくる臨床データ、
　論文の考察、
　論文を書くときの正当性、
　というのが、臨床統計の今までの目的の大きい部分でしたが、
　
　AI=機械学習の基本も、結局は統計学と確率に支配されます。
　そういう雑多な話をするスレです。
　

**卵の名無しさん** · 2018/09/17(月) 20:23:03.24

診療所から病院に患者を救急搬送する。
病院から救急車が診療所に向かっており１０時到着予定と連絡が入った。
患者が９時に急変したため診療所の普通車で病院に向かって救急車と出会ったら救急車に患者を移して搬送し病院到着を早めることになった。救急車の方が速く走れる。
９時５０分に救急車に乗り移ることができた。
病院到着は予定より何分早まるか述べよ。
車は定速走行とし、乗り換えに要する時間は考慮しない。

**卵の名無しさん** · 2018/09/17(月) 22:45:36.98

# How it works
dec2n <- function(num, N, digit = 3){ # decimal to 0,1,..,n-1 vector
if(num <= 0 & digit <= 0) return(NULL)
else{
append(dec2n(num%/%N, N ,digit-1), num%%N)
}
}
dec2n(42,5,4) # 0 1 3 2
N=5
num=42 ; digit=4
num%%N # 2 last=2
num%/%N # 8 (digit=3)
num%/%N%%N # 2nd last= 3
num%/%N%/%N # 1 (digit=2)
num%/%N%/%N%%N # 3rd last =1
num%/%N%/%N%/%N # 0 (digit=1)
num%/%N%/%N%/%N%%N # 4th last = 0
num%/%N%/%N%/%N%/%N # 0 (digit=0)

system.time(for(i in 0:10^5) dec2n(i,7))
system.time(for(i in 0:10^5) dec2nw(i,7))

**卵の名無しさん** · 2018/09/17(月) 22:46:00.21

# while loop version for dec2n
dec2nw <- function(num, N, digit = 4){
r=num%%N
q=num%/%N
while(q > 0 | digit > 1){
r=append(q%%N,r)
q=q%/%N
digit=digit-1
}
return(r)
}

**卵の名無しさん** · 2018/09/18(火) 01:40:12.09

https://www.tutorialspoint.com/execute_r_online.php

**卵の名無しさん** · 2018/09/18(火) 06:21:30.53

>>701
http://tpcg.io/QCYcf4

**卵の名無しさん** · 2018/09/18(火) 07:43:40.10

#include<stdlib.h>
#define N 5

int b[N];
//char * u[] = {"0","1","2","3","4"};
char * u[] = {"志学","而立","不惑","知命","耳順","従心"};
void fill(int k){
int i,j;

if(k < N){
for(i=0;i<N;i++){
for(j=0; j < k && b[j]!=i; j++);
;
if(j==k){
b[k] = i;
fill(k+1);
}
}
}
else{
for(i = 0;i < N; i++){
printf("%s ", u[b[i]] );
}
printf("\n");
}
}
int main(){
printf("all around\n");
fill(0);
return 0;
}

**卵の名無しさん** · 2018/09/19(水) 22:38:29.64

draft

is.sorted <- function(x){
if(length(x)==1) return(TRUE){
else{
if(x[1]<=x[2]) is.sort(x[-1]){
else return(FALSE)
}
}
}

**卵の名無しさん** · 2018/09/19(水) 22:44:21.96

is.sorted <- function(x){
if(length(x)==1){ return(TRUE)
}else{
if(x[1]<=x[2]) is.sorted(x[-1])
else return(FALSE)
}
}

**卵の名無しさん** · 2018/09/20(木) 20:45:13.64

ニューロン治療できる医療大麻オイルの紹介
https://plaza.rakuten.co.jp/denkyupikaso/diary/201809180000/

ニューロン=人工知能のモデルとなっている神経細胞のやりとりするところ
この細胞の伝達する場所に大麻受容体がありますこの受容体を通って治療効果を得る大麻を医療大麻と呼びます

人間に大麻受容体があったなんて不思議ですよね

**卵の名無しさん** · 2018/09/22(土) 08:25:11.10

>>711 は国試に２３回も落ちてて
くるくるぱーの裏口バカに
なっちゃってるのらぁあぁぁ
Ｆラン事務員の濃ゆぅぅい生ガキ汁
ド底辺の臭いが落ちないよぉ
んほおぉぉぉおぉぉ

**卵の名無しさん** · 2018/09/22(土) 08:33:45.20

>>713
>711って再帰呼び出ししているけど何をするスクリプトかわかる？

**卵の名無しさん** · 2018/09/22(土) 13:39:33.84

#include<stdio.h>
#include<stdlib.h>
#define N 8
#define W 2*N-1
int b[N],row[N],up[W],down[W],count=0;
void fill(int k){
int i,j;
if(k < N){
for(i=0;i<N;i++){
if(!row[i] && !up[i+k] && !down[i-k+N-1]){
b[k] = i;
row[i]=up[i+k]=down[i-k+N-1]=1;
fill(k+1);
row[i]=up[i+k]=down[i-k+N-1]=0;
}
}
}
else{
printf("# %2d\n",++count);
for(i =0; i < N; i++){
for(j = 0; j < N; j++){
if(b[i] == j) printf(" Ｑ ");
else printf(" ・ ");
}
printf("\n");
}
}
}
int main(){
fill(0);
return 0;　}

**卵の名無しさん** · 2018/09/22(土) 21:34:00.39

#include <stdio.h>

int arr[] ={1,2,3,4,5};
int total = (sizeof(arr) / sizeof(arr[0]));

void heapify(int arr[], int i)
{
int lft = i * 2;
int rgt = lft + 1;
int grt = i;
if (lft <= total && arr[lft] > arr[grt]) grt = lft;
if (rgt <= total && arr[rgt] > arr[grt]) grt = rgt;
if (grt != i) {
int tmp=arr[grt];
arr[grt]=arr[i];
arr[i]=tmp;
heapify(arr, grt);
}
}

int main(){
heapify(arr,0);
for (int i = 0; i < total; i++) printf("%d ", arr[i]);
printf("\n");
}

**卵の名無しさん** · 2018/09/23(日) 18:58:09.75

library(gtools)
k=3
perm=unique(permutations(3*k,3*k,rep(1:3,k),set=F,rep=F))
is3 <- function(x){
n=length(x)
y=c(x,x[1])
re=NULL
for(i in 1:(n-1)) re[i]=all((1:3) %in% y[i:(i+2)])
any(re)
}
i=(!apply(perm, 1, is3))
perm1=perm[i,]

r1 <- function(x){ # rotate by one bead
n=length(x)
y=numeric(n)
y[1]=x[n]
for(i in 1:(n-1)){
y[i+1]=x[i]
}
return(y)
}
rn <- function(x){ # every rotation
n=length(x)
mat=matrix(rep(NA,n^2),ncol=n)
mat[,1]=x
for(i in 1:(n-1)){
mat[,i+1]=r1(mat[,i])
}
return(t(mat))
}

**卵の名無しさん** · 2018/09/23(日) 18:58:40.06

same <- function(x,base){
f=function(a,b=base){ # is equal to base
all(a==b)
}
mat=rbind(rn(x),rn(rev(x))) # rotation and/or symmetric conversion
any(apply(mat,1,f))
}

mat=perm1
tail(perm1)
exvar <- function(row,mat=perm1){
var=rn(mat[row,])
n.var=nrow(var)
idx=NULL
for(i in 1:n.var){
idx=c(idx,which(apply(mat,1,function(x) same(x,var[i,]))))
}
min(idx)
}

index=sapply(1:nrow(perm1),exvar) # several minutes for k=3
re=perm1[unique(index),]
re
print(matrix(c("赤","白","青")[re],ncol=3*k),quote=F)

**卵の名無しさん** · 2018/09/24(月) 16:11:21.71

same <- function(x,base){
f=function(a,b=base){ # is equal to base
all(a==b)
}
mat=rbind(rn(x),rn(rev(x))) # rotation and/or symmetric conversion
any(apply(mat,1,f))
}

mat=perm1
core=NULL　# 基本解を集める行列
while(nrow(mat)){
var=rn(mat[1,])　# matの先頭行のバリアント
n.var=nrow(var)
idx=NULL # 先頭行とそのバリアントのindex
for(i in 1:n.var){ # 先頭行およびバリアント行のindexを返す
idx=c(idx,which(apply(mat,1,function(x) same(x,var[i,]))))
}
i=sort(unique(idx))
core=rbind(core,mat[i[1],])
mat=mat[-i,] # 一度、検討した行と同一相当行を除いたmatで再評価に回す
}
core

print(matrix(c("赤","白","青")[core],ncol=3*k),quote=F)

**卵の名無しさん** · 2018/09/24(月) 20:26:43.18

library(compiler)
cuniq=cmpfun(unique)
jperm <- function(n,r,v) {
if (r == 1)
matrix(v, n, 1)
else if (n == 1)
matrix(v, 1, r)
else {
X <- NULL
for (i in 1:n) X <- rbind(X, cbind(v[i], jperm(n - 1, r - 1, v[-i])))
cuniq(X)
}
}
cjperm=cmpfun(jperm)
v=rep(1:3,4)
n=r=length(v)
perm4=cjperm(n,r,v)

**卵の名無しさん** · 2018/09/24(月) 20:44:38.00

#include<stdio.h>
#define N 6

int b[N];
char * u[] = {"P","E","P","P","E","R"};

void p_table(int k){
int i;

if(k < N){
for(i=0;i<N;i++){
b[k] = i;
p_table(k+1);
}
}
else{
for(i = 0;i < N; i++){
printf("%s ", u[b[i]] );
}
printf("\n");
}
}
int main(){
p_table(0);
return 0;
}

**卵の名無しさん** · 2018/09/24(月) 21:05:06.98

// https://www.geeksforgeeks.org/print-all-permutations-of-a-string-with-duplicates-allowed-in-input-string/

// Program to print all permutations of a string in sorted order.
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <stdbool.h>

/* Following function is needed for library function qsort(). */
int compare(const void *a, const void * b)
{
return ( *(char *)a - *(char *)b );
}

// A utility function two swap two characters a and b
void swap(char* a, char* b)
{
char t = *a;
*a = *b;
*b = t;
}

**卵の名無しさん** · 2018/09/24(月) 21:05:27.97

// This function finds the index of the smallest character
// which is greater than 'first' and is present in str[l..h]
int findCeil(char str[], char first, int l, int h)
{
// initialize index of ceiling element
int ceilIndex = l;

// Now iterate through rest of the elements and find
// the smallest character greater than 'first'
for (int i = l+1; i <= h; i++)
if (str[i] > first && str[i] < str[ceilIndex])
ceilIndex = i;

return ceilIndex;
}

**卵の名無しさん** · 2018/09/24(月) 21:05:44.12

// Print all permutations of str in sorted order

void sortedPermutations(char str[])
{
// Get size of string
int size = strlen(str);

// Sort the string in increasing order
qsort(str, size, sizeof( str[0] ), compare);

// Print permutations one by one
bool isFinished = false;
while (!isFinished)
{
// print this permutation
static int x = 1;
printf("%d %s \n", x++, str);

// Find the rightmost character which is smaller than its next
// character. Let us call it 'first char'
int i;
for (i = size - 2; i >= 0; --i)
if (str[i] < str[i+1])
break;

**卵の名無しさん** · 2018/09/24(月) 21:06:30.92

// If there is no such chracter, all are sorted in decreasing order,
// means we just printed the last permutation and we are done.
if (i == -1)
isFinished = true;
else
{
// Find the ceil of 'first char' in right of first character.
// Ceil of a character is the smallest character greater than it
int ceilIndex = findCeil(str, str[i], i + 1, size - 1);

// Swap first and second characters
swap(&str[i], &str[ceilIndex]);

// Sort the string on right of 'first char'
qsort(str + i + 1, size - i - 1, sizeof(str[0]), compare);
}
}
}

// Driver program to test above function
int main()
{
char str[] = "";
sortedPermutations( str );
return 0;
}

**卵の名無しさん** · 2018/09/26(水) 11:07:08.91

b=0
w=1
bwb=function(x){
n=length(x)
if(n<3) return(F)
else x[n-2]==b & x[n-1]==w & x[n]==b
}
f = function(){
x=NULL
while(!bwb(x)){
x=c(x,sample(c(b,w),1,prob=c(0.2,0.8)))
}
c(W=sum(x==w),B=sum(x==b))
}

re=replicate(1e4,f())
mean(re[1,])
mean(re[2,])

**卵の名無しさん** · 2018/09/26(水) 22:08:29.29

library(Rmpfr)
f <- function(NN,prec=1000){
for(k in 1:NN){
for(m in k:NN){
for(n in m:NN){
K=mpfr(factorial(k),prec)
M=mpfr(factorial(m),prec)
N=mpfr(factorial(n),prec)
if(K+M==N) cat(k,'! + ',m,'! = ',n,'!\n')
}
}
}
}
try(f(100))

**卵の名無しさん** · 2018/10/02(火) 21:03:43.40

print $ foldr (-) 0 [1,2,3]
print $ foldr (-) 0 ( 1:(2:(3:[])) )
print $ 1-(2-(3-0))

**卵の名無しさん** · 2018/10/02(火) 22:17:49.03

elem_r y ys = foldr (\x acc -> if x==y then True else acc) False ys

elem_r 39 [3.6..]

elem_l y ys = foldl(\acc x -> if x==y then True else acc) False ys

{- ERROR!!
elem_l 39 [3,6..]

**卵の名無しさん** · 2018/10/03(水) 13:45:16.24

map' f xs = foldr (\x acc -> f x:acc) [] xs

elem_r y ys = foldr (\x acc -> if x==y then True else acc) False ys
elem_l y ys = foldl(\acc x -> if x==y then True else acc) False ys
main = do
print $ zipWith (\x y -> (x*30+3)/y) [5,4,3,2,1] [1..5]
print $ map (\(a, b) -> a+b) [(1,2),(3,5),(6,3),(2,6),(2,5)]
print $ zipWith (flip (++)) ["love you","is an angel"] ["I ","Mei "]
print $ map (flip subtract 20) [1..5]

print $ foldl (\a b -> a+b) 0 [3,5,2,1]
print $ foldl (-) 0 [1,2,3]
print $ ((0-1)-2)-3

print $ foldr max 0 [1,2,3]
print $ 1 `max` (2 `max`(3 `max` 0))
print $ foldr (-) 0 [1,2,3]
print $ foldr (-) 0 ( 1:(2:(3:[])) )
print $ 1-(2-(3-0))
print $ foldr (\a b -> a-b) 0 ( 1:(2:(3:[])) )

**卵の名無しさん** · 2018/10/08(月) 17:55:51.64

import Data.List
divisor n = find (\m -> n `mod` m ==0 )[2..floor.sqrt.fromIntegral $ n]
main = do
print $ divisor $ (10^19-1) `div` 9

**卵の名無しさん** · 2018/10/08(月) 18:06:39.22

import Data.List
divisor n = find (\m -> n `mod` m ==0 )[2..floor.sqrt.fromIntegral $ n]
divisors n = filter (\m -> n `mod` m ==0 )[2..n-1]
main = do
print $ divisor $ (10^19-1) `div` 9

**卵の名無しさん** · 2018/10/08(月) 20:25:22.38

１００万個の○が円形に並んでいます。
図のように、まず１つの○に色をぬり、
次にその●から時計回りに１０８個進んで止まり、そこにある○をぬります。
さらにその●から
時計回りに１０８個進んで止まり、
そこにある○をぬり、以下同じ作業をくり返していきます。
すでに色がぬられた●に止まったときに終了とするとき、
何個の○をぬることができますか？

foo jump spot = length $ takeWhile (\x -> x*jump `mod` spot /=0) [1..]

1 + foo 108 1000000

**卵の名無しさん** · 2018/10/08(月) 20:30:18.22

ｍ個の○が円形に並んでいます。
図のように、まず１つの○に色をぬり、
次にその●から時計回りにｎ個進んで止まり、そこにある○をぬります。
さらにその●から
時計回りにｎ個進んで止まり、
そこにある○をぬり、以下同じ作業をくり返していきます。
すでに色がぬられた●に止まったときに終了とするとき、
何個の○をぬることができますか？

ｍ　/ gcd(m.n)

m,nが互いに素なら全部塗れる。

**卵の名無しさん** · 2018/10/08(月) 20:48:49.62

gcd a b
| a `mod` b == 0 = b
| otherwise = gcd b $ a `mod` b

**卵の名無しさん** · 2018/10/10(水) 08:15:31.85

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

**卵の名無しさん** · 2018/10/12(金) 06:12:09.09

知り合いから教えてもらった自宅で稼げる方法
興味がある人はどうぞ
みんながんばろうねぇ『羽山のサユレイザ』で

K5K

**卵の名無しさん** · 2018/10/12(金) 20:06:27.81

ある医大で合格率の男女比が1.2で男子有意という結果だったという。
定員100で男子800人女子200人が受験して合格率の男女比が
1.2であったときに統計的には有意差があると言えるか？

**卵の名無しさん** · 2018/10/13(土) 08:03:38.90

man=900
woman=100
pass=100
total=man+woman
p=pass/total
i=0:100
#
sum((choose(woman,i)*choose(man,pass-i)/choose(total,pass))*(i/woman > (pass-i)/man))
sum((choose(woman,i)*choose(man,pass-i)/choose(total,pass))*(i/woman == (pass-i)/man))
sum((choose(man,i)*choose(woman,pass-i)/choose(total,pass))*(i/man > (pass-i)/woman))

sum((choose(man,i)*choose(woman,pass-i)/choose(total,pass))*(i/man > 1.2*(pass-i)/woman))

**卵の名無しさん** · 2018/10/13(土) 08:04:10.37

> sum((choose(woman,i)*choose(man,pass-i)/choose(total,pass))*(i/woman > (pass-i)/man))
[1] 0.4160339
> sum((choose(woman,i)*choose(man,pass-i)/choose(total,pass))*(i/woman == (pass-i)/man))
[1] 0.1389853
> sum((choose(man,i)*choose(woman,pass-i)/choose(total,pass))*(i/man > (pass-i)/woman))
[1] 0.4449808
> sum((choose(man,i)*choose(woman,pass-i)/choose(total,pass))*(i/man > 1.2*(pass-i)/woman))
[1] 0.3090334

**卵の名無しさん** · 2018/10/13(土) 08:28:07.92

ある医大で合格率の男女比が1.2で男子有意という結果だったという。
定員100で男子800人女子200人が受験して合格率の男女比が
1.2であったときに統計的には有意差があると言えるか？

f <- function(ratio=1.2,man=800,woman=200,pass=100){
total=man+woman
i=0:pass
sum((choose(man,i)*choose(woman,pass-i)/choose(total,pass))*(i/man > ratio*(pass-i)/woman))
}
f()

**卵の名無しさん** · 2018/10/13(土) 10:57:49.86

draft
{-ある医大で合格率の男女比が1.2で男子有意という結果だったという。
定員100で男子800人女子200人が受験して合格率の男女比が
1.2であったときに統計的には有意差があると言えるか？
-}
choose n r = [1..n] `div` [1..n-r] `div` [1..r]
ratio=1.2
pass=100
man=800
woman=200
total=man+woman
ii=[0..pass]
male_fem i =
| choose(man,i)*choose(woman,pass-i)/choose(total,pass) <ratio 0
| otherwise choose(man,i)*choose(woman,pass-i)/choose(total,pass
sum $ map male_fem ii

**卵の名無しさん** · 2018/10/13(土) 12:59:15.69

>>736
ド底辺シリツ医大卒って医療従事者にマウントするのは不可能だから事務員や国試浪人認定するしかないんだなぁ。

こういうのにサクッと答えてド底辺シリツ医大卒でも高卒の基礎学力くらいあるのを示せばいいのに。

ある医大で合格率の男女比が1.2で男子優位という結果だったという。
定員100で男子800人女子200人が受験して合格率の男女比が
1.2であったときに統計的には有意差があると言えるか？

**卵の名無しさん** · 2018/10/13(土) 15:21:07.70

>>743
Rで数えると

f <- function(ratio=1.2,man=800,woman=200,pass=100){
total=man+woman
i=0:pass
sum((choose(man,i)*choose(woman,pass-i)/choose(total,pass))*(i/man > ratio*(pass-i)/woman))
}
f()

Haskellで数えると

ratio=1.2
pass=100
man=800
woman=200
total=man+woman
ii=[0..pass]
male_fem i= if r < ratio then 0 else p
where a = (choose man i)*(choose woman (pass-i))
b = choose total pass
m=fromIntegral i/ fromIntegral man
w=fromIntegral (pass-i)/fromIntegral woman
r = m/w
p = (fromIntegral a) / (fromIntegral b)

main = do
print $ sum ( map male_fem ii)

**卵の名無しさん** · 2018/10/14(日) 06:36:53.74

f <- function(ratio=1.2,wm=0.2,total=1000,pass=100){
woman=total*wm
man=total-woman
i=0:pass
sum((choose(man,i)*choose(woman,pass-i)/choose(total,pass))*(i/man > ratio*(pass-i)/woman))
}
f()
r=seq(1,3,length=1000)
p=sapply(r,f)
plot(r,p,pch=19,type='n',bty='l',xlab='M/F ratio',ylab='probability')
abline(h=0.01,lty=3)
g <- function(wm) lines(r,sapply(r,function(r)f(r,wm)),col=10*wm,lwd=10*wm)
for(i in 1:5/10) g(i)
legend('center',bty='n',lwd=1:5,col=1:5,
legend=c('10%','20%','30%','40%','50%'))

pass=100
wm=0.2
curve( (x/(1-wm))/((pass-x)/wm),0,pass)
uniroot(function(x,u0=1.2) x*wm/(1-wm)/(pass-x)-u0, c(0,pass))
(83/800)/(17/200)
prop.test(c(83,17),c(800,200))

**卵の名無しさん** · 2018/10/14(日) 07:14:30.48

pass=100
wm=0.2
curve( (x/(1-wm))/((pass-x)/wm),0,pass)
uniroot(function(x,u0=1.2) x*wm/(1-wm)/(pass-x)-u0, c(0,pass))
(83/800)/(17/200)
prop.test(c(83,17),c(800,200))

m_suc = function(rate) (pass* rate * (wm - 1))/(rate* (wm - 1) - wm)
m_suc(1.2)
curve(m_suc(x),1,5)

total=1000
wm=0.5
h <- function(mf){
M=ceiling(m_suc(mf))
F=pass-M
prop.test(c(M,F),c(total*(1-wm),total*wm))$p.value
}
r=seq(1,3,length=1000)
plot(r,sapply(r,h),type='l')
abline(h=0.05,lty=3)
uniroot(function(x)h(x)-0.05,c(1,2))$root
m_suc(1.43)
59/41
prop.test(c(59,41),c(500,500),correct=F)
prop.test(c(60,40),c(500,500),correct=F)

**卵の名無しさん** · 2018/10/17(水) 18:14:18.05

　　　　　　　　,イ　　　　　　　　　　　　　　　　　　　　ヽ
　　　　　 /　　　　　　　　　　ヽ　ヽ　　　　　　　　　ヽ
　　　　　　 /　/　　　　|　　　{.　ハハ　　　　　　　　　　　ヽ
　　　　 /　 :!　　　　│ j　八　　 |　 |　　　　　　　　　　　　)　　　　　
　　　　　,′ i |　｛. 　 i |　{＼ﾊ　| i|　ﾘ　　　 )､　　　　　　　/　　　国試浪人事務員は薄汚いジジイ
　　　　ｉ　i |　八　ｉ Vx炒¨＾jﾉ}八/　　　　( ヽ　　　　　/　　　　敗者は敗者らしく
　　　　 | ｉ　从 { { ＼八{　////　ｲ　　　　　)ヾﾉ　　　　　(　　　　　この世界のなるべく闇いところで
　　　 Viﾊ＼|x炒＾j　　　　　　　(　　　　　(ｲ　　　　 /　　　　　　謝りながら生きていくといいのよ
　　　　＼＼_V//く　_　　　　 )　　　　　)　ヾ　　　ゝ,
　　　　　　　 `ｰ小.　　　　 _-ｧ　　(　　　　　(　　　　　　ヾ
　　　　　　 |i∧ 　く:::_ノ　　 )　　　　　i　　　　　丶､ﾞゝ
　　　　　　 |ｉ| {丶、　　´　　　(　　　　　| 　　　　iヾﾞ
　　　　　　　ﾉ八{ { {＞　.　　／)　　　　　ｉ　　　　ヾ
　　　　　　　　〃{｛　`ｰ=ミ　｀T爪八　i　 ,　ヾ　　　　八　
　　　　　　　　ヾ人＼　　｀ヽ j{__j ｼ｀ヾﾐﾐﾙ'ﾊヽ.ト＜7＼＿＿＿___
　　　　　　　　　／⌒￢=-ミ#ヾ/　　　　　　　　＼　＼/####/　　｀ヽ
　　　　　　　　/　　　　　　　ヽ##｀ヽ､ノ⌒ヽ　　　　V⌒ゝ､#;;|　　　 '.
　　　　　　　　{　　　　　　　　　ﾞi####ヽ　　八　　　 |'⌒´##ヾ　　　　　|
　　　　　　　　|　　　　　　　　　 |#####＼/:i ＼'ー‐'i"######ﾞ)　　　　 |
　　　　　　　 │　　　　　　　　　|########| :|　( ⌒ )#######|　　　　|

**卵の名無しさん** · 2018/10/21(日) 21:20:01.39

# 縦3マス、横4マスの12マスのうちランダムに選ばれた2マスにそれぞれ宝が眠っている。
# AEIBFJ…の順で縦に宝を探していく方法をとるP君と、
# ABCDEFGH…の順で横に宝を探していく方法をとるQ君が、同時に地点Aから探索を開始した。
# どっちの方が有利？
#
# ABCD
# EFGH
# IJKL

y=LETTERS[1:12]
z=matrix(y,ncol=4,byrow=T)
as.vector(z) # P
as.vector(t(z)) # Q(=y)

x=c(1,1,rep(0,10))
is.P1st <- function(){
Q=sample(x)
z=matrix(Q,ncol=4,byrow=T)
P=as.vector(z)
which.max(P) < which.max(Q)
}
k=1e3
re=replicate(k,mean(replicate(k,is.P1st())))
summary(re)

**卵の名無しさん** · 2018/10/21(日) 21:45:36.09

> x=c(1,1,rep(0,10))
> PQ <- function(){
+ Q=sample(x)
+ z=matrix(Q,ncol=4,byrow=T)
+ P=as.vector(z)
+ c( even=which.max(P) == which.max(Q),
+ p1st=which.max(P) < which.max(Q),
+ q1st=which.max(P) > which.max(Q))
+
+ }
> k=1e6
> re=replicate(k,PQ())
> mean(re['even',]) ; 13/(26+27+13)
[1] 0.197025
[1] 0.1969697
> mean(re['p1st',]) ; 26/(26+27+13)
[1] 0.393803
[1] 0.3939394
> mean(re['q1st',]) ; 27/(26+27+13)
[1] 0.409172
[1] 0.4090909

**卵の名無しさん** · 2018/10/22(月) 02:17:14.19

>>748
一般解

treasure <- function(m=3,n=4,k=2){
y=1:(m*n)
(z=matrix(y,ncol=n,byrow=T))
(P=as.vector(z))
(Q=as.vector(t(z)))
PQ <- function(x){
p=q=numeric(k)
for(i in 1:k){
p[i]=which(P==x[i])
q[i]=which(Q==x[i])
}
min(p)-min(q)
}
tre=combn(m*n,k)
re=apply(tre,2,PQ)
c(P1st=sum(re<0),Q1st=sum(re>0),even=sum(re==0))
}
treasure(3,4,2)
treasure(5,10,3)

**卵の名無しさん** · 2018/10/22(月) 22:25:36.36

|￣￣￣￣￣￣|　..--――--　 __
|　　　国　　 | .:.:.:.:.:.:.:.:..:.:.:.:.:.:.:.:.:.｀.
|　　　試　　 |.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:＼
| 　　　浪　　 |:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:ヽ
| 　　　人　　|.:.:./:.:.:.:.:.:.:.:.:.:.:ヽ:.:.:.:.:.:.:..:..:.:ﾊ
|　　！？　　 |.:.:.|:.:.:.:.:.:.::.:.:.:.:.:.|:.:.:.:.:.:.:.:.:.:.:.:.}
|＿＿＿＿＿＿|:.:.:l.:.:l:.:.:.:.:l.:.:.:.:.:ｌ:.:.:.:.:.:ｌ:.:.:.:.:.:j
　　　　l:.:.:.:.i:.:.:.:ヽl∧ﾊ＼:.:.l:.:/∨ハ:.:/l:.:.:.:ｌ:.:.ｌ
　　　 |:.:.:.:.l:.:.:.:l斗―ミ､ヽ!/ ,. -―ﾐ､ !:.:./:.:.:ｌ
　　　　iﾄ､:.:.|:.:.:.|（・）` 　　 '´（・）` l:.ノ:.:.:.:.:l
　　　　　ﾉ:.:ﾊ:.:.:l ` 一'　　, 　ｰ一'´ﾉ:.:.:.:.|.:.:人
　　　／:.:.:.У:人ι　　　　　　　　　ﾑｲ:.:.:l:.:.:.:./ヽ
　　　(/i:.:{::.ｲ:.:.:ﾊ.丶　　 -―-　　／:.:l:.:.:.:l.:.:.:八
.　　　　ﾚ'´}:.:.:./:.:.}ﾘ}＞　＿_　 .イ人:.:ﾊ:.:.:|:.:ﾘ}
　　　　　 __)ノﾄソ／　＼_　 _／　＼ﾘ　乂へ
　　　イ　　　/##｛　　＿/∨＼＿　 }### i　　＼
　　/　　　 /###レイ.　ヽ人／　ヽ/###.ｌ　　　　ヽ
.　ｌ　　　　/####l　　　 /　ヽ　　 /####ｉ　　　　 }

**卵の名無しさん** · 2018/10/22(月) 22:27:47.24

国試浪人の事務員については専用スレがあるのでそちらへ

https://egg.5ch.net/test/read.cgi/hosp/1537867846/29

**卵の名無しさん** · 2018/10/22(月) 22:28:47.06

妬みマンチクリン　がどういう意味なんだろうと思ってググってみたら
こんなスレが出てきました

自称医科歯科卒が専門医をねたむスレ [転載禁止]©2ch.net
https://egg.5ch.net/test/read.cgi/hosp/1419213994/l50

このスレを辿ってみたら白黒コピーの医師免許の一部が
大部分を隠してアップされてましたが
https://egg.5ch.net/test/read.cgi/hosp/1419213994/437

これって病院で事務用に保管しているコピーですよね

**卵の名無しさん** · 2018/10/23(火) 03:47:27.97

あれ？
朝鮮殺戮殺人学会が警察にテロ工作員を送り込み、

なんと
警察から犯罪ライセンスを与えられ、

監禁罪、薬物大量投与テロ、
などやり、
テロ工作拠点だとバレた

福山友愛病院については？
福山友愛病院　事件　でググれば？

http://sp.nicovideo.jp/watch/sm7997483?ss_id=09863826-7a23-4d06-8791-91e89aac58b7&;ss_pos=19

**卵の名無しさん** · 2018/10/23(火) 07:05:26.69

>>751
スクリプトの最適化とか示せばいいのに
どれでもいいから掲載してあるRのコードをCで高速化してみ。

**卵の名無しさん** · 2018/10/23(火) 07:11:17.37

>>751

専門医スレで誰も正解出せなかった、
これを解くコードを書いて仮想国試浪人程度の頭脳があることを示してくれ。

door to balloon timeが短縮するのに必要な計算。

診療所から病院に患者を救急搬送する。
病院から医師搭乗の救急車が診療所に向かっており１０時到着予定と連絡が入った。
患者の病態が悪化したら、診療所の普通車で病院に向かい救急車と出会ったら
救急車に患者を移して搬送し病院到着を急ぐという計画を立てた。
普通車から救急車への患者の乗り換えで１０分余分に時間がかかる。
道路事情から救急車は病院から診療所への道は
平均時速６０ｋｍで、逆方向は平均時速４５ｋｍで定速走行する。診療所の普通車は信号待ちもあり平均時速３０kmで定速走行する。

何時以降の病態悪化は診療所の車を使わずに救急車の到着を待つ方が病院に早く着くか？

**卵の名無しさん** · 2018/10/23(火) 07:42:35.45

国立卒ならこれくらい答えてシリツとの差をみせてほしいね。

インフルエンザの迅速キットは特異度は高いが感度は検査時期によって左右される。
ある診断キットが開発されたとする。
このキットは特異度は９９％と良好であったが、
感度については確かな情報がない。
事前確率分布として一様分布を仮定する。
５０人を無作為抽出してこの診断キットで診断したところ４０人が陽性であった。
この母集団の有病率の期待値と９５％ＣＩはいくらか？
またこの診断キットの感度の期待値と９５％ＣＩはいくらか

**卵の名無しさん** · 2018/10/23(火) 11:55:53.89

digi = function(x){ # 1000 -> 4 , 999 -> 3
n=ceiling(log10(x))
ifelse(10^n==x,n+1,n)
}

n2a <- function(num){ # nmu to array 122 -> c(1,2,2)
N=10
r=num%%N
q=num%/%N
while(q>0){
r=append(q%%N,r)
q=q%/%N
}
return(r)
}

one2n <- function(x){ # 121 -> 13
a=n2a(x)
k=digi(x)
p=2^((k-1):0)
sum(a*p)
}

x=1212121212
one2n(x)

**卵の名無しさん** · 2018/10/23(火) 13:18:37.27

このスレは事務員が日がな一日妄想を垂れ流し
見物人たちがそれを見てフルボッコにするスレである

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

**卵の名無しさん** · 2018/10/23(火) 14:23:32.78

-- 不定長整数が扱えて遅延IOなのでメモリ不足になりにくいな --

main = do
putStr $ "Type number: "
str <- getLine
let xs = concat $ iterate (\x->[1:n| n<-x] ++ [2:n|n<-x]) [[1],[2]]
n = read str
putStrLn $ concat $ map show (xs !! (n-1))

**卵の名無しさん** · 2018/10/23(火) 14:54:10.02

>>760
コンパイルして実行したけど、１億めはメモリ不足に陥った。

**卵の名無しさん** · 2018/10/23(火) 15:44:29.51

{1, 0, 1, 1, 1, 1, 1, 0, 1, 0, 1, 1, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0}

https://www.wolframalpha.com/input/?i=IntegerDigits+%5B10%5E8,+2%5D

**卵の名無しさん** · 2018/10/23(火) 16:49:32.31

>>761

> digit12(10^8) # 1億め
12222212122221111211111112
> digit12(10^12) # 1兆め
221211122121211212112121112111111111112

Rのコードはここ

http://tpcg.io/D2sseW

**卵の名無しさん** · 2018/10/23(火) 18:48:55.68

>>763 は国試に２３回も落ちた挙句
婚活にも失敗してる素人童貞で
くるくるぱーの裏口バカに
なっちゃってるのらぁあぁぁ
Ｆラン事務員の濃ゆぅぅい生ガキ汁
ド底辺の臭いが落ちないよぉ
んほおぉぉぉおぉぉ

**卵の名無しさん** · 2018/10/23(火) 21:27:28.15

>>763
１０の６８乗を無量大数というらしい

無量大数＋１を２進数表示できるかやってみた。

Prelude> :main
Input integer : 100000000000000000000000000000000000000000000000000000000000000000001
1110110101100011101000100011000111010100110001001111101100100111010011001010011110101010101010000110001111101110010010111101110101001000010101101100010111000100000000000000000000000000000000000000000000000000000000000000000001

さすが不定長整数を扱えるHaskell。

**コンブ薬屋** · 2018/10/23(火) 23:11:54.77

>>765
んで、あんた、どこ卒？
さては国試浪人の事務員だな

**卵の名無しさん** · 2018/10/24(水) 07:56:03.66

>>766
俺、医科歯科卒。二期校時代に入学。

国立卒なのでこういう問題にも答が出せる。

インフルエンザの迅速キットは特異度は高いが感度は検査時期によって左右される。
ある診断キットが開発されたとする。
このキットは特異度は９９％と良好であったが、
感度については確かな情報がない。
事前確率分布として一様分布を仮定する。
５０人を無作為抽出してこの診断キットで診断したところ４０人が陽性であった。
この母集団の有病率の期待値と９５％信用区間はいくらか？
またこの診断キットの感度の期待値と９５％信用区間はいくらか

**卵の名無しさん** · 2018/10/24(水) 13:25:08.39

>>767は国試に２３回も落ちた挙句
婚活にも失敗してる素人童貞で
くるくるぱーの裏口バカに
なっちゃってるのらぁあぁぁ
Ｆラン事務員の濃ゆぅぅい生ガキ汁
ド底辺の臭いが落ちないよぉ
んほおぉぉぉおぉぉ

**卵の名無しさん** · 2018/10/24(水) 14:57:07.89

無職ではまともな女は相手しないだろうw

**卵の名無しさん** · 2018/10/24(水) 18:00:20.55

学習動画をあげといてやろう

http://sp.nicovideo.jp/watch/sm34064890?cp_in=watch_watchRelatedContents

**卵の名無しさん** · 2018/10/24(水) 18:00:40.97

>>768
>国立卒なのでこういう問題にも答が出せる。

を受けて正解を投稿すれば、

ド底辺シリツ医大卒でも国立卒レベルの数学ができることを示す機会なのに

馬鹿を晒す機会にするとはド底辺シリツ医大卒の裏口ガイジらしいぜｗ

**卵の名無しさん** · 2018/10/24(水) 18:11:39.40

>>769

無職のナマポでさえ避けたがるのがド底辺シリツ医大卒

ナマポの症例報告です。

実例↓

517 卵の名無しさん 2018/02/25(日) 11:36:00.56 ID:gq76tAvs
福岡のあの歯科大かな？
歯科口腔外科で抜歯依頼したら爺ちゃんが〇〇歯科大卒の先生は避けて下さいねがあった
けど。
この爺ちゃん、聖マリ卒の先生もよけて下さいと初診時に言ってた札付き爺さん。
生保受給者のくせにね。

**卵の名無しさん** · 2018/10/25(木) 07:04:57.74

裁判所事務官採用試験（大阪高裁管轄）の結果
http://www.courts.go.jp/saiyo/vcms_lf/3all-h30s4-di.pdf
男880　1次通過者224　合格者43
女736 　1次通過者159　合格者87

調整の出来ない1次試験の通過率　　　　男25.5%　女21.6％
調整し放題の面接と称する試験合格率　　男19.2%　女54.7%

f <- function(ratio=1.2,man=800,woman=200,pass=100){
total=man+woman
i=0:pass
sum((choose(man,i)*choose(woman,pass-i)/choose(total,pass))*(i/man > ratio*(pass-i)/woman))
}

f((87/736)/(43/880),736,880,43+87)
f((87/159)/(43/224),159,224,43+87)

prop.test(c(43,87),c(880,736))
prop.test(c(224,159),c(880,736))
prop.test(c(43,87),c(224,159))

**卵の名無しさん** · 2018/10/25(木) 07:38:38.53

> binom.test(43,43+87,224/(224+159),alt="less")

**卵の名無しさん** · 2018/10/25(木) 14:30:04.77

さあシリツのじかんだ

http://cdn-ak.f.st-hatena.com/images/fotolife/h/hanbunningen/20110905/20110905095207.jpg

**卵の名無しさん** · 2018/10/25(木) 18:32:09.06

この方がすっきり。

男880　1次通過者224　合格者43
女736 　1次通過者159　合格者87

x=c(880,736)
y=c(224,159)
z=c(43,87)

fisher.test(cbind(x-y,y))
fisher.test(cbind(y-z,z))
fisher.test(cbind(x-z,z))

prop.test(y,x)
prop.test(z,y)
prop.test(z,x)

**卵の名無しさん** · 2018/10/25(木) 18:35:26.32

>>772
>初診時に言ってた
シリツ医大卒＝裏口バカと認識しているわけだな。

>>
東京医大、本来合格者入学許可へ　今年の受験生５０人

2018年10月25日 02時06分
　東京医科大＝８月、東京都新宿区

　東京医科大が今年の入試で本来合格ラインを上回っていたのに、不正の影響で不合格となった受験生５０人に対し、来年４月の入学を認める方針を固めたことが２４日、関係者への取材で分かった。
昨年の本来合格者１９人については、難しいとの意見が出ているもようだ。東京医大は５０人のうち入学希望が多数に上った場合は、来年の一般入試の募集人員減も検討。
<<
https://www.nishinippon.co.jp/sp/nnp/national/article/460101/

入学者数って100人程度のだよな？
半分が裏口ってこと？

**卵の名無しさん** · 2018/10/25(木) 19:31:40.10

treasure <- function(m=3,n=4,k=2){
y=1:(m*n)
(z=matrix(y,ncol=n,byrow=T))
(P=as.vector(z))
(Q=as.vector(t(z)))
PQ <- function(x){
p=q=numeric(k)
for(i in 1:k){
p[i]=which(P==x[i])
q[i]=which(Q==x[i])
}
min(p)-min(q)
}
tre=combn(m*n,k)
re=apply(tre,2,PQ)
print(c(P1st=sum(re<0),Q1st=sum(re>0),even=sum(re==0)))
invisible(list(P1st=tre[,which(re<0)],Q1st=tre[,which(re>0)],even=tre[,which(re==0)]))
}
t342=treasure(3,4,2)
print(matrix(LETTERS[t342$P1st],nrow=2),quote=F)
print(matrix(LETTERS[t342$Q1st],nrow=2),quote=F)
print(matrix(LETTERS[t342$even],nrow=2),quote=F)

t452=treasure(4,5,2)
print(matrix(LETTERS[t452$P1st],nrow=2),quote=F)
print(matrix(LETTERS[t452$Q1st],nrow=2),quote=F)
print(matrix(LETTERS[t452$even],nrow=2),quote=F)

**卵の名無しさん** · 2018/10/25(木) 19:32:05.84

縦nマス、横n+1マスのn(n+1)マスのうちランダムに選ばれた2マスにそれぞれ宝が眠っている。
縦1列を探し終えたらすぐ右の1列に移って宝を探していく方法をとるP君と、横1行を探し終えたらすぐ下の1行に移って宝を探していく方法をとるQ君が、同時に左上の地点から探索を開始した。
例えば、n=3の時はP君はAEIBFJCGKDHLの順で探す。Q君はABCDEFGHIJKの順で探すことになる。

ABCD
EFGH
I JK L

1つの地点を捜索するのにかかる時間は同じで、相手が1度探し終えた地点を重複して調べることも当然ある。
相手より先に宝を見つけた方を勝者とする。同時の場合は引き分けとする。
どちらの方が有利になるだろうか？

**卵の名無しさん** · 2018/10/25(木) 19:48:03.86

>>779
> t(sapply(1:20,treasure1))
P1st Q1st even
[1,] 0 0 1
[2,] 4 5 6
[3,] 26 27 13
[4,] 84 83 23
[5,] 203 197 35
[6,] 413 398 50
[7,] 751 722 67
[8,] 1259 1210 87
[9,] 1986 1910 109
[10,] 2986 2875 134
[11,] 4320 4165 161
[12,] 6054 5845 191
[13,] 8261 7987 223
[14,] 11019 10668 258
[15,] 14413 13972 295
[16,] 18533 17988 335
[17,] 23476 22812 377
[18,] 29344 28545 422
[19,] 36246 35295 469
[20,] 44296 43175 519

**卵の名無しさん** · 2018/10/25(木) 20:43:08.89

どれ、たまには算数の問題を投稿してやろう

http://seiga.nicovideo.jp/seiga/im6350659

**卵の名無しさん** · 2018/10/25(木) 21:18:10.86

>>780
4以上だと短軸探索が有利になるようにみえるな。

**卵の名無しさん** · 2018/10/25(木) 22:33:55.27

>>778

n=2

ABC
DEF

の場合

短軸方向探索Pが先に宝を発見する埋め方：４通り
> print(matrix(LETTERS[t232$P1st],nrow=2),quote=F)
[,1] [,2] [,3] [,4]
[1,] C D D E
[2,] D E F F

長軸方向探索Qが先に宝を発見する埋め方：５通り
> print(matrix(LETTERS[t232$Q1st],nrow=2),quote=F)
[,1] [,2] [,3] [,4] [,5]
[1,] B B B C C
[2,] C E F E F

同時に宝を発見する埋め方：６通り
> print(matrix(LETTERS[t232$even],nrow=2),quote=F)
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] A A A A A B
[2,] B C D E F D

**卵の名無しさん** · 2018/10/25(木) 23:25:11.32

算数の問題です

http://sp.seiga.nicovideo.jp/seiga/#!/im8552821

**卵の名無しさん** · 2018/10/25(木) 23:37:44.75

算数の問題です

http://sp.seiga.nicovideo.jp/seiga/#!/im8632142

**卵の名無しさん** · 2018/10/26(金) 00:02:56.11

>>779
数学板で２つの宝を先に発見した方を勝者にしたらどうなるかというレスがあったのでプログラムを改変してみた。

treasure2 <- function(m=3,n=4,k=2){
y=1:(m*n)
(z=matrix(y,ncol=n,byrow=T))
(P=as.vector(z))
(Q=as.vector(t(z)))
PQ <- function(x){
p=q=numeric(k)
for(i in 1:k){
p[i]=which(P==x[i])
q[i]=which(Q==x[i])
}
max(p)-max(q)
}
tre=combn(m*n,k)
re=apply(tre,2,PQ)
print(c(P1st=sum(re<0),Q1st=sum(re>0),even=sum(re==0)))
invisible(list(P1st=tre[,which(re<0)],Q1st=tre[,which(re>0)],even=tre[,which(re==0)]))
}

t342=treasure2(3,4,2)
#短軸方向探索Pが先に２つの宝を発見する埋め方
print(matrix(LETTERS[t342$P1st],nrow=2),quote=F)
#長軸方向探索Qが先に２つの宝を発見する埋め方
print(matrix(LETTERS[t342$Q1st],nrow=2),quote=F)
#同時に２つめの宝を発見する埋め方
print(matrix(LETTERS[t342$even],nrow=2),quote=F)

**卵の名無しさん** · 2018/10/26(金) 00:04:20.43

> t342=treasure2(3,4,2)
P1st Q1st even
27 26 13
> #短軸方向探索Pが先に２つの宝を発見する埋め方
> print(matrix(LETTERS[t342$P1st],nrow=2),quote=F)
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19]
[1,] A A A A A B B B B B C C C D E E E E F
[2,] E F I J K E F I J K I J K K F I J K I
[,20] [,21] [,22] [,23] [,24] [,25] [,26] [,27]
[1,] F F G G G I I J
[2,] J K I J K J K K
> #長軸方向探索Qが先に２つの宝を発見する埋め方
> print(matrix(LETTERS[t342$Q1st],nrow=2),quote=F)
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19]
[1,] A A A A A B B B B C C C C C D D D D D
[2,] B C D G H C D G H D E F G H E F G H I
[,20] [,21] [,22] [,23] [,24] [,25] [,26]
[1,] E E F F G H H
[2,] G H G H H I J
> #同時に２つめの宝を発見する埋め方
> print(matrix(LETTERS[t342$even],nrow=2),quote=F)
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
[1,] A B C D D E F G H H I J K
[2,] L L L J L L L L K L L L L

やはり、直観と違ってイーブンにはならないな。

**卵の名無しさん** · 2018/10/26(金) 00:14:13.73

縦５マス、横１０マスで宝が３マスに埋まっているときに

全部の宝を発見した方が勝者とすると、縦方向探索Pと横方向探索Qとして

勝者となる埋め方の場合の数は

> treasure2(5,10,3)
P1st Q1st even
9142 8832 1626

と算出できた。

部屋割りカウントできない底辺頭脳でもこれくらいは数えられるかな？

**卵の名無しさん** · 2018/10/26(金) 08:17:44.20

算数の問題です

http://seiga.nicovideo.jp/seiga/im8552821

中卒国試浪人の事務員には解けないだろうな

**卵の名無しさん** · 2018/10/26(金) 08:21:49.38

昼休みの薬の説明会はフィブラストだった。
熱傷診療ガイドライン(www.jsbi-burn.org/members/guideline/pdf/guideline2.pdf)に引用された論文としてこんなのが商用パンフに掲載されていた。
http://imagizer.imageshack.com/img923/3955/cIV6t9.jpg
早期群: n1=17 m1=36.2 sd1=21.5
後期群: n2=12 m2=56.1 sd2=30.1
で早期群の方が上皮化期間が有意(t検定でp<0.05)に短いと主張している。

生データがないので標本数・平均・標準偏差から検定してみる。

ttest2 = function(n1,m1,s1, n2,m2,s2){ # ni:標本数 mi:平均 si:標準偏差
n=n1+n2-2
u=((n1-1)*s1^2+(n2-1)*s2^2)/n
t=(m1-m2)/sqrt(u/n1+u/n2)
pe=2*pt(-abs(t),n)
t=(m1-m2)/sqrt(s1^2/n1 + s2^2/n2)
n=(s1^2/n1+s2^2/n2)^2 / ( (s1^2/n1)^2/(n1-1) + (s2^2/n2)^2/(n2-1))
pu=2*pt(-abs(t),n)
p.values=c(T.test=pe,Welch.test=pu)
return(p.values)
}

> ttest2(n1,m1,sd1, n2,m2,sd2)
T.test Welch.test
0.04701090 0.06461679
等分散を仮定したｔ検定でやっと0.05未満になるが、等分散を仮定しないWelchの方法だと有意差がでない。

**卵の名無しさん** · 2018/10/26(金) 08:22:44.14

http://imagizer.imageshack.com/img923/3955/cIV6t9.jpg
のグラフをみると後期群の最大値は外れ値に見える。

後期群が n2=12 m2=56.1 sd2=30.1
を満たすようなシミュレーションデータを作成して
そのデータの最大値を除いた標本群を修正後期群として
早期群と検定（等分散を仮定するｔ検定と仮定しないWelchの検定）してｐ値を出す。
これを1万回繰り返してｐ値がどんな値になるかやってみた。
その結果は、

等分散ｔ検定でｐ値は平均 0.105(95%CI 0.0916 0.116)で有意とは言えない。
Welochの検定でｐ値は平均 0.126(95%CI 0.114 0.129)で有意とは言えない。

シミュレーションデータによる検定のRスクリプトはここに書いた。
http://egg.2ch.net/test/read.cgi/hosp/1483665995/82

ガイドラインに引用されている論文すら、あんまり信用できないと思った。

“Statistics are like bikinis. What they reveal is suggestive, but what they conceal is vital.”

**卵の名無しさん** · 2018/10/26(金) 08:39:18.67

>>789
中卒で受けられる国試って何だろね？
その一文だけで頭が悪いのがわかる。

**卵の名無しさん** · 2018/10/26(金) 11:29:32.95

こういうのの実行結果を数学板に投稿すると
レスが返って来て勉強になるなぁ。

# ABCDEFGに対してなら
# BCDEFGAが一番勝率高い気がする

library(gtools)
n=7
k=2
perm=permutations(n,n)
Q=perm[1,]
np=nrow(perm)
p1st=numeric(np)
for(i in 1:np){
P=perm[i,]
tre=combn(n,k)
nt=ncol(tre)
re=numeric()
for(j in 1:nt){
re[j]=min(which(tre[1,j]==P),which(tre[2,j]==P))-
min(which(tre[1,j]==Q),which(tre[2,j]==Q))
}
p1st[i]=sum(re<0)
}
plot(p1st)
p1st[which.max(p1st)]
(p.max=which(p1st==15))
print(matrix(LETTERS[perm[p.max,]],ncol=7),quote=F)

#

**卵の名無しさん** · 2018/10/26(金) 11:48:48.32

宝を２個先に～い見つけた方が勝者になるとき

library(gtools)
n=8
k=2
perm=permutations(n,n)
Q=perm[1,]
np=nrow(perm)
p1st=numeric(np)
for(i in 1:np){
P=perm[i,]
tre=combn(n,k)
nt=ncol(tre)
re=numeric()
for(j in 1:nt){
re[j]=max(which(tre[1,j]==P),which(tre[2,j]==P))-
max(which(tre[1,j]==Q),which(tre[2,j]==Q))
}
p1st[i]=sum(re<0)
}
plot(p1st)
mv=p1st[which.max(p1st)]
(p.max=which(p1st==mv))
print(matrix(LETTERS[perm[p.max,]],ncol=8),quote=F)

**卵の名無しさん** · 2018/10/26(金) 12:35:57.72

数学板の、分からない問題はここに書いてね　スレで話題のテーマで問題を作ってみた。

ド底辺シリツ医大生が８人いてそのうち３人は裏口入学であることが判明している。
８人をA,B,C,D,E,F,G,Hとする。
A,B,C,D,E,F,G,Hの順に調査して３人みつかれば調査は終了とする。
A,B,C,D,E,F,G,Hの順の調査と比べて早くみつかる可能性の探索順を列挙せよ。

[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8]
[1,] A B D E F G H C
[2,] A C D E F G H B
[3,] A D B E F G H C
[4,] A D C E F G H B
[5,] B A D E F G H C
[6,] B C D E F G H A
[7,] B D A E F G H C
[8,] B D C E F G H A
[9,] C A D E F G H B
[10,] C B D E F G H A
[11,] C D A E F G H B
[12,] C D B E F G H A
[13,] D A B E F G H C
[14,] D A C E F G H B
[15,] D B A E F G H C
[16,] D B C E F G H A
[17,] D C A E F G H B
[18,] D C B E F G H A

おい、ド底辺頭脳よ、これであっているか検証してみろ。

**卵の名無しさん** · 2018/10/26(金) 12:56:44.55

treasure <- function(m=3,n=4,k=2){
y=1:(m*n)
(z=matrix(y,ncol=n,byrow=T))
(P=as.vector(z))
(Q=as.vector(t(z)))
PQ <- function(x){
p=q=numeric(k)
for(i in 1:k){
p[i]=which(P==x[i])
q[i]=which(Q==x[i])
}
min(p)-min(q)
}
tre=combn(m*n,k)
re=apply(tre,2,PQ)
return(c(短軸有利=sum(re<0),長軸有利=sum(re>0),同等=sum(re==0)))
}
sapply(1:12,function(k) treasure(3,4,k))
sapply(1:20,function(k) treasure(4,5,k))

**卵の名無しさん** · 2018/10/26(金) 13:32:10.72

>789　や他の算数の問題の答え　まだぁ

宅建国試浪人の中卒事務員さん

**卵の名無しさん** · 2018/10/26(金) 15:19:53.27

>>797
算数の問題ですらないじゃん、

おまえは問題作成すらできないドアホなのか？

最近のトピックをつかったこういう問題面白いだろ？

答える頭はお前にはないだろうけど。

ド底辺シリツ医大生が８人いてそのうち３人は裏口入学であることが判明している。
８人をA,B,C,D,E,F,G,Hとする。
A,B,C,D,E,F,G,Hの順に調査して３人みつかれば調査は終了とする。
A,B,C,D,E,F,G,Hの順の調査と比べて早くみつかる可能性の探索順を列挙せよ。

ある医大で合格率の男女比が1.2で男子優位という結果だったという。
定員100で男子800人女子200人が受験して合格率の男女比が
1.2であったときに統計的には有意差があると言えるか？

**卵の名無しさん** · 2018/10/26(金) 16:34:55.99

国立大学をでていれば、これくらいの計算ができて当然。

できない椰子はド底辺シリツ医大卒と同じ頭脳レベルといえる。

東京医大、本来合格者入学許可へ　今年の受験生５０人

2018年10月25日 02時06分
　東京医科大＝８月、東京都新宿区
　東京医科大が今年の入試で本来合格ラインを上回っていたのに、不正の影響で不合格となった受験生５０人に対し、来年４月の入学を認める方針を固めたことが２４日、関係者への取材で分かった。
昨年の本来合格者１９人については、難しいとの意見が出ているもようだ。東京医大は５０人のうち入学希望が多数に上った場合は、来年の一般入試の募集人員減も検討。

https://www.nishinippon.co.jp/sp/nnp/national/article/460101/

https://www.tokyo-med.ac.jp/med/enrollment.htmlによると

学年第1学年第2学年

在学者数 133 113

昨年入学者の留年者や退学者が０として、

大学が公式認定した裏口入学者が少なくとも今年は１３３人中５０人、昨年が１１３人中１９人ということになる。

裏口入学率の期待値、最頻値、およびその９５％信頼区間を求めよ。

**卵の名無しさん** · 2018/10/26(金) 16:36:52.25

library(rjags)
y=c(50,19)
N=c(133,113)
Ntotal=length(y)
a=1
b=1
dataList=list(y=y,N=N,Ntotal=Ntotal,a=a,b=b)
# JAGS model
modelString ="
model{
for(i in 1:Ntotal){
y[i] ~ dbin(theta,N[i])
}
theta ~ dbeta(a,b)
}
"
writeLines(modelString,'TEMPmodel.txt')
jagsModel=jags.model('TEMPmodel.txt',data=dataList)
codaSamples=coda.samples(jagsModel,var=c('theta'),n.iter=20000,na.rm=TRUE)
summary(codaSamples)
js=as.vector(as.matrix(codaSamples))

BEST::plotPost(js,,xlab="裏口確率")
BEST::plotPost(js,showMode = TRUE,xlab="裏口確率")

**卵の名無しさん** · 2018/10/26(金) 16:43:26.58

>>799
昨年と今年で裏口入学率に違いがあったといえるか？

裏口入学率の差の期待値とその９５％信頼区間を求めよ

**卵の名無しさん** · 2018/10/26(金) 16:44:57.64

>>801

解法はここに置いた

https://rio2016.5ch.net/test/read.cgi/math/1490734993/208

**卵の名無しさん** · 2018/10/26(金) 16:56:31.17

>>800
裏口入学率の事前分布は一様分布よりJefferey分布の方がいいかなと思ってやってみたが、殆ど数値に影響なかったな。

**卵の名無しさん** · 2018/10/26(金) 20:01:23.00

>789　や他の算数の問題の答え　まだぁ

宅建国試浪人の中卒事務員さん