如何不进行模型拟合即可查看各类Stan/RStan操作的输出结果

Stan提供了大量可作用于数组、矩阵与混合数据类型的操作，实际开发时很多人会遇到拿不准函数实际执行效果、想快速看运行输出的需求。
比如拿到一个向量，想确认传入diag_matrix()后的返回结果，哪怕官方文档明确说明这个函数会返回以输入值为对角线的对角矩阵，开发时直接看实际输出的参考价值还是很高。不管是在R环境调用Stan函数，还是直接在Stan程序里，都有非常轻量的实现方式，不需要跑完整的模型采样流程，就能把结果直接打印到控制台，具体方案如下：

R 环境实现 #

借助cmdstanr包的函数暴露功能，可以直接把Stan的内置函数映射到R环境里调用，全程不需要走MCMC采样流程，步骤非常简单：

• 写一个最简的Stan测试文件，不需要写完整的data、parameters、model块逻辑
• 编译这个最简模型后，调用$expose_functions()把Stan函数挂载到R的全局环境
• 直接传参调用目标函数，运行一次就能在控制台看到输出

对应可运行代码：

library(cmdstanr)
# 生成最简测试模型
test_model_file <- write_stan_file("
model {}
")
# 编译模型
mod <- cmdstan_model(test_model_file)
# 暴露Stan内置函数到R全局环境
mod$expose_functions(global = TRUE)
# 测试目标函数
x <- c(1, 2, 3, 4, 5)
diag_matrix(x)

运行后直接就能在控制台看到生成的5阶对角矩阵，整个过程最快几秒就能跑完。

Stan 程序内实现 #

如果要直接在Stan代码运行时输出结果，用Stan内置的print()函数即可。如果只想让结果打印一次、不要每次采样迭代都重复输出，把打印逻辑放在transformed data块就行——这个块只会在模型初始化阶段执行一次，刚好匹配需求。
对应可运行代码示例：

data {
  int K;
  vector[K] x;
}
transformed data {
  matrix[K, K] test_output = diag_matrix(x);
  print("diag_matrix函数返回结果：", test_output);
}
model {
  // 其余模型正常逻辑即可，哪怕留空也能正常运行输出
}

注意不要把打印语句放在model、transformed parameters这类每次采样迭代都会执行的块里，不然会重复打印刷满控制台。要测试其他数组、矩阵操作的话，只需要把示例里的diag_matrix(x)替换成你要验证的目标函数即可，逻辑完全通用。

内容的提问来源于stack exchange，提问作者socialscientist