Question

我正在编写一个简短的Sweave文档，输出到Beamer演示文稿中，我使用sagetex包来解决beta二项分布中两个参数的等式，我需要将参数值分配给R会话所以我可以对这些值进行额外的处理。以下代码摘录显示了我如何与sage交互：

<<echo=false,results=hide>>=
mean.raw <- c(5, 3.5, 2)
theta <- 0.5  
var.raw <- mean.raw + ((mean.raw^2)/theta)
@

\begin{frame}[fragile]
\frametitle{Test of Sage 2}
\begin{sagesilent}
var('a1, b1, a2, b2, a3, b3')
eqn1 = [1000*a1/(a1+b1)==\Sexpr{mean.raw[1]}, ((1000*a1*b1)*(1000+a1+b1))/((a1+b1)^2*(a1+b1+1))==\Sexpr{var.raw[1]}]
eqn2 = [1000*a2/(a2+b2)==\Sexpr{mean.raw[2]}, ((1000*a2*b2)*(1000+a2+b2))/((a2+b2)^2*(a2+b2+1))==\Sexpr{var.raw[2]}]
eqn3 = [1000*a3/(a3+b3)==\Sexpr{mean.raw[3]}, ((1000*a3*b3)*(1000+a3+b3))/((a3+b3)^2*(a3+b3+1))==\Sexpr{var.raw[3]}]
s1 = solve(eqn1, a1,b1)
s2 = solve(eqn2, a2,b2)
s3 = solve(eqn3, a3,b3)
\end{sagesilent}

Solutions of Beta Binomial Parameters:
\begin{itemize}
\item $\sage{s1[0]}$
\item $\sage{s2[0]}$
\item $\sage{s3[0]}$
\end{itemize}
\end{frame}

一切都编译得很好，在那张幻灯片中，我能够看到该项目化列表中三个方程式各个参数的解决方案（例如，该beamer幻灯片的逐项列表中的第一个项目输出为{{1} （我无法发布beamer幻灯片的图像，但我希望你能得到这个想法）。

我想将参数值[a1=(328/667), b1=(65272/667)]保存到R对象中，以便我可以在模拟中使用它们。我在sagetex包中找不到任何关于如何将sage命令的输出保存到变量以供其他程序使用的文档（在本例中为R）。关于如何将这些值纳入R的任何建议？

Answer 1

哇，你真的在混合两个世界;）

我能给你的唯一想法是solve命令的“solution_dict = True”参数。然后你会得到一个Python字典，它可能会帮助你输出值。但我不知道Sweave究竟做了什么以及该过程的哪一步重写了什么。

一般来说，如果你只在sagetex中编写它并通过rpy2 Python包装器调用R可能会更好。但这可能对你来说太多了 - 可能只是为了一张幻灯片，然后通过一些pdf合并将它们拼接在一起？

通过sagetex和Sweave将sage变量值分配给R对象

1 个答案: