在 RcppArmadillo 中生成序列

Question

我找不到生成序列并在 RcppArmadillo 中获取列平均值的简单方法？

有什么简单的解决办法吗

colMeans(G[(N_burn+1):N_Gibbs,])

Answer 1

这是你的第三个问题，我告诉你去 Rcpp Gallery. The post on Armadillo subsetting has pretty much all you need, the Armadillo API documentation and the Rcpp documentation 剩下的：

• 犰狳使用 unsigned int 类型索引
• 所以在uint
中创建一个索引向量
• 子集原始向量（或矩阵列）
• 使用例如 Rcpp 的 "Sugar" 函数 mean() 将在输入中获取向量

大功告成。