我现在想用一个legend命令输出一个scatterplot的correlation coefficient 的ρ和p。输出想要分2行，一行是ρ，一行是p。
使用了Hmisc。命令如下：
r0=round(rcorr(x,y)[[1]][1,2],4)
plot(x,y,main="x vs. y", xlab="x",ylab="y")
legend("topright", legend=c(paste("rho =",r0), paste("p =",round(rcorr(x,y)$P[1,2],4))),col=c(1,2),inset=0.02)
我想把那个rho变成希腊字母ρ，但是试了很多次都不对。如果rho正确了，后面的那个数字输不出来了。
我试过：
legend("topright", legend=c(expression(paste(rho," =",r0)), paste("p =",round(rcorr(x,y)$P[1,2],4))),col=c(1,2),inset=0.02)
这样输出，前面的ρ是出来了，但是r0只是字母，不是数字。
还试过：
  legend("topright", legend=c(paste(expression(paste(rho," =")),r0), paste("p =",round(rcorr(x,y)$P[1,2],4))),col=c(1,2),inset=0.02)
这样后面数字出来了，但是前面的ρ不出来。想请教各位，如何写才能让2个都出来。而且分成2行输出呢？谢谢！