我用的是knitr+RStudio组合，下面是一段代码：

<<echo=TRUE>>=<br />
library(DAAG)<br />
data(ironslag)<br />
attach(ironslag)<br />
n <- nrow(ironslag)<br />
e1 <- e2 <- e3 <- e4 <- numeric(ncol(combn(n,2)))</p>
<p>for (k in 1:ncol(combn(n,2))) {<br />
y <- magnetic[-combn(53,2)[,k]]<br />
x <- chemical[-combn(53,2)[,k]]<br />
J1 <- lm(y ~ x)<br />
yhat1 <- J1$coef[1] + J1$coef[2] * chemical[combn(53,2)[,k]]<br />
e1[k] <- sum((magnetic[combn(53,2)[,k]] - yhat1)^2)</p>
<p>J2 <- lm(y ~ x + I(x^2))<br />
yhat2 <- J2$coef[1] + J2$coef[2] * chemical[combn(53,2)[,k]] +<br />
J2$coef[3] * chemical[combn(53,2)[,k]]^2<br />
e2[k] <- sum((magnetic[combn(53,2)[,k]] - yhat2)^2)</p>
<p>J3 <- lm(log(y) ~ x)<br />
logyhat3 <- J3$coef[1] + J3$coef[2] * chemical[combn(53,2)[,k]]<br />
yhat3 <- exp(logyhat3)<br />
e3[k] <- sum((magnetic[combn(53,2)[,k]] - yhat3)^2)</p>
<p>J4 <- lm(log(y) ~ log(x))<br />
logyhat4 <- J4$coef[1] + J4$coef[2] * log(chemical[combn(53,2)[,k]])<br />
yhat4 <- exp(logyhat4)<br />
e4[k] <- sum((magnetic[combn(53,2)[,k]] - yhat4)^2)<br />
}<br />
print(c(mean(e1), mean(e2), mean(e3), mean(e4)))<br />
@

</p>

最后那个print函数还没打印出来程序运行的结果，knitr就帮我把pdf文档编译出来了。。结果就是pdf中print下面没有输出。。但是我用Sweave编译的话，又可以输出结果了，但是鉴于文档不太美观，所以想向大家请教问题的解决方案。