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

1 名前：卵の名無しさん mailto:sage [2020/06/15(月) 19:09:47 ID:jnY7bdq5.net] タイトルまんまです すなわち、タイトルの正しさを検証するスレです。 ※前スレ 底辺私立医大を卒業した頭悪いよね？ egg.2ch.net/test/read.cgi/hosp/1472273747/ 底辺私立医大を卒業した医者って頭悪いよね？ Part2 egg.2ch.net/test/read.cgi/hosp/1520807375/ 底辺私立医大を卒業した医者って頭悪いよね？ Part3 egg.2ch.net/test/read.cgi/hosp/1522967067/ 底辺私立医大を卒業した医者って頭悪いよね？ Part4 egg.2ch.net/test/read.cgi/hosp/1526677794/ 底辺私立医大を卒業した医者って頭悪いよね？ Part5 egg.2ch.net/test/read.cgi/hosp/1529666055/ 底辺私立医大を卒業した医者って頭悪いよね？ Part6 egg.5ch.net/test/read.cgi/hosp/1531305684/ 底辺私立医大を卒業した医者って頭悪いよね？ Part7 egg.2ch.net/test/read.cgi/hosp/1532472024/ 底辺私立医大を卒業した医者って頭悪いよね？ Part8 egg.2ch.net/test/read.cgi/hosp/1533510399/ 底辺私立医大を卒業した医者って頭悪いよね？ Part9 egg.2ch.net/test/read.cgi/hosp/1535547275/ 底辺私立医大を卒業した医者って頭悪いよね？ Part10 https://egg.2ch.net/test/read.cgi/hosp/1537867846/ 底辺私立医大を卒業した医者って頭悪いよね？ Part11 https://egg.5ch.net/test/read.cgi/hosp/1544320163/ 底辺私立医大を卒業した医者って頭悪いよね？ Part12 https://egg.5ch.net/test/read.cgi/hosp/1550969602/ 底辺私立医大を卒業した医者って頭悪いよね？ Part13 https://egg.5ch.net/test/read.cgi/hosp/1575242106/ 底辺私立医大を卒業した医者って頭悪いよね？ Part14 https://egg.5ch.net/test/read.cgi/hosp/1578616997/ 底辺私立医大を卒業した医者って頭悪いよね？ Part15 https://egg.5ch.net/test/read.cgi/hosp/1584050953/ 底辺私立医大を卒業した医者って頭悪いよね？ Part16 https://egg.5ch.net/test/read.cgi/hosp/1587907116/ 901 名前：卵の名無しさん mailto:sage [2020/11/08(日) 13:06:06.92 ID:oC4O4pgk.net] ある病院に内科医A,B,C、外科医D,Eがいて１週間(日〜土）の当直と呼び出し待機の割り当てをする。 以下の条件を満たすように割りあてる。 (1) 誰も少なくとも1回は日〜土の間で当直および待機に割り当てられる (2)　誰も続けて勤務（当直または待機）してはならない。 (3)　誰においても１週間の当直総数および待機総数の上限はどちらも２日である (4)　内科医が当直のときは待機は外科医、外科医が当直の時は内科医が待機する。 (5) 当直は２日以上の間隔を空ける 今週の勤務割当は次の通りとする。 日 月 火 水 木 金 土 当直 A B D E C B D 待機 D E C B D E A 週を跨いで(1)〜(5)の条件を満たす来週の割りあて方は何通りあるか？ 902 名前：卵の名無しさん mailto:sage [2020/11/08(日) 13:47:42.36 ID:3UafLdM8.net] 割り当ては132通り 例 > print(rei,quote=F) 日 月 火 水 木 金 土 当直 E B A C B D E 待機 C D E D E A B > print(rei,quote=F) 日 月 火 水 木 金 土 当直 C B E D A B E 待機 E D A B E D C 903 名前：卵の名無しさん mailto:sage [2020/11/08(日) 17:32:58.85 ID:3UafLdM8.net] " 1から6までの数字を使って5ケタの数字を作ります。 同じ数字を何回使ってもいいとすると、6の倍数は1296通りある。 小さい順に並べたときに1000番目にくる数字を述べよ。 " rm(list=ls()) source('toolmini.R') library(gtools) pm=gtools::permutations(6,5,rep=TRUE) f <- function(x) sum(x * 10^(4:0))%%6==0 y=pm[apply(pm,1,f),] head(y,10) ; tail(y,10) cat(paste(y[777,],collapse='')) # 小さい順に並べたときにn番目にくる数字 order2num <- function(n){ # cat(paste(y[n,],collapse=''),'\n') a=10*as.numeric(paste(dec2nw(n,6),collapse=''))+11100 a+(1:6)[(a+1:6)%%6==0] } n=1:(6^4-1) plot(n,sapply(n,order2num)) order2num(333) order2num(777) 904 名前：卵の名無しさん mailto:sage [2020/11/08(日) 20:10:46.07 ID:3UafLdM8.net] ある病院に内科医A,B,C、外科医D,Eがいて１週間(日〜土）の当直と呼び出し待機の割り当てをする。 以下の条件を満たすように割りあてる。 (1) 誰も少なくとも1回は日〜土の間で当直および待機に割り当てられる (2)　誰も続けて勤務（当直または待機）してはならない。 (3)　誰においても１週間の当直総数および待機総数の上限はどちらも２日である (4)　内科医が当直のときは待機は外科医、外科医が当直の時は内科医が待機する。 (5) 当直は２日以上の間隔を空ける 今週の勤務割当は次の通りとする。 日 月 火 水 木 金 土 当直 A B D E C B D 待機 D E C B D E A 週を跨いで(1)〜(5)の条件を満たす来週の割りあて方は何通りあるか？ library(gtools) pm=permutations(5,7,rep=T) # 当直可能配列 f <- function(x){ # (1) & 連続当直不可 & (3) & (5)（週跨ぎ） int=NULL for(i in 1:5){ a=which(i==x) # a : i の当直日 int[i]=ifelse(length(a)==1, 905 名前：TRUE,diff(a)>2) # 当直間隔が２を越えるか？ } all(1:5 %in% x) & all(diff(x)!=0) & all(table(x)<=2) & all(int) & all(x[1]!=c(1,2,4)) & # 日曜日A,B,D当直不可 x[2]!=4 # 月曜日D当直不可 } i=apply(pm,1,f) sum(i) y=pm[i,]　# 当直可能配列 [] [ここ壊れてます] 906 名前：卵の名無しさん mailto:sage [2020/11/08(日) 20:11:10.62 ID:3UafLdM8.net] # 待機可能配列 f1 <- function(x){ all(1:5 %in% x)　& all(diff(x)!=0) & all(table(x)<=2) & # 最低1回　& 連続待機不可 & 待機日総数<=2 all(x[1]!=c(1,4)) # 日曜日A,D待機不可 } y1=pm[apply(pm,1,f1),] g <- function(x){ # x:当直候補配列 idx=NULL for(i in 1:nrow(y1)){ # y1=待機候補配列 if(all(x %in% 1:3 == !(y1[i,] %in% 1:3))){ #4 当直医属性!=待機医属性 idx=append(idx,i) } } re=NULL　# 条件を満たす待機配列のy1での行index for(i in idx){ cow=FALSE # 連続勤務continual work? for(j in 1:5){ a=which(j==x) # 当直曜日 b=which(j==y1[i,])# 待機曜日 if(any(diff(sort(c(a,b)))<=1)){ # 連続勤務ならbreak cow=TRUE break } } if(!cow) re=append(re,i) # 連続勤務でなければ採用 } re } 907 名前：卵の名無しさん mailto:sage [2020/11/08(日) 20:11:49.71 ID:3UafLdM8.net] gg <- function(x,print=TRUE){ re=g(x) if(!is.null(re)){ n=length(re) ans=vector(mode='list',length=n) for(i in 1:n){ rei=rbind(LETTERS[x],LETTERS[y1[re[i],]]) colnames(rei)=c('日','月','火','水','木','金','土') rownames(rei)=c('当直','待機') ans[[i]]=rei if(print){print(rei,quote=F) cat('\n')} } invisible(ans) } } zz=apply(y,1,gg) nzz=length(zz) zzz=NULL for(i in 1:nzz){ if(!is.null(zz[[i]])) zzz=c(zzz,zz[[i]]) } print(head(zzz),quote=F) print(zzz,quote=F) 908 名前：卵の名無しさん mailto:sage [2020/11/08(日) 20:12:12.93 ID:3UafLdM8.net] > print(head(zzz),quote=F) [[1]] 日 月 火 水 木 金 土 当直 C A E B A D E 待機 E D C D E B A [[2]] 日 月 火 水 木 金 土 当直 C A E B C D E 待機 E D C D E A B [[3]] 日 月 火 水 木 金 土 当直 C A E B C D E 待機 E D C D E B A [[4]] 日 月 火 水 木 金 土 当直 C A E C B D E 待機 E D B D E A C [[5]] 日 月 火 水 木 金 土 当直 C A E C B D E 待機 E D B D E C A [[6]] 日 月 火 水 木 金 土 当直 C A E D A B E 待機 E D B C E D A 909 名前：卵の名無しさん mailto:sage [2020/11/08(日) 20:12:38.07 ID:3UafLdM8.net] [[101]] 日 月 火 水 木 金 土 当直 E C B D C A E 待機 B D E A E D C [[102]] 日 月 火 水 木 金 土 当直 E C B D E A B 待機 B D E A C D E [[103]] 日 月 火 水 木 金 土 当直 E C B D E A C 待機 B D E A C D E [[104]] 日 月 火 水 木 金 土 当直 E C B D E B A 待機 B D E A C D E [[105]] 日 月 火 水 木 金 土 当直 E C B D E B A 待機 B D E C A D E [[106]] 日 月 火 水 木 金 土 当直 E C B D E C A 待機 B D E C A D E 910 名前：卵の名無しさん mailto:sage [2020/11/08(日) 20:41:39.89 ID:oC4O4pgk.net] >>837 これは間違い、106通りだと思う。 バグがなければだけど 911 名前：卵の名無しさん mailto:sage [2020/11/09(月) 05:32:00.75 ID:NqCZIoxd.net] >>839 更に (6) ２週続けて週２回の当直はしてはならない。 を付け加えてみよう。 912 名前：卵の名無しさん mailto:sage [2020/11/09(月) 05:40:31.38 ID:NqCZIoxd.net] ある病院に内科医A,B,C、外科医D,Eがいて１週間(日〜土）の当直と呼び出し待機の割り当てをする。 以下の条件を満たすように割りあてる。 (1)　誰も少なくとも1回は日〜土の間で当直および待機に割り当てられる (2)　誰も続けて勤務（当直または待機）してはならない。 (3)　誰においても１週間の当直総数および待機総数の上限はどちらも２日である (4)　内科医が当直のときは待機は外科医、外科医が当直の時は内科医が待機する。 (5)　当直は２日以上の間隔を空ける (6)　２週続けて週２回の当直はしてはならない。 今週の勤務割当は次の通りとする。 日 月 火 水 木 金 土 当直 A B D E C B D 待機 D E C B D E A 週を跨いで(1)〜(6)の条件を満たす来週の割りあて方は何通りあるか？ library(gtools) pm=permutations(5,7,rep=T) # 当直可能配列 f <- function(x){ # (1) & 連続当直不可 & (3) & (5)（週跨ぎ） int=NULL for(i in 1:5){ a=which(i==x) # a : i の当直日 int[i]=ifelse(length(a)==1,TRUE,diff(a)>2) # 当直間隔が２を越えるか？ } all(1:5 %in% x) & all(diff(x)!=0) & all(table(x)<=2) & all(int) & all(x[1]!=c(1,2,4)) & # 日曜日A,B,D当直不可 x[2]!=4 &# 月曜日D当直不可 sum(x==2)<2 & sum(x==4) # B,D の当直は2回未満 } i=apply(pm,1,f) sum(i) y=pm[i,]　# 当直可能配列 913 名前：卵の名無しさん mailto:sage [2020/11/09(月) 08:25:50.50 ID:Wee9Qn0s.net] lwa=rbind(c(1,2,4,5,3,2,4),c(4,5,3,2,4,5,1)) # last week allocation rei=rbind(LETTERS[c(1,2,4,5,3,2,4)],LETTERS[c(4,5,3,2,4,5,1)]) colnames(rei)=c('日','月','火','水','木','金','土') rownames(rei)=c('当直','待機') print(rei,quote=F) 914 名前：卵の名無しさん mailto:sage [2020/11/09(月) 14:17:47.89 ID:Z0+EXym3.net] 米国が移民政策成功してるのって弱者が病気になろうと 気にせずゴミのように放置できるから社会コストが跳ね上がらないからだろう。 日本は多民族共生できるほど民度が高くないんだよなぁ。 日本だと気の毒だからと国保にタカらせ集中治療室まで使わせるからパンクする 915 名前：卵の名無しさん mailto:sage [2020/11/10(火) 05:18:34.17 ID:WoaXhYeo.net] https://i.imgur.com/v51hvTf.jpg 916 名前：卵の名無しさん mailto:sage [2020/11/10(火) 14:32:59.72 ID:WoaXhYeo.net] 旅行にいく余裕もない人の税金を横取りして旅にでる強盗キャンペーン。 それでコロナが広がるならまさに死神政策。 917 名前：卵の名無しさん mailto:sage [2020/11/11(水) 19:37:44.50 ID:G3dq4Z8e.net] 緊急提言ではさらに、これらの対策を取った上でも、感染状況ごとに分類した四つのステージのうち、感染者が急増している「ステージ3」相当以上と判断された場合には 「社会経済活動に一定の制約を求めるような強い対策を行う必要があることから、国民が一丸となって対策を進めていく必要がある」とした。 https://mainichi.jp/articles/20201109/k00/00m/040/249000c >国民が一丸となって 一丸となったらクラスターが発生ｗ 918 名前：卵の名無しさん mailto:sage [2020/11/13(金) 09:56:44.03 ID:wZuIWy8w.net] 理論を認めてもその通りに行動できない暗愚が学術会議の任命権者ておかしくない？ https://i.imgur.com/1vY3hqw.jpg https://youtu.be/fVosZS6 919 名前：-UJw?t=498 8分18秒 [] [ここ壊れてます] 920 名前：卵の名無しさん mailto:sage [2020/11/13(金) 21:20:33.70 ID:Sh6hhYou.net] Keltonの新著を読み始めた。 The Deficit Myth: Modern Monetary Theory and the Birth of the People's Economy (English Edition) Kindle版 この記述は意外だったな。 But even Senator Bernie Sanders has echoed Reagan, saying, “I am concerned about the debt. It’s not something we should be leaving to our kids and our grandchildren.” 921 名前：卵の名無しさん mailto:sage [2020/11/13(金) 21:22:53.38 ID:Sh6hhYou.net] The fifth myth is that deficits make the United States dependent on foreigners. This myth would have us believe that countries like China and Japan have enormous leverage over us because they hold large quantities of US debt. 米国債を売ろうという誘惑に言及した総理は不審死したよあなぁ。 922 名前：卵の名無しさん mailto:sage [2020/11/14(土) 17:11:40.19 ID:Fnd8Onk+.net] zを-1でない複素数として ω= （z+2i) / (z+1) とおく。 |z|=1のとき、偏角が3π/4となるωを求めよ。 options(digits=5) f <- function(x){ z=cos(x)+1i*sin(x) omega=(z+2i)/(z+1) Arg(omega) } curve(f,0,2*pi) ; abline(h=3/4*pi,lty=3) low=optimise(f,c(4,5),maximum = T)$maximum (x=uniroot(function(x)f(x)-3/4*pi,c(low,5))$root) (z=cos(x)+1i*sin(x)) (omega=(z+2i)/(z+1)) 923 名前：卵の名無しさん mailto:sage [2020/11/14(土) 22:56:05.54 ID:Fnd8Onk+.net] z=a+bi ω=(z+2i)/(z+1) =(a+bi+2i)/(a+bi+1) 分母子にa+1-biをかけると 分子=(a+bi+2i)(a+1-ib) 分母=(a+1)^2+b^2 ωの偏角が3/4πということは ωの虚部/ωの実部=tan(3/4π)=-1 分母は実数なので約分されて 分子の虚部/分子の実部=-1 これを変形して 分子の虚部＋分子の実部=0 複素数α+iβをα+βに変形するにはi=1と置換すればいいので 分子の式(a+bi+2i)(a+1-ib)でi=1と置換して (a+b+2)(a+1-b)=0 ∴ a+b+2=0 またはa+1-b=0 |z|=1ゆえa^2+b^2=1 a+b+2=0 かつ　a^2+b^2=1に実数解なし a+1-b=0　かつ　a^2+b^2=1 を解くと(a,b)=(-1,0),(0,1) z≠-1なのでa=0,b=1が残る 924 名前：卵の名無しさん mailto:sage [2020/11/15(日) 01:51:01.75 ID:L0SCzfIz.net] z=a+bi (a,bは実数) ω=(z+2i)/(z+1) 分母子に(a-bi+1)をかけると 分子=(z+2i)(a-bi+1)=(a+bi+2i)(a-bi+1)= (a^2+a+b^2+2b) +i(2a+b+2) 分母=(z+1)(a-bi+1)=(a+bi+1)(a-bi+1)=(a+1)^2+b^2 (実数) ωの虚部/ωの実部=分子の虚部/分子の実部=tan(4π/3)=-1 これを変形すれば、 分子の実部+分子の虚部=0 複素数α+iβをα+βに変形するにはこの形でi=1と置換すればいいので 分子のiを1に置換して a^2+a+b^2+2b +(2a+b+2) = a^2+b^2+3a+3b+2=0 |z|=1ゆえa^2+b^2=1なので 連立方程式 a^2+b^2+3a+3b+2=0 a^2+b^2=1 を解くと a=-1,b=0 または a=0,b=-1 z≠-1ゆえa=-1,b=0は不適。 候補はa=0,b=-1 z=a+bi=-i このとき ω=(z+2i)/(z+1)=i/(1-i) = -0.5+0.5i 偏角=atan(0.5/(-0.5))=(3/4)π 925 名前：卵の名無しさん mailto:sage [2020/11/15(日) 09:54:17.96 ID:L0SCzfIz.net] △ABC (∠A＝50°, ∠B＝70°) をAを原点、Bをx軸上に配置する。 △ABCの内部の点をPとしてPA=2 PB=3, PC=4 のとき、 (1)B,C,Pの座標を求めよ。 (2)△PAB,△PBC,△PCAの面積を求めよ > re$B [1] 4.9368+0i > re$C [1] 3.4432+4.1035i > re$P [1] 1.9619+0.3 926 名前：879i > with(re,ABC2S(P,A,B)) #△PAB [1] 0.95738 > with(re,ABC2S(P,B,C)) #△PBC [1] 5.8139 > with(re,ABC2S(P,C,A)) #△PCA [1] 3.3576 [] [ここ壊れてます] 927 名前：卵の名無しさん mailto:sage [2020/11/15(日) 09:56:10.43 ID:L0SCzfIz.net] > re$B [1] 4.9368+0i > re$C [1] 3.4432+4.1035i > re$P [1] 1.9619+0.3879i > with(re,ABC2S(P,A,B)) #△PAB [1] 0.95738 > with(re,ABC2S(P,B,C)) #△PBC [1] 5.8139 > with(re,ABC2S(P,C,A)) #△PCA [1] 3.3576 > re$S　# △ABC [1] 10.129 https://i.imgur.com/PFt8IBP.png 928 名前：卵の名無しさん mailto:sage [2020/11/15(日) 14:55:20.09 ID:L0SCzfIz.net] [1] 4.652337六角形AB,CD,EFの辺の長さをn,BC,DE,FAの長さをm、 ∠A=degA° ∠C=degC°　∠E=degE°、とする。 m=1,n=2, degA=60,degC=90,degE=120のとき この六角形の面積を求めなさい。 作図するプログラムを書くのが面白かっただけ。 https://i.imgur.com/m045k8H.png > Hexagon(1,2,60,90,120) [1] 4.652337 929 名前：卵の名無しさん mailto:sage [2020/11/15(日) 18:05:37.35 ID:L0SCzfIz.net] よくも悪くも、中国共産党の対日外交には虚々実々の駆け引きというか、ある程度の枠がはまっているのに対して、かりに将来、中国の民主化が進んで「反日デモ」的な大衆感情に熱狂的に支持されて議会の多数派を占めるような政権が誕生した場合、日本に対する姿勢が今までと比較にならないくらい強硬になるのではないかという不安は、多くの日本人が感じたことがあるのではないか 中国化する日本　より 930 名前：卵の名無しさん mailto:sage [2020/11/15(日) 20:41:06.93 ID:L0SCzfIz.net] >>859 > re=DV2A(50,70,2,3,4) > attach(re) > A=0i > B [1] 4.936762+0i > C [1] 3.443221+4.103471i > P [1] 1.961936+0.387857i > ABC2S(P,A,B) #△PAB [1] 0.9573792 > ABC2S(P,B,C) #△PBC [1] 5.813914 > ABC2S(P,C,A) #△PCA [1] 3.357635 > S　# △ABC [1] 10.12893 > D=intsect(A,P,B,C) ; pt(D,'D') > seg(P,D,col=2,lty=3) > E=intsect(A,C,P,B) ; pt(E,'E') > seg(P,E,col=2,lty=3) > F=intsect(A,B,C,P) ; pt(F,'F') > seg(P,F,col=2,lty=3) > D [1] 4.605388+0.910444i > E [1] 0.4868298+0.5801812i > F [1] 1.807311+0i > detach(re) 931 名前：卵の名無しさん mailto:sage [2020/11/16(月) 15:10:00.97 ID:FmGCMvZ1.net] こういう毛針に騙されたんだなぁ https://livedoor.blogimg.jp/rock1963roll/imgs/f/b/fb36dca5-s.jpg 932 名前：卵の名無しさん mailto:sage [2020/11/17(火) 11:39:26.43 ID:zTNDOLQn.net] # 小数点付きの数numをN進法で表示する(62進法まで対応 0-9,a-z,A-Z) rm(list=ls()) num=1/3;N=3;kmin=5 decN <- function(num, N, kmin = 5){ # kmin：最小小数点後桁 int=floor(num) r=int%%N q=int%/%N while(q > 0){ r=append(q%%N,r) q=q%/%N } # rに整数のN進法表示数列を格納 k=max(nchar(num)-nchar(floor(num))-1,kmin) # 同長もしくはkminの長さの小数表示 a=numeric(k) x=round(num-floor(num),k) # e.g. 7.28-floor(7.28)!=0.280に対応 for(i in 1:k){ y=round(x*N,k) # e.g. 0.728*5-3.64 !=0 に対応 a[i]=floor(y) x=y-a[i] # r . a[1] a[2] a[3] ... a[k] } b=list(integer=r,decimal=a,num=sum(c(int,a)*(1/N)^(0:k))) fig=c(0:9,letters,LETTERS)[1:N] if(N<=62){ # Nが62以下は数値として表示 cat(paste(fig[b$integer+1],sep='',collapse=''), '.',paste(fig[b$decimal+1],sep='',collapse=''),sep='') cat('\n') } else{ # Nが63以上は整数部分と小数部分を数列で表示 print(b[1:2]) } invisible(b) # b$num:検証用 } 933 名前：卵の名無しさん mailto:sage [2020/11/17(火) 11:39:34.53 ID:zTNDOLQn.net] del=sapply(1:16, function(x) decN(0.1,x,16)) decN(pi,16) decN(0.1,13,10) # １３進法 0 1 2 3 4 5 6 7 8 9 A J Q options(digits=22) decK <- function(PI=pi,K=10){ x=decN(PI,13,kmin=K) K <- function(x){ if(x==0) return('0') else{ c(as.character(1:9),'A','J','Q')[x] } } K=Vectorize(K) re=paste(c(K(x$integer),'.',K(x$decimal)),collapse = '') print(re,quote=F) invisible(re) } decK(pi) decK(exp(1)) decK(sqrt(2)) 934 名前：卵の名無しさん mailto:sage [2020/11/17(火) 18:07:44.62 ID:annbVQoU.net] １４進法 で使う数字を小さい方から0 1 2 3 4 5 6 7 8 9 A J Q Kとするときπ、ネイピア数、√２を１４進法で小数１０桁まで表示せよ。 935 名前：卵の名無しさん mailto:sage [2020/11/17(火) 18:08:49.66 ID:annbVQoU.net] 62進法 で使う数字を小さい方から0 1 2 3 4 5 6 7 8 9 a b c ... x y z A B C ... X Y Z とするときπ、ネイピア数、√２を１４進法で小数１０桁まで表示せよ。 936 名前：卵の名無しさん mailto:sage [2020/11/17(火) 21:03:35.41 ID:annbVQoU.net] > decN(pi,62) 3.8MhuCIRym8egOl > decN(exp(1),62) 2.Ix4FDGbb8J4wtH > decN(sqrt(2),62) 1.pGeGLWSyyECdM 937 名前：卵の名無しさん mailto:sage [2020/11/17(火) 23:26:12.97 ID:annbVQoU.net] 0と1の2個の数を使いマイナス二進法で1から10を数えよ　（マイクロソフトの入社試験） 答と解説は https://i.imgur.com/F5AdX1Q.jpg https://www.youtube.com/watch?v=CTzw-zi7Zro 20まで計算させてみた ans 1 1 2 110 3 111 4 100 5 101 6 11010 7 11011 8 11000 9 11001 10 11110 11 11111 12 11100 13 11101 14 10010 15 10011 16 10000 17 10001 18 10110 19 10111 20 10100 それを面倒にした問題 (1) 2021をマイナス二進法で表示せよ。 (2) 3, 33, 333, 3333を各々マイナス三進法で表示せよ。 938 名前：卵の名無しさん mailto:sage [2020/11/17(火) 23:26:38.75 ID:annbVQoU.net] MS <- function( # Minus N-ary System num, N=-2, d=16, print=TRUE ){ base=N^((d-1):0) pm=expand.grid(replicate(d,0:(-N-1),simplify = FALSE)) colnames(pm)=1:d f <- function(x) sum(x*base) == num re=pm[apply(pm,1,f),] i=1 flg <- re[1]==0 while(flg){ i=i+1 flg <- re[i]==0 } ans=paste(as.character(re[i:d]),collapse='') if(print) print(ans,quote=F) invisible(ans) } data.frame(ans=sapply(1:20,function(n)MS(n,p=F))) MS(2021,-2) MS(3,-3,d=3) MS(33,-3,d=5) MS(333,-3,d=7) MS(3333,-3,d=9) 939 名前：卵の名無しさん mailto:sage [2020/11/18(水) 00:40:03.43 ID:vfLJTrWt.net] 3を0,1,2の3つの数字をつかったマイナス三進法で表すと3 = 1*(-3)^2 + 2*(-3) + 0*(-3)^0なので　120と表せる。 1から10までをマイナス三進法で表示すると 1 1 2 2 3 120 4 121 5 122 6 110 7 111 8 112 9 100 10 101 問題　：　円周率をマイナス三進法で小数１０桁まで表示せよ。 940 名前：卵の名無しさん mailto:sage [2020/11/18(水) 01:43:54.17 ID:vfLJTrWt.net] MS <- function( # Minus N-ary System num, N=-2, d=16, print=TRUE ){ base=N^((d-1):0) pm=expand.grid(replicate(d,0:(-N-1),simplify = FALSE)) colnames(pm)=1:d f <- function(x) sum(x*base) == num re=pm[apply(pm,1,f),] i=1 flg <- re[1]==0 while(flg){ i=i+1 flg <- re[i]==0 } ans=paste(as.character(re[i:d]),collapse='') if(print) print(ans,quote=F) invisible(ans) } piN <- function(N=-2){ # pi for minus N-ary system Pi=unlist(strsplit(as.character(pi),'.',fixed=T)) seisu=MS(as.numeric(Pi[1]),N,d=3,print=F) # 3(mod10) = 120(mod-3) pm=expand.grid(replicate(10,0:(abs(N)-1),simplify = F)) base=(1/N)^(1:10) f <- function(x) (pi-floor(pi)- sum(x*base))^2 dif=apply(pm,1,f) shosu=paste(as.character(pm[which.min(dif),]),collapse='') print(paste0(seisu,'.',shosu),quote=F) } 941 名前：卵の名無しさん mailto:sage [2020/11/18(水) 01:44:08.48 ID:vfLJTrWt.net] > piN(-2) [1] 111.0110010001 > piN(-3) [1] 120.0220210200 > piN(-5) [1] 3 942 名前：.0434333021 [] [ここ壊れてます] 943 名前：卵の名無しさん mailto:sage [2020/11/18(水) 12:25:29.26 ID:EQ55+bnh.net] # e+0*e^0+1*e^(-1)+0*e^(-2)+1*e^(-3) # [1] 3.135948 n=10 e=exp(1) 1*e^1+0*e^0 < pi pi < 1*e^1 + 1*e^1 a=numeric() b=numeric() b[1]=e+0*e^0+e^(-1) a[1]=as.numeric(b[1]<pi) for(i in 1:n){ a[i+1]=as.numeric((b[i]+e^(-i-1))<pi) b[i+1]=b[i]+a[i+1]*e^(-i-1) } paste(as.character(a),collpse="") 944 名前：卵の名無しさん mailto:sage [2020/11/18(水) 21:59:58.98 ID:+VMfzcIp.net] 武漢からの飛行機は乗客降ろさず機内で検査。 海外からの入国は全て隔離検査。 電車バスはマスク無しは乗車拒否かつ検温を改札で実施。 店はマスク無しは入店拒否。 マスクはコンビニで無料配布かつスマホで一人一人の枚数管理するアプリを政府が開発。 学校は閉鎖しオンライン授業に切り替え。 ざっと台湾の対策はこんなもの。 945 名前：卵の名無しさん mailto:sage [2020/11/19(木) 08:20:55.90 ID:LJQTN5dj.net] # (x-r*cos(θ))^2 + (y-r*sin(θ))^2=r^2 # y = r*sin(θ) + sqrt(r^2*(sin(θ))^2 + 2*r*x -x^2) area <- function(r,θ){ integrate(function(x) r*sin(θ) + sqrt(r^2*(sin(θ))^2 + 2*r*x -x^2),0,1)$value } 946 名前：卵の名無しさん mailto:sage [2020/11/19(木) 08:27:00.14 ID:LJQTN5dj.net] # (x-r*cos(θ))^2 + (y-r*sin(θ))^2=r^2 # y = r*sin(θ) + sqrt(r^2*(sin(θ))^2 + 2*r*x -x^2) area <- function(r,θ){ f = function(x) r*sin(θ) + sqrt(r^2*(sin(θ))^2 + 2*r*x -x^2) O=c(0,0) P=c(r*cos(θ),r*sin(θ)) Q=c(1,f(1)) integrate(f,0,1)$value } 947 名前：卵の名無しさん mailto:sage [2020/11/19(木) 13:16:21.52 ID:AV9OAvON.net] rm(list=ls()) source('toolmini.R') plot(NULL,xlim=c(0,1.5),ylim=c(0,1.5),asp=F,ann=F,axes=F) abline(h=0,v=0,col=8) r=1.015 ci <- function(x,R=r) sqrt(R^2-x^2) cir(0,0,r,lty=4) x=seq(ci(1),r,0.01) lines(x,ci(x),type='h',col=4) Q=ci(1)+1i seg(Q,ci(1)+0i) seg(Q,Q+1) ; seg(Q+1,Q+1-1i) ; seg(Q-1i,Q-1i+1) seg(0i,Q,col=8) pt(0.05+0.05i,'θ'); pt(Q/2,'r',cex=1.5) 948 名前：卵の名無しさん mailto:sage [2020/11/19(木) 13:17:26.24 ID:AV9OAvON.net] 会食が原因になっていることが多いと言っておいて、Gotoイートで会食を促進するという謎の政策を進める政府。 デフレ下に増税するという世襲低学歴の暗愚が最長政権の土人の国になったなぁ。 一方、文明国では 武漢からの飛行機は乗客降ろさず機内で検査。 海外からの入国は全て隔離検査。 電車バスはマスク無しは乗車拒否かつ検温を改札で実施。 店はマスク無しは入店拒否。 マスクはコンビニで無料配布かつスマホで一人一人の枚数管理するアプリを政府が開発。 学校は閉鎖しオンライン授業に切り替え。 949 名前：卵の名無しさん mailto:sage [2020/11/19(木) 15:46:43.48 ID:AV9OAvON.net] source('toolmini.R') plot(NULL,xlim=c(0,2),ylim=c(-1,1),asp=1,ann=F,axes=T) O=0i;A=1+0i;B=1i;C=1+1i pt(O,'O'); pt(A,'A') ; pt(B,'B') ; pt(C,'C') seg(O,A);seg(0,B);seg(C,A);seg(C,B) ; seg(O,C,lty=3) r=1.070437 th=-0.3647908 ; pt(0.15-0.025i,'θ') # Pの偏角 P=r*cos(th)+r*1i*sin(th) ; pt(P,'P(rcosθ,rsinθ)') ; seg(O,P,col=8) Cir(P,r,col=8) ci <- function(x) r*sin(th) + sqrt(r^2 - (x-r*cos(th))^2) Q=1+1i*ci(1) ; pt(Q,'Q') ; seg(P,Q,col=8) x=seq(0,1,0.01) lines(x,ci(x),type='h',col=4) pt(P/2,'r',cex=1.25) # area <- function(r,th){ # 円盤片の面積を算出 ci <- function(x) r*sin(th) + sqrt(r^2 - (x-r*cos(th))^2) integrate(ci,0,1)$value } area=Vectorize(area) # demo area(1.146384,-pi/8) 950 名前：卵の名無しさん mailto:sage [2020/11/19(木) 15:46:48.74 ID:AV9OAvON.net] th2r <- function(theta){ # θから円弧長を決定 r=uniroot(function(x) area(x,th=theta) - 1/2, c(1,2),tol=1e-24)$root #面積が1/2になるr ci <- function(x) r*sin(theta) + sqrt(r^2 - (x-r*cos(theta))^2) Q=1+1i*ci(1) P=r*cos(theta)+1i*r*sin(theta) r*bac(Q,P,0i) } th2r=Vectorize(th2r) # demo th2r(-pi/6) x=seq(-pi/7,-pi/9, 0.0005) plot(x,th2r(x),type='l',lwd=2,xlab=bquote(theta),ylab='円弧 OQ',cex.lab=1.25) (opt=optimize(th2r,c(-pi/7,-pi/9 951 名前：),tol=1e-24)) uniroot(function(x) area(x,th=opt$minimum) - 1/2, c(1,2),tol=1e-24)$root # r [] [ここ壊れてます] 952 名前：卵の名無しさん mailto:sage [2020/11/19(木) 15:48:01.54 ID:AV9OAvON.net] # 正方形の紙を頂点から切り始め、円周の一部（円弧）で面積を二等分するとき、 # 切り取り線の長さを出来るだけ短くするにはどうしたらよいか？ 953 名前：卵の名無しさん mailto:sage [2020/11/19(木) 17:41:10.39 ID:AV9OAvON.net] " ガチャの確率問題 とあるソーシャルゲームのガチャにはN種類のキャラがいて、全て等確率で出てくる。 ｎ連ガチャを回したとき、どのキャラでもいいので同じキャラがｍ回以上出てくる確率を求めよ。" " mean(replicate(1e7,any(table(sample(N,n,rep=T))>=m))) 954 名前：卵の名無しさん mailto:sage [2020/11/19(木) 18:51:57.01 ID:AV9OAvON.net] rm(list=ls()) options(digits=10) source('toolmini.R') plot(NULL,xlim=c(0,2),ylim=c(-1,1),asp=1,ann=F,axes=T) O=0i;A=1+0i;B=1i;C=1+1i pt(O,'O'); pt(A,'A') ; pt(B,'B') ; pt(C,'C') seg(O,A);seg(0,B);seg(C,A);seg(C,B) ; seg(O,C,lty=3) r=1.070437 th=-0.3647908 ; pt(0.15-0.025i,'θ') # Pの偏角 P=r*cos(th)+r*1i*sin(th) ; pt(P,'P(rcosθ,rsinθ)') ; seg(O,P,col=8) Cir(P,r,col=8) ci <- function(x) r*sin(th) + sqrt(r^2 - (x-r*cos(th))^2) Q=1+1i*ci(1) ; pt(Q,'Q') ; seg(P,Q,col=8) x=seq(0,1,0.005) lines(x,ci(x),type='h',col=4) pt(P/2,'r',cex=1.25) ; seg(O,C,lty=3) " integral r*sin(θ) + sqrt(r^2 - (x-r*cos(θ))^2) dx = " f <- function(x,r=1.070436,θ=-0.3647908) r*sin(θ) + sqrt(r^2 - (x-r*cos(θ))^2) curve(f(x),0,1,add=T) integrate(f,0,1)$value 955 名前：卵の名無しさん mailto:sage [2020/11/19(木) 18:52:02.15 ID:AV9OAvON.net] Fn <- function(x,r=1.070436,θ=-0.3647908){ (1/2)*(x-r*cos(θ))*sqrt(r^2 - (x-r*cos(θ))^2) - (1/2)*r^2*atan((r*cos(θ)-x)/sqrt(r^2-(x-r*cos(θ)^2))) + r*x*sin(θ) } Fn(1)-Fn(0) Area <- function(r,θ){ (1/2)*(1-r*cos(θ))*sqrt(r^2-(1-r*cos(θ))^2) - (1/2)*atan((r*cos(θ)-1)/sqrt(r^2+r*cos(θ)^2-1)) + r*sin(θ) + (1/2)*r*cos(θ)*sqrt(r^2-r^2*cos(θ)^2) + (1/2)*r^2*atan(r*cos(θ)/sqrt(r^2+r*cos(θ)^2)) } Area(r=1.070436,θ=-0.3647908) 956 名前：卵の名無しさん mailto:sage [2020/11/21(土) 09:04:39.17 ID:YdGgHDOx.net] 田村議員の質疑には、民主主義と法治国家を死守するという熱意に溢れている。 これに対する閣僚、官僚の答弁は法治国家を破壊しかねない、デタラメなものばかり。 https://www.youtube.com/watch?v=E9kPDgbczXk https://www.youtube.com/watch?v=bQOwuAiKseU 957 名前：卵の名無しさん mailto:sage [2020/11/21(土) 14:05:50.26 ID:vNIkaki5.net] 「感染拡大は、Go to が原因とする科学的エビデンスはない」と自民党の連中が言うけどさ ゴールデンウィーク前は「外出するな、旅行するな、帰省もするな」の大合唱だったんだよ。「人が移動したらウィルスも移動する」からね それがいつの間にか、マスクして感染対策すりゃOK！ってなってるの何で？ 958 名前：卵の名無しさん mailto:sage [2020/11/21(土) 18:12:40.62 ID:yHIGlwPH.net] LET'S 振込め詐欺見張隊 新117 L-FSM-N117 www.lets-co.jp/lets/security/furikomenew117.html 959 名前：卵の名無しさん mailto:sage [2020/11/22(日) 11:00:36.59 ID:m/vS68nE.net] # https://www3.nhk.or.jp/news/special/coronavirus/data-all/ # 2019年の交通事故による死亡者の数は3215人、https://www.nippon.com/ja/japan-data/h00625/ rm(list=ls()) x=2596 # new infection r=1981 # culmutive death toll n=130951 # culmutive infection r/n # morality rate x*r/n # estimated deaths in new infecton ci=binom.test(r,n)$conf data.frame(lower=ci[1]*x,upper=ci[2]*x) d=365.2425 x*r/n*d # annual estimated as.vector(x*ci*d) x*r/n*d/3215 # ratio to traffic accident casul as.numeric(x*ci*d/3215) 960 名前：卵の名無しさん mailto:sage [2020/11/22(日) 11:09:15.63 ID:m/vS68nE.net] # https://www3.nhk.or.jp/news/special/coronavirus/data-all/ # 2019年の交通事故による死亡者の数は3215人、https://www.nippon.com/ja/japan-data/h00625/ rm(list=ls()) x=2596 # new infection r=1981 # culmutive death toll n=130951 # culmutive infection r/n # morality rate x*r/n # estimated deaths in new infection ci=binom.test(r,n)$conf # 95% confidence interval(C.I.) data.frame(lower=ci[1]*x,upper=ci[2]*x) d=365.2425 # days in one-year x*r/n*d # annual estimated as.vector(x*ci*d) # 95%C.I. x*r/n*d/3215 # ratio to traffic accident death casualities as.numeric(x*ci*d/3215) # 95%C.I. 961 名前：卵の名無しさん mailto:sage [2020/11/22(日) 13:07:16.02 ID:m/vS68nE.net] # https://www3.nhk.or.jp/news/special/coronavirus/data-all/ # 2019年の交通事故による死亡者の数は3215人、https://www.nippon.com/ja/japan-data/h00625/ rm(list=ls()) x=2596 # new infection r=1981 # culmutive death toll n=130951 # culmutive infection r/n # morality rate x*r/n # estimated deaths in new infection ci=binom.test(r,n)$conf # 95% confidence interval(C.I.) data.frame(lower=ci[1]*x,upper=ci[2]*x) d=365.2425 # days in one-year x*r/n*d # annual estimated as.vector(x*ci*d) # 95%C.I. ta=3215 # annual traffic accident fatality x*r/n*d/ta # ratio to traffic accident fatality as.vector(x*ci*d/ta) # 95%C.I. 962 名前：卵の名無しさん mailto:sage [2020/11/22(日) 17:43:25.52 ID:m/vS68nE.net] # https://www3.nhk.or.jp/news/special/coronavirus/data-all/ # 2019年の交通事故による死亡者の数は3215人、https://www.nippon.com/ja/japan-data/h00625/ rm(list=ls()) x=2596 # new infection r=1981 # culmulative death toll n=130951 # culmulative infection r/n # morality rate x*r/n # estimated deaths in new infection ci=binom.test(r,n)$conf # 95% confidence interval(C.I.) data.frame(lower=ci[1]*x,upper=ci[2]*x) d=365.2425 # days in one-year x*r/n*d # annual estimated as.vector(x*ci*d) # 95%C.I. ta=3215 # annual traffic accident fatality x*r/n*d/ta # ratio to traffic accident fatality as.vector(x*ci*d/ta) # 95%C.I. 963 名前：卵の名無しさん mailto:sage [2020/11/23(月) 12:41:47.57 ID:jv+E87O2.net] 小選挙区制で（与党）議員どおしのcheck & balanceが消滅し、 内閣人事局で政治家と官僚のcheck & balanceが消滅。 制度化したのが民主党、制度を悪用したのが自民党。 その結果が 　デフレ下に増税、 　疫病下に旅行推進 という歴史に残る愚策。 964 名前：卵の名無しさん mailto:sage [2020/11/23(月) 23:23:57.53 ID:jv+E87O2.net] 小選挙区制で（与党）議員どおしのcheck & balanceが消滅し、 内閣人事局で政治家と官僚のcheck & balanceが消滅。 制度化に動いたのが民主党を作った政治家達、制度化して悪用したのが自民党。 その結果が 　デフレ下に増税、 　疫病下に旅行推進 という歴史に残る愚策。 965 名前：卵の名無しさん mailto:sage [2020/11/25(水) 16:52:40.72 ID:1mPkdWnl.net] 財政破綻論を唱える連中には、この一言で全ての片が付く。 「日銀のFTF買い入れの財源教えて」 966 名前：卵の名無しさん mailto:sage [2020/11/26(木) 08:32:39.74 ID:/Ibshp5i.net] 財政破綻論を唱える連中には、この一言で全ての片が付く。 「日銀のETF買い入れの財源教えて」 967 名前：卵の名無しさん [2020/11/26(木) 19:00:26.74 ID:Ok784nx+.net] 医師版になぜ統計？ 俺、統計できるんだぜ お前らと違って あたまいいだろ ホルホル、つうこと？ 968 名前：卵の名無しさん mailto:sage [2020/11/26(木) 22:10:36.52 ID:LFiRhrQe.net] IMO <- function(x='IMO'){ y=unlist(strsplit(x,'')) n=length(y) z=numeric(n) for(i in 1:n) z[i]=which(y[i]==LETTERS) w=z-1 sum(w*26^((n-1):0))+1 } IMO('KINTAMA') IMO('KKK') 969 名前：卵の名無しさん mailto:sage [2020/11/26(木) 22:36:55.28 ID:LFiRhrQe.net] pm=expand.grid(replicate(3,LETTERS,simplify=F)) pm=t(apply(pm,1,rev)) f <- function(x){ x[1]=='I' & x[2]=='M' & x[3]=='O' } which(apply(pm,1,f)) pm[5735,] revLET <- function(x='IMO'){ y=unlist(strsplit(x,'')) n=length(y) z=numeric(n) for(i in 1:n) z[i]=which(y[i]==LETTERS) z } revLET('IMO') (9-1)*26^2+(13-1)*26+15 IMO <- function(x='IMO'){ y=unlist(strsplit(x,'')) n=length(y) z=numeric(n) for(i in 1:n) z[i]=which(y[i]==LETTERS) w=z-1 sum(w*26^((n-1):0))+1 } 970 名前：卵の名無しさん mailto:sage [2020/11/26(木) 22:37:08.75 ID:LFiRhrQe.net] IMO('KINTAMA') IMO('SEX') IMO('YMO') source('toolmini.R') pm[12296,] revLET("SEX") dec2nw(12296,26,0) n2IMO <- function(x){ y=dec2nw(x,26,0)+1 n=length(y) y[n]=y[n]-1 print(paste(LETTERS[y],collapse = ''),quote=F) } n2IMO(2020) 971 名前：卵の名無しさん mailto:sage [2020/11/27(金) 00:04:49.17 ID:MEMGUU7d.net] １レスいくらだか知らねーけど、そんなバイトしなきゃならんほど困窮させた張本人擁護とかよくやるわ 972 名前：卵の名無しさん mailto:sage [2020/11/27(金) 07:56:33.46 ID:MEMGUU7d.net] IMO <- function(x){ y=unlist(strsplit(x,'')) n=length(y) z=numeric(n) for(i in 1:n) z[i]=which(y[i]==LETTERS) w=z-1 sum(w*26^((n-1):0))+1 } > IMO('CZR') [1] 2020 n2IMO <- function(x){ r=x%%26 q=x%/%26 while(q > 0){ r=append(q%%26,r) q=q%/%26 } y=r+1 n=length(y) y[n]=y[n]-1 print(paste(LETTERS[y],collapse = ''),quote=F) } > n2IMO(2020) [1] CZR 973 名前：卵の名無しさん mailto:sage [2020/11/27(金) 10:40:53.44 ID:0gxQ/Sbs.net] 中国の数字は信用できないと言われるが公文書改竄に賄賂五輪誘致GDPやらコロナやら数字誤魔化してる日本は言える立場に無い 974 名前：卵の名無しさん [2020/11/27(金) 11:10:02.25 ID:DvhL2xEx.net] ナースステーションの端で泣いている男性研修医(38)を見かけた時の適切な対応を１つ選べ。 a. 蹴り飛ばす b. 抗うつ薬を処方する c. 退職を勧める d. 「サボっている研修医が居る」と院長へ報告する これ本当に国試に出された問題なんだって。 ふざけるのもいい加減にしろと言いたいわ。 出題者のレベルが下がってるのは、やはり 学生の質が低下してるのと関係があるんだろうな。 しかし精神的に幼稚な出題者だよなｗ 出題者が一番、医師としての適性が欠けてるじゃんｗ 975 名前：卵の名無しさん mailto:sage [2020/11/27(金) 12:54:33.84 ID:0gxQ/Sbs.net] >>897 976 名前：英語と統計は現代の読み書きソロバン [] [ここ壊れてます] 977 名前：卵の名無しさん [2020/11/27(金) 13:27:18.61 ID:LgYm5vL7.net] そりゃ全員裏口かブーストかどっちかだもんな 978 名前：卵の名無しさん mailto:sage [2020/11/27(金) 15:36:24.87 ID:MEMGUU7d.net] >>897 こういうのができないと金を毟り取られちゃう。 COVID19の潜伏期間の論文 https://www.nejm.org/doi/full/10.1056/NEJMoa2001316 結論は #--- incubation period --- # from Li et al NEJM 2020 # lognormal mean = 5.2 ln.par1 = 1.434065 ln.par2 = 0.6612 あるシリツ医が新型コロナ肺炎に罹患したとする。 行動調査によって発症前にキャバクラに行っており接客したキャバ嬢がシリツ医発症の２日後に発症していたことがわかった。 キャバ嬢はシリツ医から移されたと主張して１億円の賠償を求めている。 潜伏期間には幅がありキャバ嬢から移された可能性もあると主張してその確率を計算して賠償金を値切りたい。 いくら値切れるか計算せよ。 979 名前：卵の名無しさん mailto:sage [2020/11/27(金) 16:10:54.68 ID:MEMGUU7d.net] >>897 学部教員の声某 会議にて某教授曰く「統計できない人は進級させなくていいですから」 別の某教授からの怒りの電話：「基本的なことがわかっていない学生・大学院生が多い．教養部のうちにちゃんと教育してもらわないと！」 www.tmd.ac.jp/artsci/math/lec/tokunaga/statistics08_01.pdf 980 名前：卵の名無しさん mailto:sage [2020/11/27(金) 16:11:58.85 ID:MEMGUU7d.net] 最終的には多くの学生が「わかってしまえば簡単だった」という感想を漏らします． 皆さんの先輩M君（現M４）の名言：「自転車と同じで，いったん乗れるようになればあとは楽々ですよ」 981 名前：卵の名無しさん mailto:sage [2020/11/27(金) 19:45:15.34 ID:xmdXpYdQ.net] # https://www3.nhk.or.jp/news/special/coronavirus/data-all/ # 2019年の交通事故による死亡者の数は3215人、https://www.nippon.com/ja/japan-data/h00625/ rm(list=ls()) fn <- function( x=2523, # new infection r=2086, # culmulative death toll n=142811 # culmulative infection ){ cat('mortality rate = ',r/n,'\n') # morality rate cat('estimated deaths in newly infected = ',x*r/n,'\n') # estimated deaths in new infection ci=binom.test(r,n)$conf # 95% confidence interval(C.I.) print(data.frame(lower=ci[1]*x,upper=ci[2]*x)) d=365.2425 # days in one-year cat('annual deaths estimated = ',x*r/n*d, '\n') # annual estimated CI=as.vector(x*ci*d) # 95%C.I. print(data.frame(lower=CI[1],upper=CI[2])) ta=3215 # annual traffic accident fatality cat('ratio to traffic accident fatality = ',x*r/n*d/ta,'\n') # ratio to traffic accident fatality CI=as.vector(x*ci*d/ta) # 95%C.I. print(data.frame(lower=CI[1],upper=CI[2])) } fn() 982 名前：卵の名無しさん mailto:sage [2020/11/27(金) 19:56:10.93 ID:xmdXpYdQ.net] テレビのドキュメンタリーでやってたが 看護師が、これだけ働いて感染の危険も冒してるのにボーナス無しはあり得ないと病院側に訴え。 でも病院側は、利益が出なければボーナスなどとても出せない、経営が行き詰まってると。 実際私立病院ならそうだろう。 こういう医療現場こそコロナ補助で国が救済するべき。 なんで飲み屋や観光業など娯楽業種を助けるために 体張ってる医療機関を犠牲にするんだか。 983 名前：卵の名無しさん mailto:sage [2020/11/27(金) 21:46:30.17 ID:xmdXpYdQ.net] あるドキュメンタリーでやってたが 看護師が、これだけ働いて感染の危険も冒してるのにボーナス無しはあり得ないと病院側に訴え。 でも病院側は、利益が出なければボーナスなどとても出せない、経営が行き詰まってると。 実際、私立病院ならそうだろう。 こういう医療現場こそコロナ補助で国が救済するべきと思う。 医療機関で適正な賃金を払えない現実こそ、医療崩壊だと思う。 984 名前：卵の名無しさん mailto:sage [2020/11/28(土) 07:45:58.02 ID:Vq0tXK6x.net] 生活圏を超える人の移動は、確実に接触率を上げ、これが実効再生産数を挙げて、感染者数上昇を引き起こす。 現在、GOTOでウィルスを全国に拡散中。 しばらくは、旅行を楽しむような若者や元気な人が隠れ陽性となり、これらが顕在化して爆増だな。 985 名前：卵の名無しさん mailto:sage [2020/11/28(土) 08:38:33.70 ID:Vq0tXK6x.net] DR. FAUCI: You know, every single day we meet with the task force, and we take a look at what's going on. And you don't want to make a pronouncement that no one should ever go into a restaurant. I mean, I think that might be overkill right now, but everything is on the table. It may come to the situation where we strongly recommend. Right now, myself personally, I wouldn't go to a restaurant. I just wouldn't because I don't want to be in a crowded place. I have an important job to do. I don't want to be in a situation where I'm going to be all of a sudden self-isolating for 14 days. https://www.cbsnews.com/news/transcript-dr-anthony-fauci-discusses-coronavirus-on-face-the-nation-march-15-2020/ 986 名前：卵の名無しさん mailto:sage [2020/11/28(土) 09:18:40.36 ID:Vq0tXK6x.net] 安倍ぴょんまでも総合的・俯瞰的に粛清が始まってワクワクするのだが、朝敵安倍を撃ての詔が総理に下ったのかな？ 桜を見る会でなく、前夜祭に焦点をあてているのは菅が無傷で安倍を撃ち落とすためだろうな。 桜を見る会そのものを取り上げると、スガもプリズン　になってしまう。 987 名前：卵の名無しさん mailto:sage [2020/11/29(日) 12:10:27.26 ID:WqPJnRzs.net] バカはトリアージにより老人が追い出されると思って喜んでいたみたいだけど 現実には老人が優先されたというオチ まあ若者はコロナで死なないけど老人はすぐ死んでしまうから 老人優先は当然なんだよなあ 988 名前：卵の名無しさん mailto:sage [2020/11/29(日) 18:45:11.35 ID:WqPJnRzs.net] 感染者が入院希望しても入院できないのが医療崩壊だと思うのだが、現状は医療崩壊ではないのか？ 989 名前：卵の名無しさん mailto:sage [2020/11/29(日) 18:47:35.39 ID:WqPJnRzs.net] 感染者が入院希望しても入院できないのが医療崩壊だと思うのだが、現状は医療崩壊ではないのか？ 990 名前：卵の名無しさん mailto:sage [2020/11/29(日) 20:16:35.60 ID:WqPJnRzs.net] あるシリツ医大で裏口入学率を推定したい。 何の情報もないのでその確率を0から1までの>>488 一様分布とする。 調査により何人の裏口が続いたら裏口入学率の95%信頼区間の下限値が90%を超えるか？ そのときの裏口入学率の期待値はいくらか？ library(HDInterval) f <- function(x) hdi(qbeta,0.95,shape1=1+x,shape2=1) f=Vectorize(f) n=0:100 y=f(n) plot(n,(1+n)/(2+n),bty='l',ylab='Prob ',ylim=c(0,1),pch=19) segments(n,y['lower',],n,y['upper',]) # points(n,y['lower',]) # points(n,y['upper',]) abline(h=0.90,col=8) n0=uniroot(function(x,u0=0.90) f(x)[1] - u0, c(0,100))$root ; n0 (n0+1)/(n0+2) (ceiling(n0)+1) / (ceiling(n0)+2) 991 名前：卵の名無しさん mailto:sage [2020/11/29(日) 21:30:08.11 ID:WqPJnRzs.net] k=1e3 oz=0.6/0.4 game <- function(x){ OZ=P=numeric() OZ[x]=oz^(x-1) P[x]=OZ[x]/(1+OZ[x]) rbinom(k,1,P[x]) } sim <- function(a,b){ mean(game(a) > game(b)) } sim(7,3) 992 名前：卵の名無しさん mailto:sage [2020/11/30(月) 06:52:25.73 ID:QN29Tqgt.net] 海外からのコロナ流入を止めるのが遅すぎた　中国だけでなく欧州からの流入が続いたのが痛かった 検査の強化と 993 名前：隔離施設の手当が必要なのに、中国でそれが有効だったのを見ていたのに学ぶことなく折角クラスター対策で稼いだ時間を無駄にした GoToキャンペーンはコロナが終息してから開始する約束だったのに反故にして前倒し実施して傷口を広げた 科学的根拠を重視せず思い付きで無意味な施策を繰り出した（一斉休校、アベノマスク、GoToキャンペーン）それもお友達利権つきで　それらの施策は国会で審議もしていないから野党はどうすることもできないし、野党は桜を見る会だけ取り上げていたわけではない 少なくとも野党は「自民と同罪」ではない 君は結局自民のサポーターで、政府への批判が避けられないから野党も一緒だとレッテル貼りをしているのだろう [] [ここ壊れてます] 994 名前：卵の名無しさん mailto:sage [2020/11/30(月) 06:52:41.10 ID:QN29Tqgt.net] 同じ島国なのに、支那ウィルスの抑え込みに成功した台湾とニュージーランド。 失敗した日本とイギリス。 この同じ島国であるながら、違いがなぜ生じたかを考えればわかる。 アベがFラン過ぎた。 ジョンソンも支那ウィルスを当初、極東の風土病くらいにしか思わず、対策を油断した。 一方台湾もNZも観光も主要産業だったが、1月に危機感を感じトップダウンで鎖国した。 島国の利点を生かし、外国からの人の流れを遮断し鎖国することで乗り切った。 1月、2月は経済が縮小したが、その成果が表れて現在は支那ウィルスを抑え込んだ。 今は国内の経済活動も解禁した。 それに比べて日本はどうだい？ アベは1月に全く危機感も持たず、「春節の支那人観光客のインバウンドを歓迎します。」と呑気なビデオメッセージを垂れ流してた。 アベは「非常事態宣言で鎖国する法的根拠がないー」とほざいて言い逃れする。 しかしアベにやる気があれば、総理大臣の伝家の宝刀「超法規的措置」が取れた。 日本赤軍がよど号をハイジャックし、刑務所の極左メンバーの釈放と、北チョン国への亡命を要求した。 当時の総理は、刑法、入管法も停止して、法秩序を捻じ曲げて、国交のない北チョン国に赤軍メンバーを亡命させた。 「超法規的措置」とはこのように強力な権力である。 しかしFラン大卒のアベは、1月に日本の鎖国…外人の入国禁止措置を取らなかった。 馬鹿のアベに危機感と責任感と決断力がなかった。 ま、1月の国会を「桜を見る会がー！」で無駄に貴重な時間を潰した糞ボケ野党どもも、痔民党と同罪だけどな。 日本にはろくな政党がない。 日本には”真の”保守政党が必要だ。 　 995 名前：卵の名無しさん mailto:sage [2020/11/30(月) 22:24:00.68 ID:HI2s9ZOn.net] f <- function( wAB=80, # Aの対B勝数 lAB=20, # Aの対B負数 wAC=30, # Aの対C勝数 lAC=20, # Aの対C負数 k=1e6){ # 乱数発生数 b=rbinom(k,wAB+lAB,lAB/(wAB+lAB)) # Bの勝数乱数 ba=b/(wAB+lAB-b) ; hist(ba,br='sc',fr=F) # Bのオッズ乱数 c=rbinom(k,wAC+lAC,lAC/(wAC+lAC)) # Cの勝数乱数 ca=c/(wAC+lAC-c) ; hist(ca,br='sc',fr=F) # Cのオッズ乱数 bc=ba/ca # Bの対Cオッズ乱数 summary(bc) hist(bc,br='sc',fr=F) P 996 名前：bc=bc/(1+bc) hist(Pbc) summary(Pbc) pbc=Pbc[!is.na(Pbc)] print(summary(pbc)) plotPost(pbc,compVal = 0.5,showMode=T,showCurve = T,xlab='Prob[B beats C]') hdi(pbc) } [] [ここ壊れてます] 997 名前：卵の名無しさん mailto:sage [2020/11/30(月) 22:27:58.54 ID:HI2s9ZOn.net] library(gmp) f <- function(N){ a=factorialZ(N) b=as.character(a) c=as.numeric(unlist(strsplit(b,''))) d=rev(c) if(d[1]!=0) return(NA) else return(rle(d)$values[2]) } N=1:1000 y=sapply(N,f) plot(N,y,bty='n') z=data.frame(N,digit=y) head(z,10) tail(z,10) 998 名前：卵の名無しさん mailto:sage [2020/12/01(火) 08:54:10.45 ID:iy4usOEg.net] library(gmp) f <- function(N){ a=factorialZ(N) b=as.character(a) c=as.numeric(unlist(strsplit(b,''))) d=rev(c) if(d[1]!=0) return(NA) else return(rle(d)$values[2]) } N=1:10004 y=sapply(N,f) plot(N,y,bty='n') z=data.frame(N,digit=y) head(z,10) tail(z,10) a=z[-1:-4,] # non-zero digit b=vector('list',length=9) for(i in c(2,4,6,8)){ b[[i]]=a[a$digit==i,]$N } b c=sapply(c(2,4,6,8), function(i) length(b[[i]])) data.frame(digit=c(2,4,6,8),個数=c) 999 名前：卵の名無しさん mailto:sage [2020/12/01(火) 10:28:29.05 ID:VGK+tvyz.net] # 差X-Yのpdfの公式 # ∫[-∞,∞] pdf1(x+y)*pdf2(y) dy # x ~ pdf1,y ~ pdf2から x-y のpdfを返す pdf_minus <- function(pdf1,pdf2){ f <- function(x,y) pdf1(x+y)*pdf2(y) f=Vectorize(f,vectorize.args = 'y') pdf <- function(x) integrate(function(y) f(x,y),-Inf,Inf,rel.tol=1e-14)$value pdf=Vectorize(pdf) invisible(pdf) } # ex. COVID19 incubation time mu = 1.434065 ; sg = 0.6612 x=rlnorm(1e6,mu,sg) ; y=rlnorm(1e6,mu,sg) hist(x-y,xlim=c(-30,30),freq=F,breaks=250,lty=1,col='skyblue',main='') pdf_minus(function(x)dlnorm(x,mu,sg), function(x)dlnorm(x,mu,sg)) -> pdf curve(pdf(x),add=T) var=integrate(function(x) x^2*pdf(x),-Inf,Inf)$value sqrt(var) integrate(pdf,2,Inf,rel.tol=1e-12) mean(x-y>2) library(polspline) fit=logspline(x-y) curve(dlogspline(x,fit),add=T,col=4) 1-plogspline(2,fit) 1000 名前：卵の名無しさん mailto:sage [2020/12/01(火) 12:55:50.90 ID:iy4usOEg.net] 日本の誇り ★亀井静香 「簡単に言いますと・・・ 小泉・竹中改革と称する路線の逆をやれば日本の未来が開かれる、 このように私は思います。」 https://www.youtube.com/watch?v=DmiCQKBM2So&t=350s

---

---

[ 続きを読む ] / [ 携帯版 ]

read.cgi ver5.27 [feat.BBS2 +1.6] / e.0.2 (02/09/03) / eucaly.net products.
担当:undef