hi,
i am using this script on our webserver to plot data / draw ellipses
for a flex movie...
would like to add the areas of the ellipses to another legend just
like the means for
the distance.
the data looks like: x (dispersion), y (distance), Club (name of
the club for grouping)
thanks, sam
## usage: Rscript nikeScatter.R infile outfile level1 level2 minX maxX
minY maxY
## load the library
library(car)
## get the args
args = commandArgs(TRUE);
## get the data format is: x, y, group
shotData <- read.table(args[1],
sep = ",",
header=T)
## setup the pdf file name
pdf(file = args[2])
## do the plot
scatterplot(y ~ x | Club, data=shotData,
ellipse = TRUE,
levels = c(as.numeric(args[3]), as.numeric(args[4])),
robust = TRUE,
boxplots = "",
xlab = "Dispersion",
ylab = "Distance",
xlim = c(as.numeric(args[5]), as.numeric(args[6])),
ylim = c(as.numeric(args[7]), as.numeric(args[8])),
smooth = FALSE,
reg.line = FALSE,
legend.plot = FALSE)
## do the custom legends
n.groups <- length(levels(shotData$Club))
col=rep(palette(), length.out=n.groups + 1)
pch=1:n.groups
usr <- par('usr')
distance <- sapply(levels(shotData$Club),function(x)
{mean(shotData[shotData$Club==x,"y"])})
dispersion <- sapply(levels(shotData$Club),function(x)
{mean(shotData[shotData$Club==x,"x"])})
legend(usr[1] + 7, usr[4], legend=levels(shotData$Club), col=col[2:
(n.groups+1)], pch=pch, pt.cex=.85, cex=.85)
legend("topright", legend=format(distance, digits=1, nsmall=3),
col=col[2:(n.groups+1)], pch=pch, pt.cex=.85, cex=.85)
______________________________________________
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.
