高校数学の質問スレ（医者・東大卒専用） Part438

高校数学の質問スレ（医者・東大卒専用） Part438 (979ﾚｽ)
上下前次1-新
通常表示 512ﾊﾞｲﾄ分割ﾚｽ栞
抽出解除ﾚｽ栞

766: １３２人目の素数さん [sage] 2025/03/25(火) 19:45:01.28 ID:JvJVqzfB

分け方の最初の30通りを列挙
[1,] 赤   |    赤   |    赤   赤   黒   黒   黒   白    白    青
[2,] 赤   |    黒   |    赤   赤   赤   黒   黒   白    白    青
[3,] 赤   |    白   |    赤   赤   赤   黒   黒   黒    白    青
[4,] 赤   |    青   |    赤   赤   赤   黒   黒   黒    白    白
[5,] 黒   |    黒   |    赤   赤   赤   赤   黒   白    白    青
[6,] 黒   |    白   |    赤   赤   赤   赤   黒   黒    白    青
[7,] 黒   |    青   |    赤   赤   赤   赤   黒   黒    白    白
[8,] 白   |    白   |    赤   赤   赤   赤   黒   黒    黒    青
[9,] 白   |    青   |    赤   赤   赤   赤   黒   黒    黒    白
[10,] 赤   |    赤   赤   |    赤   黒   黒   黒   白    白    青
[11,] 赤   |    赤   黒   |    赤   赤   黒   黒   白    白    青
[12,] 赤   |    赤   白   |    赤   赤   黒   黒   黒    白    青
[13,] 赤   |    赤   青   |    赤   赤   黒   黒   黒    白    白
[14,] 赤   |    黒   黒   |    赤   赤   赤   黒   白    白    青
[15,] 赤   |    黒   白   |    赤   赤   赤   黒   黒    白    青
[16,] 赤   |    黒   青   |    赤   赤   赤   黒   黒    白    白
[17,] 赤   |    白   白   |    赤   赤   赤   黒   黒    黒    青
[18,] 赤   |    白   青   |    赤   赤   赤   黒   黒    黒    白
[19,] 黒   |    赤   赤   |    赤   赤   黒   黒   白    白    青
[20,] 黒   |    赤   黒   |    赤   赤   赤   黒   白    白    青
[21,] 黒   |    赤   白   |    赤   赤   赤   黒   黒    白    青
[22,] 黒   |    赤   青   |    赤   赤   赤   黒   黒    白    白
[23,] 黒   |    黒   黒   |    赤   赤   赤   赤   白    白    青
[24,] 黒   |    黒   白   |    赤   赤   赤   赤   黒    白    青
[25,] 黒   |    黒   青   |    赤   赤   赤   赤   黒    白    白
[26,] 黒   |    白   白   |    赤   赤   赤   赤   黒    黒    青
[27,] 黒   |    白   青   |    赤   赤   赤   赤   黒    黒    白
[28,] 白   |    赤   赤   |    赤   赤   黒   黒   黒    白    青
[29,] 白   |    赤   黒   |    赤   赤   赤   黒   黒    白    青
[30,] 白   |    赤   白   |    赤   赤   赤   黒   黒    黒    青

http://rio2016.5ch.net/test/read.cgi/math/1723152147/766

816: １３２人目の素数さん [sage] 2025/05/01(木) 09:09:32.28 ID:L1qIlz9/

Stanで作ったらコンパイルに時間がかかる。簡単なモデルはJAGSの方がいい。離散変数も扱えるし。
# JAGS model
library(rjags)

# Prepare the data
outcome_data <- c(rep(1, 17), rep(2, 21), rep(3, 15), rep(4, 21), rep(5, 20), rep(6, 6))
N <- length(outcome_data)
data_jags <- list(outcome = outcome_data, N = N)

# Initial values (adjust as needed)
inits_jags <- list(
list(alpha = rep(1, 6), eta = 1),
list(alpha = runif(6, 0.1, 2), eta = 5)
)

# Compile the model
model_jags <- jags.model(
file = "hierarchical_dice_model.jag",
data = data_jags,
n.chains = 2,
n.adapt = 1000
)

# Sampling
samples_jags <- coda.samples(
model = model_jags,
variable.names = c("prob_simplex", "alpha", "eta"),
n.iter = 4000
)

# Summary of the results
cat("\nJAGS Sampling Results Summary:\n")
summary(samples_jags)

# Extract posterior samples (prob_simplex)
prob_simplex_posterior_jags <- as.matrix(samples_jags[, grep("prob_simplex", varnames(samples_jags))])
head(prob_simplex_posterior_jags)

# Plotting (example: posterior distribution of probabilities for each outcome)
cat("\nPosterior Distribution Plots for Each Outcome:\n")
par(mfrow = c(2, 3))
for (i in 1:6) {
plot(prob_simplex_posterior_jags[, i], type = "l", main = paste("Prob[", i, "]"), xlab = "Iteration", ylab = "Probability")
abline(h = 1/6, col = "red", lty = 2)
}

dice_prob_mean=prob_simplex_posterior_jags
colors <- c("skyblue", "lightcoral", "lightgreen", "gold", "lightsalmon", "lightcyan")
for (i in 1:ncol(dice_prob_mean)) {
BEST::plotPost(dice_prob_mean[, i],
compVal=1/6,
xlab=paste("pip ", i),
xlim=c(0, 0.4),
main="",
col=colors[i],
border="black")
}

http://rio2016.5ch.net/test/read.cgi/math/1723152147/816

上下前次1-新書関写板覧索設栞歴

ｽﾚ情報赤ﾚｽ抽出画像ﾚｽ抽出歴の未読ｽﾚ AAｻﾑﾈｲﾙ

ぬこの手ぬこTOP 1.664s*