我现在想用一个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行输出呢?谢谢!