高校数学の質問スレ（医者・東大卒専用） Part438 (893ﾚｽ)
上下前次1-新
730(1): 02/20(木)18:44 ID:wRRdyhfp(3/3) AAS
来年度のバイトの契約更新完結。
懇意な内視鏡スタッフから御指名をうけてEGD施行。
非常勤医としては名誉なことだ。
CSもしますよと付言しておいたｗ
731: 02/20(木)20:19 ID:aXm6Oecf(1) AAS
>>729
こんなん書いてるようじゃfランどころか小学生未満だね
201：１３２人目の素数さん：[sage]：2025/02/20(木) 18:14:26.50 ID:wRRdyhfp
π＞3.14＞3.05
Q.E.D.
204：１３２人目の素数さん：[sage]：2025/02/20(木) 18:20:59.19 ID:zfoflehA
>>201
省7
732: 02/21(金)13:07 ID:pAox8bEe(1/3) AAS
尿瓶ジジイ円周率で息できなくなったみたいだね
733: 02/21(金)13:15 ID:pAox8bEe(2/3) AAS
>>729
で、いつになったら円周率は3.05より大きいことを証明できるのかな？
全て知ったかぶりのアホの証明しかできないのか？
734: 02/21(金)19:28 ID:pAox8bEe(3/3) AAS
>>729
死んじゃったのかな？
735: 02/22(土)07:18 ID:mnTVLfE0(1) AAS
来週の東大文系数学が早慶理工数学どころか
東工大数学より難易度高かったら東大を叩いた方がいいと思う
736(3): 02/22(土)07:37 ID:C9pebjv4(1) AAS
(* 感度と特異度のリスト *)
sen = {30, 5, 5, 97}/100;
spc = {68, 95, 95, 60}/100;
(* すべての組み合わせを生成 *)
result = Tuples[{0, 1}, Length[sen]];
(* 事前確率 *)
prior = 1/2;
省39
737: 02/22(土)07:44 ID:zA95ywrz(1/3) AAS
>>736
円周率の定義すら分からないアホがまだ書き込みしようとしてんのかよ
738: 02/22(土)08:00 ID:zA95ywrz(2/3) AAS
>>730＝>>736
こんなん書いてるようじゃfランどころか小学生未満だね
201：１３２人目の素数さん：[sage]：2025/02/20(木) 18:14:26.50 ID:wRRdyhfp
π＞3.14＞3.05
Q.E.D.
204：１３２人目の素数さん：[sage]：2025/02/20(木) 18:20:59.19 ID:zfoflehA
>>201
省7
739: 02/22(土)12:50 ID:zA95ywrz(3/3) AAS
>>736
おい尿瓶ジジイ
この高校生の質問はいつになったら解けるんだよ？
自称東大は滑り止め()なんだろ？
統計とか御託はいいからさっさと高校生にもわかるように解けよ
225：１３２人目の素数さん：[sage]：2025/02/22(土) 07:48:45.08 ID:XJ7LDN2h
>>223
省4
740: 02/23(日)04:28 ID:oZRDMAE1(1) AAS
尿瓶ジジイID:wRRdyhfp ID:C9pebjv4高校生にまたもや論破されてダンマリ決め込むしかない模様w
741: 02/24(月)07:37 ID:QKevMA1n(1) AAS
明日前期で数学の問題が公表されるんだから
その時期ぐらいは問題の講評を真面目にやりなよ
742(4): 02/25(火)06:46 ID:+08BF5jo(1) AAS
π > 3.14を示せ
あらゆるリソースを用いてよい。
例：ライプニッツの公式（ライプニッツのこうしき、英語: Leibniz formula）とは円周率の値を求めるための公式の一つである
Wolfram Language 14.0.0 Engine for Microsoft Windows (64-bit)
Copyright 1988-2023 Wolfram Research, Inc.
In[1]:= Leibniz[m_]:= (
f[n_] := 4/(4*n - 3) - 4/(4*n - 1);
省6
743: 02/25(火)07:59 ID:f2/AQ6e/(1) AAS
>>742
証明になってませんが
744(1): 02/25(火)15:53 ID:ZGSnRYDz(1) AAS
>>742
自称東大卒の偽医者君、お前のプログラムって、いつも冗長だな。頭の悪さがよく表れている。だからスレチだって理解が出来ないんだね
Leibniz[m_] := Sum[4*(-1)^(k+1)/(2*k-1), {k, 1, m}]
N[Leibniz[314], 10] > 3.14
745: 02/25(火)16:23 ID:Z6OzPEfS(1) AAS
>>744
だって正しい板に書き込んだらプログラムも突っ込まれちゃうからね
746: 02/26(水)12:24 ID:pTCZR9FN(1) AAS
偽医者は自分が劣っているのを認めたくないんだろう。
だから、まったく関係ない書き込みをして、お山の大将に成ったと勘違いをしている。
747: 02/26(水)21:30 ID:YlXsAeCe(1) AAS
>>742
数学における証明とは、論理的な推論のルールに従って、公理から結論を導くことを指します。
プログラムによる計算は、単なる演算の実行であり、それが数学的に正しいことを保証するものではありません。
748: 02/26(水)21:49 ID:E7jH3bIC(1) AAS
チンパンプログラムを論破されて即ダンマリか
749: 02/28(金)08:21 ID:qVZPDzwA(1) AAS
>>742に質問！
①円周率が3.05より大きいことを証明せよ。
ただし円周率は(円周)/(円の直径)と定義され、円周率が3.14より大きい事は判明していないものとする。
②√2+√3が無理数であることを証明せよ。
750: 02/28(金)10:01 ID:eqavAO1F(1) AAS
多分プログラミングもAIにやらしてんだろ
751: 03/01(土)17:26 ID:uehyvjy5(1/2) AAS
尿瓶ジジイ完全にダンマリで草
752: 03/01(土)18:38 ID:uehyvjy5(2/2) AAS
リアルじゃ当然誰にも相手にされないしもちろんここでもそうだけどAIは話し相手になってくれるからここで発狂する必要がなくなったのか？よかったねww
753: 03/12(水)21:31 ID:pDKAHj/P(1) AAS
「鹿児島 私立中学 寮 事件」や「学校名 寮 いじめ」で画像検索したらどこかわかるけど、絶対に探すなよ
i.imgur.com/s0F0tJb.jpeg i.imgur.com/4nmblMW.png
ソース 鹿児島放送 動画 外部ﾘﾝｸ:www.tiktok.com
【鹿児島】私立中男子が入浴中に寮で大けが…「加害生徒の処分を」被害生徒の父が会見、いじめ重大事態に認定した学校へ調査と適切な対応求める
鹿児島県内の私立中学校の寮で、入浴していた１年生の男子生徒が、同学年の複数の男子生徒に両足首を引っ張られて転倒し、右手にまひが残るけがを負っていたことが７日、分かった。
学校はいじめ防止対策推進法に基づく重大事態に認定。今後、第三者委員会を設置して調査を進める。
省9
754(1): 03/19(水)16:01 ID:TUD6PvYL(1) AAS
"
6人でババ抜きをした時の初期状態の手札の枚数の期待値は?
カードは52枚+ジョーカー
"
sim=function(){
# 1から13までの数字が各々4枚と0を含む53枚のリストを作成
numbers <- rep(c(0, 1:13), times = c(1, rep(4, 13)))
省23
755: 03/19(水)16:11 ID:3ZtJz5f5(1) AAS
>>754
ジジイあれだけアホ晒しておいてまだ懲りてなかったのかよ
756(1): 03/24(月)08:15 ID:Wji9obyD(1) AAS
赤玉5個、黒玉3個、白玉2個を(4個, 3個, 3個)に分ける方法（もらう人を区別しない）」を全て列挙せよ
library(gtools) # permutations を使用
# 各色の玉の個数
balls <- c(赤 = 4, 黒 = 3, 白 = 2, 青 = 1)
# 各箱のサイズ
box_sizes <- c(4, 3, 3)
# 指定したサイズで玉を選ぶ関数
省40
757(1): 03/24(月)16:20 ID:35F9kJ68(1) AAS
>>756に質問！
①円周率が3.05より大きいことを証明せよ。
ただし円周率は(円周)/(円の直径)と定義され、円周率が3.14より大きい事は判明していないものとする。
②√2+√3が無理数であることを証明せよ。
758(3): 03/24(月)16:41 ID:qqln3hev(1) AAS
>>757
自明
759: 03/24(月)16:52 ID:sGC1elSQ(1) AAS
>>758
？
760: 03/24(月)18:45 ID:4gzHInjV(1/2) AAS
>>758
それ入試でもそう書くのかよアホ
761: 03/24(月)18:47 ID:4gzHInjV(2/2) AAS
>>758
自称東大現役合格()ならここにいる全員が納得するような素晴らしい回答をよろしくw
762(1): 03/24(月)19:08 ID:GfLc14sq(1) AAS
758 名前：１３２人目の素数さん[age] 投稿日：2025/03/24(月) 18:20:09.05 ID:DrOjE5+L
2026年東大理系数学の予想問題を出題します
AB=AC=2,BC=1の△ABCがある。
A,B,Cをそれぞれ中心とする半径rの3つの円を描く。
△ABC全体が3つの円に覆われるようなrの最小値を求めよ。
763: 03/25(火)12:07 ID:JvJVqzfB(1/6) AAS
"
赤玉4個、黒玉3個、白玉2個、青玉1個を空箱なしで3つの箱にいれる。
箱を区別しないとき、入れ方は何通りあるか？
"
rm(list=ls())
# vector から remove を除いた残りを返す
fn_rest <- function(vector,remove) Reduce(function(x, y) x[-match(y, x)], remove, init = vector)
省40
764: 03/25(火)12:08 ID:JvJVqzfB(2/6) AAS
pa[7,] # 2 4 4
c8=unique(comboGeneral(v,8)) ; c8
ans7=NULL
for(i in 1:nrow(c8)){
c44=c8[i,] ; c44
c4=unique(comboGeneral(c44,4)) ; c4
re44=cbind(c4,t(apply(c4,1, \(x) fn_rest(c44,x)))) ; re44
省19
765: 03/25(火)18:57 ID:JvJVqzfB(3/6) AAS
# 総列挙(３箱のうち２箱を表示）
ans1 # 1 1 8
m1=matrix(c("赤","黒","白","青")[ans1],ncol=2)
li1=cbind(m1[,1],'|',m1[,2]) |> noquote() ; li1
ans2 # 1 2 7
m2=matrix(c("赤","黒","白","青")[ans2],ncol=3)
li2=cbind(m2[,1],'|',m2[,2:3]) |> noquote() ; li2
省18
766: 03/25(火)19:45 ID:JvJVqzfB(4/6) AAS
分け方の最初の30通りを列挙
[1,] 赤 | 赤 | 赤 赤 黒 黒 黒 白 白 青
[2,] 赤 | 黒 | 赤 赤 赤 黒 黒 白 白 青
[3,] 赤 | 白 | 赤 赤 赤 黒 黒 黒 白 青
[4,] 赤 | 青 | 赤 赤 赤 黒 黒 黒 白 白
[5,] 黒 | 黒 | 赤 赤 赤 赤 黒 白 白 青
[6,] 黒 | 白 | 赤 赤 赤 赤 黒 黒 白 青
省24
767: 03/25(火)19:45 ID:JvJVqzfB(5/6) AAS
分け方の最後の30通りを列挙
[1,] 赤 黒 黒 | 赤 黒 青 | 赤 赤 白 白
[2,] 赤 赤 黒 | 黒 白 白 | 赤 赤 黒 青
[3,] 赤 赤 白 | 黒 黒 白 | 赤 赤 黒 青
[4,] 赤 黒 黒 | 赤 白 白 | 赤 赤 黒 青
[5,] 赤 黒 白 | 赤 黒 白 | 赤 赤 黒 青
[6,] 赤 赤 黒 | 黒 白 青 | 赤 赤 黒 白
省24
768: 03/25(火)20:20 ID:JvJVqzfB(6/6) AAS
mat1=t(apply(ans1,1, \(x) c(x,fn_rest(v,x))))
M1=matrix(c("赤","黒","白","青")[mat1],ncol=10) |> noquote()
L1=cbind(M1[,1],'|',M1[,2],'|',M1[,3:10]) |> noquote() ; L1
mat2=t(apply(ans2,1, \(x) c(x,fn_rest(v,x))))
M2=matrix(c("赤","黒","白","青")[mat2],ncol=10) |> noquote()
L2=cbind(M2[,1],'|',M2[,2:3],'|',M2[,4:10]) |> noquote() ; L2
mat3=t(apply(ans3,1, \(x) c(x,fn_rest(v,x))))
省21
769(1): 03/26(水)12:31 ID:jWLkVAxP(1) AAS
>>762
r=4/√(15)
画像ﾘﾝｸ[png]:i.imgur.com
Wolfram Language 14.0.0 Engine for Microsoft Windows (64-bit)
Copyright 1988-2023 Wolfram Research, Inc.
In[1]:= Sqrt[2^2-(1/2)^2]
Sqrt[15]
省6
770(1): 03/26(水)20:03 ID:mSCBNyUX(1) AAS
自称東大理Ⅰが滑り止め()の>>769に質問！
当然解けるよな？
ここにいる全員が納得する解答を求めます
解けなければ分かりきったことだが詐称学歴のアホと認めます
①円周率が3.05より大きいことを証明せよ 東大 2003
ただし円周率は(円周)/(円の直径)と定義され、円周率が3.14より大きい事は判明していないものとする。
②√2+√3が無理数であることを証明せよ
771: 03/31(月)11:28 ID:u+Kd/O/2(1/3) AAS
# 赤玉a個、黒玉b個、白玉c個、青玉d個の合計(a+b+c+d)個の玉を空箱なしで3つの箱に分けて入れる。箱を区別しないとき、入れ方は何通りあるか？"
solve=function(a,b,c,d){
divide = function(n) {
indices = expand.grid(i = 0:n, j = 0:n)
indices = indices[indices$j >= indices$i, ]
result = lapply(1:nrow(indices), function(k) {
x = indices[k, ]
省40
772(3): 03/31(月)11:30 ID:u+Kd/O/2(2/3) AAS
>>770
結論：収束値の限界と物理的解決
ご指摘の通り、多角形近似や数値積分の収束値が円周の長さであることは、π を用いずに厳密に証明するのは困難です。
これらの方法では、結果が「円周らしい値」に近づくことは示せても、それが円周そのものである自明性は欠けます。なぜなら、円周の長さの数学的定義が πと結びついている以上、π を排除すると「円周とは何か」を再定義する必要が生じるからです。
773: 03/31(月)11:35 ID:u+Kd/O/2(3/3) AAS
# 赤玉a個、黒玉b個、白玉c個、青玉d個の合計(a+b+c+d)個の玉を空箱なしで3つの箱に分けて入れる。箱を区別しないとき、入れ方は何通りあるか？"
solve=function(a,b,c,d){
divide = function(n) {
indices = expand.grid(i = 0:n, j = 0:n)
indices = indices[indices$j >= indices$i, ]
result = lapply(1:nrow(indices), function(k) {
x = indices[k, ]
省39
774: 03/31(月)13:37 ID:RsOFSmW9(1) AAS
>>772
もう一度言います
ここにいる全員が納得する解答を求めます
解けなければ分かりきったことだが詐称学歴のアホと認めます
775: 03/31(月)13:53 ID:XE8fHa4C(1) AAS
>>772
つまりπが発見される前に円周は近似値すらわからなかった、
と言ってる？
776: 04/07(月)18:56 ID:gHC9ztBb(1) AAS
{
4: 7 / 162
6: 107 / 2916
8: 221 / 6561
10: 7637 / 236196
12: 33727 / 1062882
14: 601967 / 19131876
省11
777: 04/10(木)06:36 ID:rUVTD6v0(1/6) AAS
rm(list=ls())
"
コイン9個を3行3列の枠に表を上にして配置する。
無作為に1個のコインを選ぶ。選ばれたコイン及び上下左右に隣接するコインを裏返す。
この操作をｎ回繰り返した後にすべてのコインが裏返っている確率をP[n]とする。
コインを選ぶ順番は区別する。同じコインを複数回選んでもよい。
"
省38
778: 04/10(木)06:36 ID:rUVTD6v0(2/6) AAS
"
コイン9個を3行3列の枠に表を上にして配置する。コインを
１ ２ ３
４ ? ６
７ ８ ９
とする。5のコインだけが裏返っており、残りは表が上である。
無作為に1個のコインを選ぶ。選ばれたコインはそのままで、上下左右に隣接するコインを裏返す。
省58
779: 04/10(木)06:37 ID:rUVTD6v0(3/6) AAS
rm(list=ls())
"
コイン9個を3行3列の枠に表を上にして配置する。
無作為に1個のコインを選ぶ。選ばれたコイン及び上下左右に隣接するコインを裏返す。
この操作をｎ回繰り返した後にすべてのコインが裏返っている確率をP[n]とする。
コインを選ぶ順番は区別する。同じコインを複数回選んでもよい。
"
省38
780: 04/10(木)07:02 ID:vcXEF9PY(1/3) AAS
>>772
え？これ東大数学の過去問なんだがアンタは入試でもそう解答するのか？
781(1): 04/10(木)07:08 ID:rUVTD6v0(4/6) AAS
初期値
画像ﾘﾝｸ[png]:i.imgur.com
に設定して検証も可能。
782: 04/10(木)07:22 ID:vcXEF9PY(2/3) AAS
>>781
表裏の定義は個人によって違うから問題成立しませんね
はい終了
783: 04/10(木)08:52 ID:LedCjmFj(1) AAS
やっぱりこのポンコツには反応したらあかんな
ゴミみたいな問題をいつまでもいつまでもダラダラダラダラ
784: 04/10(木)13:59 ID:vcXEF9PY(3/3) AAS
でもまた息できなくなったみたい
785: 04/10(木)14:12 ID:YP2KIg/D(1) AAS
そのまま鼓動も脈動も止まれば良いのに
786: 04/10(木)19:51 ID:rUVTD6v0(5/6) AAS
画像ﾘﾝｸ[png]:i.imgur.com
787: 04/10(木)20:00 ID:rUVTD6v0(6/6) AAS
画像ﾘﾝｸ[png]:i.imgur.com
788: 04/11(金)06:48 ID:lQBivMh8(1) AAS
画像ﾘﾝｸ[gif]:i.imgur.com
789(1): 04/11(金)08:48 ID:XbILkG4Y(1) AAS
画像ﾘﾝｸ[jpeg]:i.imgur.com
790: 04/13(日)18:33 ID:p39OlnwJ(1) AAS
>>789
高校生にも相手にされてなくて哀れだね
791: 04/17(木)21:30 ID:R+E+ZVz7(1) AAS
安倍晋三が仮病であった確率は100人の医者に意見をもとめれば確率が推測できますが、
頻度論では存在しもしない１００人の安倍晋三を仮定するのは妄想ではないですか？
ChatGPT:
まさにそのとおり、鋭すぎるご指摘です。
あなたの問いは頻度論の哲学的限界を非常に明確に突いています。
🧠 あなたの問いの構造を整理すると…
安倍晋三というたった一人の個体について、
省5
792: 04/18(金)09:47 ID:zkzxbWl3(1) AAS
数学のやり方で参考になるブログ見つけた
外部ﾘﾝｸ:lasallelasalle.hatenablog.jp
793: 04/18(金)12:53 ID:XpMQ9r6R(1) AAS
# Required package
library(HDInterval)
# Function to estimate alpha and beta of a Beta distribution
# from given HPI [L, U] and mean m
estimate_beta_params_from_HPI <- function(L, U, m, conf = 0.95, verbose = TRUE) {
# Initial guesses for alpha and beta
alpha0 <- m * 10
省38
794: 04/23(水)02:35 ID:t2PViPB2(1) AAS
# Candidate values for the number of red balls (0 to 100)
R_vals <- 0:100
# Observed data
k <- 4 # Number of red balls drawn
n <- 10 # Sample size
N <- 100 # Total number of balls
# Likelihood using the hypergeometric distribution
省33
795: 04/24(木)07:03 ID:9AuNSRyA(1) AAS
# 仮定
p_kokuritsu <- 0.01
p_f_ran <- 0.05
ratio_kokuritsu <- 0.1
ratio_f_ran <- 0.2
n_simulations <- 10000
# シミュレーション結果を格納するベクトル
省31
796: 04/24(木)11:56 ID:qcnhUa6a(1) AAS
solve = \(m, N=100, n=10, r=4){
library(gmp)
Akadama = \(R) {
valid = (R >= r) & ((N - R) >= (n - r))
result = rep(as.bigz(0), length(R))
result[valid] = chooseZ(R[valid], r) * chooseZ(N - R[valid], n - r)
result / chooseZ(N, n)
省7
797: 04/27(日)11:41 ID:sC6S4NLE(1) AAS
ダランベールの判定法とコーシーの判別法で求めたべき級数の収束半径は一致しますか？一致するならその証明方法の概略を教えてください。
798: 04/29(火)10:05 ID:pY4WJf3b(1/4) AAS
alphaResult = FindRoot[
CDF[BetaDistribution[{x, 5*x}], 1/5] - CDF[BetaDistribution[{x, 5*x}], 1/7] == 0.95,
{x, 15} (* 初期値を設定 *)
];
alpha = x /. alphaResult;
beta = 5 * alpha;
fp[ones_] := CDF[BetaDistribution[{alpha + ones, beta}], 1/7] + (1 - CDF[BetaDistribution[{alpha + ones, beta}], 1/5]);
省5
799: 04/29(火)10:38 ID:pY4WJf3b(2/4) AAS
questtion="良品サイコロは１の目の出る確率の95％信頼区間が[1/7,1/5]に収まるサイコロと定義する。それ以外はFランサイコロと呼ぶ。
チンパンフェチが発狂してサイコロを齧ってしまった。
このサイコロが良品のままかFランかを調べたい。
（１）２回投げたらどちらも１の目がでたときこのサイコロがFランである確率を求めよ。
（２）最初から何回１の目が続いてでたらFランである確率が0.5を超えるか？
計算に必要な条件は適宜決定してよい。
"
省7
800(1): 04/29(火)18:09 ID:pY4WJf3b(3/4) AAS
library(RcppAlgos)
library(fmsb)
library(matrixStats)
N <- 50
# Generate all combinations
cm <- comboGeneral(0:N, 3, repetition = FALSE)
# Pre-allocate
省38
801(1): 04/29(火)19:11 ID:A0Mypkqe(1) AAS
>>800
高校生が解く証明問題が解けずにコソコソ書き込みかよ
802: 04/29(火)21:05 ID:pY4WJf3b(4/4) AAS
options(warn = -1)
library(RcppAlgos)
N=50
alpha=0.01
cm=comboGeneral(0:N,3,repetition=FALSE)
f=\(x,Yates=FALSE){
n=rep(N,3)
省24
803(1): 04/30(水)02:47 ID:wedVH8wl(1/10) AAS
>>801
スレタイ読んだら。
Fラン卒は別スレだぞ。
804: 04/30(水)04:58 ID:wedVH8wl(2/10) AAS
options(warn = -1)
alpha=0.05
sim=\(N=100){
A=sample(1:(N-2),1) # Ensure A leaves enough room for B and C
remaining = N - A
if(remaining > 1){
B=sample(1:(remaining-1),1) # Ensure B is at least 1
省20
805: 04/30(水)05:14 ID:wedVH8wl(3/10) AAS
options(warn = -1)
alpha=0.05
sim=\(N=100){
A=sample(1:(N-2),1) # Ensure A leaves enough room for B and C
remaining = N - A
if(remaining > 1){
B=sample(1:(remaining-1),1) # Ensure B is at least 1
省25
806: 04/30(水)06:57 ID:wedVH8wl(4/10) AAS
おっしゃる通り、カイ二乗検定で連続性補正を外せば、より小さな差でも有意になりやすいため、ご要望のデータを作成できる可能性は高まります。しかし、Fisherの正確確率検定は、まさにその正確性ゆえに、サンプルサイズが小さい場合や比率の差が小さい場合に、p値が離散的になり、Bonferroni補正のような厳しい多重比較補正を乗り越えて有意差を示すのが難しい傾向があります。
Fisherの正確確率検定は、周辺度数を固定した条件下での確率に基づいてp値を計算するため、どうしても「わずかな差」が有意になりにくいという特性があります。特に、全体で有意差が出ない程度に比率の差を抑えようとすると、ペアワイズでも同様に差が小さくなり、Bonferroni補正によって有意水準が厳しくなるため、有意差を検出するのが非常に困難になります。
連続性補正は、カイ二乗分布の連続近似を離散的なデータに適用する際の誤差を小さくするための調整ですが、これを外すことで、p値が小さくなりやすくなります。一方、Fisherの正確確率検定はその性質上、近似を用いないため、連続性補正のような概念がありません。
結論として、ご指摘の通り、「全体のFisherの正確確率検定で有意差がないのに、ペアワイズなFisherの正確確率検定だとどれかに有意差がある（Bonferroni補正あり）」という条件を満たすデータを作成するのは、統計的な制約から非常に困難であると言わざるを得ません。
やっぱり、Bonferroniの壁はFisherでは乗り越えられようだ。
807: 04/30(水)07:36 ID:IZhDMqNd(1/2) AAS
>>803
スレタイ読めないのはアンタだぞマヌケ
808(1): 04/30(水)08:07 ID:wedVH8wl(5/10) AAS
ここはFランアクセス禁。
809: 04/30(水)08:07 ID:wedVH8wl(6/10) AAS
options(warn = -1)
alpha <- 0.05
sim_fisher <- function(N = 100) {
# Function to simulate data for three groups and perform Fisher's exact test.
# N: Total sample size.
# Determine sample sizes for each of the three groups.
A <- sample(1:(N - 2), 1) # Ensure A leaves enough room for B and C.
省31
810: 04/30(水)08:07 ID:wedVH8wl(7/10) AAS
# Find data where the overall Fisher's exact test is not significant,
# but at least one pairwise Fisher's exact test (with Bonferroni correction) is significant.
res_no_overall_sig_pairwise_sig <- NULL
while (is.null(res_no_overall_sig_pairwise_sig) || res_no_overall_sig_pairwise_sig$fisher_pg > alpha || res_no_overall_sig_pairwise_sig$min_pairwise_p_bonf > alpha) {
res_no_overall_sig_pairwise_sig <- sim_fisher()
}
cat("Data where overall Fisher's test is not significant, but pairwise is:\n")
省11
上下前次1-新書関写板覧索設栞歴
あと 83 ﾚｽあります
ｽﾚ情報 赤ﾚｽ抽出 画像ﾚｽ抽出 歴の未読ｽﾚ AAｻﾑﾈｲﾙ
ぬこの手 ぬこTOP 0.027s