Hello, I think you are trying to combine too many things into one.
for(i in 1:100){ > points(cor(v$logA[v$id>i], v$logB[v$id>i], use="complete.obs"), i)) > } > (i created a plot object already) For points() you need to have x and y, but you are giving the output of cor() This will plot your data and print the correlations on the plot (lathough maybe not in the best place): for(i in 1:100){ points(v$logA[v$id>i], v$logB[v$id>i], pch=i) text(0, i, cor(v$logA[v$id>i], v$logB[v$id>i], use="complete.obs")) } Hopefully that gives you an idea of how to get the output you actually want. Sarah -- Sarah Goslee [[alternative HTML version deleted]] ______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html