Hello,

I am trying to generate a plot legend that contains calculated summary 
statistics, one statistic is R^2.  I have tried several variations using 
the commands "expression" and "bqoute" as stated on the R help pages.  I 
have not been able to get the R^2 super script correct along with the 
calculated statistics.

I provided an example below, what I want is the legend (wdt_leg and 
spas_leg) as below but with the R^2 as superscript.

# Example Data
x=c(1,2,3,4);  y=c(1,2,3,4); z=c(1.25,1.5,2.5,3.5)

# first stats based on data, used to populate legend
wdt_n = 50;  wdt_mbias = 0.58
wdt_mae = 2.1;  wdt_R2 = 0.85
# second stats based on data, used to populate legend
spas_n = 50; spas_mbias = 0.58
spas_mae = 2.1; spas_R2 = 0.85

# create legend
wdt_leg <- paste("WDT (N = ", wdt_n,", Bias = ",wdt_mbias,", MAE = 
",wdt_mae, ", R2 = ", wdt_R2, ")", sep="")
spas_leg <- paste("SPAS (N = ", spas_n,", Bias = ",spas_mbias,", MAE = 
",spas_mae, ", R2 = ", spas_R2, ")", sep="")

# create plot
plot(x,y, col="red1", pch=1); lines(x,z, type="p", col="green4",pch=3)
leg.txt <- c(spas_leg, wdt_leg)
legend("topleft", legend = leg.txt, col=c("red1","green4"), pch=c(1,3), 
cex=0.85)

Any help/suggestions on this would be great.  Thanks for your time in 
advance.
Doug

-- 
---------------------------------
Douglas M. Hultstrand, MS
Senior Hydrometeorologist
email: dmhul...@metstat.com
web: http://www.metstat.com
---------------------------------


        [[alternative HTML version deleted]]

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to