Re: [R] adding statistical output to a plot

2010-05-30 Thread Jim Lemon
On 05/29/2010 10:11 PM, dunner wrote: ... Perhaps "apply" simply doesn't deal with the colname and I have to try a completely different method? Or will apply pass enough information to the included function (mtab.norm) to allow me to annotate the plot? It doesn't seem to, as I tried this on a d

Re: [R] adding statistical output to a plot

2010-05-29 Thread Tal Galili
Hi Dunner, The way I usually resolve this is by extracting the names with colnames and then reinsert them to the new object. BTW, why are you trying to reproduce the image from minitab? what features in it do you find useful ? Best, Tal Contact Details:-

Re: [R] adding statistical output to a plot

2010-05-29 Thread David Winsemius
On May 29, 2010, at 8:11 AM, dunner wrote: Thanks to both of you, I used the below code instead, which has been cleaned and pruned a little. It produces a trellis of qq plots with the proper annotations as per Minitab from a dataframe passed to mtab.matrix(). However, I can't get indivi

Re: [R] adding statistical output to a plot

2010-05-29 Thread dunner
Thanks to both of you, I used the below code instead, which has been cleaned and pruned a little. It produces a trellis of qq plots with the proper annotations as per Minitab from a dataframe passed to mtab.matrix(). However, I can't get individual variable names (these are column names in the pa

Re: [R] adding statistical output to a plot

2010-05-29 Thread Jim Lemon
On 05/29/2010 07:39 PM, dunner wrote: I have written a function to emulate minitab's QQ plotting output (with SW test and AD test results on the graph): ... However, I get the following message when I try to apply it: ab<-rnorm(1:100, 2) aa<-log(ab) Warning message: In log(ab) : NaNs produce

Re: [R] adding statistical output to a plot

2010-05-29 Thread Tal Galili
Two things: 1) doing log on rnorm, will get you NaN's because you are doing log on a negative number (some of the time). 2) Running > traceback() 5: sort(y) 4: qqplot(x, pch = 16, cex = pt.ht, col = palette()[c(1, 2)], xlab = paste(colnames(x)), ylab = "Normal Quantile-Quantile Plot", main

[R] adding statistical output to a plot

2010-05-29 Thread dunner
I have written a function to emulate minitab's QQ plotting output (with SW test and AD test results on the graph): mtab.norm<-function(x) { library(nortest) library(lattice) x<-as.numeric(x) x<-as.vector(x) plot.ht<-4.6 plot.wd<-4.6 pt.ht=plot.ht/5 txt.sz<-(plot.ht/7.5) X11(width=plot.wd, height=