高校数学の質問スレ（医者・東大卒専用） Part438 (991ﾚｽ)
上下前次1-新
抽出解除 ﾚｽ栞
192(2): 2024/08/29(木)17:42:28.62 ID:Mdm5M0Mq(1/2) AAS
>>190
算出するための関数を作るのが楽しいんだね。
それをシミュレーションで検証。
224: 2024/09/17(火)07:36:31.62 ID:X3jkfzLo(1) AAS
vonNeuman <- function(PDF,xmin=0,xmax=1){
N=1e6
ymax=max(PDF(seq(xmin,xmax,length=N+1)))
Ux=runif(N,xmin,xmax)
Uy=runif(N,0,ymax)
Rand=Ux[which(Uy<=PDF(Ux))]
hist(Rand,xlim=c(xmin,xmax),freq=FALSE,breaks=30,col=sample(colors(),1),main='')
curve(PDF,add=TRUE,lwd=2)
invisible(Rand)
}
vonNeuman(dnorm,-3,3)
vonNeuman(\(x)sin(x)/2,0,pi)
397: 2024/12/08(日)06:35:45.62 ID:+/92f6rf(1) AAS
論破されてることすら理解出来ない低知能クンだからね
医者東大コンプなのに高校数学すら分からなくて、頑張って勉強したwolframを使いたいけど、同じ問題でパラメーター変えることしか出来ない低脳なんでしょ。
689: 02/06(木)17:26:31.62 ID:r+g7n4SK(2/5) AAS
隔離病棟だよ
スレタイ読めない上に出題と質問の違いもわからずに発狂してる尿瓶ジジイID:rViXyt0wの
703: 02/08(土)06:03:19.62 ID:SqsAAYfd(1) AAS
毎日書き込むほど暇でないだけでは
休日も来てないのはPCを開かないからでしょ
高校数学の本スレ
誰も次を立ててないけど作っとく？
単発の質問スレが増えてきてるし
709: 02/08(土)12:31:03.62 ID:1THAJYf7(1) AAS
友達どころか精神科連れて行ってくれる家族もいなさそう
812: 04/30(水)08:12:39.62 ID:wedVH8wl(9/10) AAS
# Find data where the overall Chi-squared test is not significant (p > alpha),
# but at least one pairwise proportion test (with Bonferroni correction) is significant (p <= alpha).
res_no_overall_sig_pairwise_sig <- NULL
while (is.null(res_no_overall_sig_pairwise_sig) || res_no_overall_sig_pairwise_sig$chisq_pg > alpha || res_no_overall_sig_pairwise_sig$min_pairwise_p_bonf > alpha) {
res_no_overall_sig_pairwise_sig <- sim_chisq() # Keep simulating until the condition is met.
}
cat("Data where overall Chi-squared test is not significant, but pairwise proportion test is:\n")
print(res_no_overall_sig_pairwise_sig)
cat("\n")
# Find data where the overall Chi-squared test is significant (p <= alpha),
# but all pairwise proportion tests (with Bonferroni correction) are not significant (p > alpha).
res_overall_sig_no_pairwise_sig <- NULL
while (is.null(res_overall_sig_no_pairwise_sig) || res_overall_sig_no_pairwise_sig$chisq_pg < alpha || res_overall_sig_no_pairwise_sig$min_pairwise_p_bonf < alpha) {
res_overall_sig_no_pairwise_sig <- sim_chisq() # Keep simulating until the condition is met.
}
cat("Data where overall Chi-squared test is significant, but pairwise proportion test is not:\n")
print(res_overall_sig_no_pairwise_sig)
options(warn = 0)
857: 05/31(土)05:25:02.62 ID:jzcOJBMt(1) AAS
#' @title ベイズ事後確率計算関数
#' @description 帰無仮説と対立仮説の事後確率を計算
#' @param s 観測成功数（1の目が出た回数）
#' @param n 総試行回数
#' @param p0 帰無仮説の確率（例: 1/6）
#' @param prior 帰無仮説の事前確率（0～1）
#' @param alpha 対立仮説のベータ分布αパラメータ
#' @param beta 対立仮説のベータ分布βパラメータ
#' @return list(bf01=ベイズファクター, post=事後確率, method=使用手法)
calculate_posterior <- function(s, n, p0, prior, alpha, beta) {
# 入力検証
stopifnot(
s >= 0, n > 0,
p0 > 0 && p0 < 1,
prior >= 0 && prior <= 1,
alpha > 0, beta > 0
)
# 帰無仮説の尤度計算
m0 <- dbinom(s, n, p0)
# 対立仮説の周辺尤度計算（extraDistr有無で自動切替）
if (requireNamespace("extraDistr", quietly = TRUE)) {
m1 <- extraDistr::dbbinom(s, n, alpha, beta)
method <- "extraDistr::dbbinom()"
} else {
integrand <- function(p) dbinom(s, n, p) * dbeta(p, alpha, beta)
m1 <- integrate(integrand, 0, 1)$value
method <- "数値積分"
}
# ベイズファクターと事後確率計算（指定された式を使用）
bf01 <- m0 / m1
post <- m0 * prior / (m0 * prior + (1 - prior) * m1)
# 結果をリストで返す
list(
bf01 = bf01,
post = post,
method = method,
inputs = list(s = s, n = n, p0 = p0, prior = prior, alpha = alpha, beta = beta),
likelihoods = list(m0 = m0, m1 = m1)
)
}
903(1): 06/21(土)18:30:36.62 ID:o0NLVsTG(1) AAS
>>902
理由www
何処が傑作なんだよwww
簡単な総当たりで解ける問題なんか傑作とは言えないな
907: 06/21(土)20:12:22.62 ID:iVrSlwPC(2/2) AAS
>>906
話通じてない
日本語も不自由みたいだね
944: 06/23(月)21:38:33.62 ID:3REuqFh0(1) AAS
>>943
想定www
独りよがり思い込み評価とかゴミクズそのものwwww
上下前次1-新書関写板覧索設栞歴
ｽﾚ情報 赤ﾚｽ抽出 画像ﾚｽ抽出 歴の未読ｽﾚ
ぬこの手 ぬこTOP 0.056s