paste()、bquote()、as.expression()、c()、...のほぼすべての組み合わせを試しました...
plot(d[,"y"], type="l",xlim=c(1,n), ylim=c(min(d[,"y"]),max(d[,"y"])),
ylab="Y", xlab="T", main="ARMA(1,1)",
sub=c(as.expression(bquote(phi == .(coef_ar)),
as.expression(bquote(theta == .(coef_ma))))))
これは "phi = 0.5" (この場合はギリシャ記号) をプロットするだけで、2 番目の部分 (シータ) はプロットしません。誰でも私を助けてください!
ありがとう!