search for: qqnorm

I want to qqnorm every column in a table. When I try the first column using qqnorm(table$column1), it worked. But when I use qqnorm(table[1]), it tells me "Error in stripchart(x1, ...) : invalid plotting method". What happen? How can I make a function that qqnorms every column? thanks a lot. -jion...

I know this topic has had plenty of discussion in the last couple of days, but.... I've been trying to compare the effects of different fitted methods for systems of equations (OLS, SUR, 2SLS, 3SLS ) and would like to compare the residual plots (easy) and the qqnorm/qqplot of the fits for the different fitted methdos. For example, qqnorm( residuals( lm( q ~ p + f + a ) ) ) par( new = TRUE ) qqnorm( residuals( tsls( q ~ p + d, ~ d + f + a ) ) ) The resulting plot contains the two qqnorm plots, and the x axis stays the same. The y axis ("Sample Quantiles&...

Hi all, I have just started using R to produce qqnorm plots. I am trying to switch the x and y axes so that the theoretical values are plotted on the y axis and my data on the x axis. Can anyone help me with this? Thanks -- View this message in context: http://www.nabble.com/Flipping-axes-of-qqnorm-tp23248007p23248007.html Sent from the R help maili...

Referring to "Using R for Data Analysis and Graphics" by J H Maindonald, and available from the R site, I found the example on p.30 non-working: > stem(qqnorm(possum$hdlngth)) Error in stem(qqnorm(possum$hdlngth)) : 'x' must be numeric Since qqnorm(possum$hdlngth) plots, and > class(possum$hdlngth) [1] "numeric" , the problem must be here: > class(qqnorm(possum$hdlngth)) [1] "list" Does 'stem' require numerical...

# I can not draw a 45 degree line on a qqnorm() plot, jj <- sample(c(1:100), 10) qqnorm(jj) abline() don't work. Thank you.

Hello R friends, I'm wondering why I get funny qqnorm() results. It seems that they should all be reflected in the normal qqline(). For instance: if I qqnorm() bimodal or uniform data I get a sigmoidal in which the qqnorm() points lie above the qqline() at -ve theoretical quantiles, and the qqnorm() points lie below the qqline() at +ve theoretical q...

Hello, I am trying to do some plotting to check random effect assumptions for a model I fit using lme. I want to use qqnorm and pairs (similarly to examples given in Pinheiro & Bates p. 188), but it's not working. Here's some relevant code and the error message: library(nlme) data(Machines) m1 <- lme(fixed=score~Machine,random=~1|Worker/Machine, data=Machines) qqnorm(m1, ~ranef(.)) Error in rep(names(f...

Dear List, invoking qqnorm-plots in Version 63.3 produces funny things: using the option `type="s"´ on qqnorm should give a nice *line* of observed quantiles. Now, the line is walking along in order to the points index instead from lowest to highest, wich makes funny slopes. try x <- table(rnorm(1000) # or sim...

Hi all, I am an R newbie, and I have a question about scripting. I have the following lines which I want to put int\ o a script which I can call from the shell of a Mac/Linux machine : myrns <- read.csv(file="/Users/vihan/test.csv",sep="",header=FALSE) qqnorm(myrns) qqline(myrns) This works fine on an interactive R session. However, as I understand the basic framework of an R script for executing on the shell to get a pdf output it is something like : #!/usr/bin/Rscript ? list inputs? pdf("plotexample.pdf", width=7, height=9) plot(...)...

Dear list-reader, by running the following script: library(zoo) sessionInfo() search() packageDescription("zoo") data(EuStockMarkets) dax <- as.zoo(EuStockMarkets[1:10, "DAX"]) daxr <- diff(log(dax)) identical(as.vector(qnorm(daxr)), qnorm(coredata(daxr))) qqnorm(coredata(daxr)) qqnorm(daxr) qqnorm() produces an error: > qqnorm(daxr) Fehler in if (xi == xj) 0L else if (xi > xj) 1L else -1L : Argument hat L?nge 0 Calls: qqnorm ... xtfrm.default -> as.vector -> rank -> switch -> .gt Ausf?hrung angehalten However, this is not witnesse...

Hi all, Does anyone know how to plot overlapping qqnorm plots on the same window? Suppose I have data in the vector x and y: qqnorm(x) lines(qqnorm(y)) I though these two lines will do the job... However, lines doesn't seem to work. Anyways, thanks in advance! -- View this message in context: http://www.nabble.com/Plotting-two-qqnorm-plots%3A-tp2...

Hello, I have a doubt about the plotting position of the theoretical quantile for the qqnorm command in R. Let F be the theoretical distribution of Y, we observed a sample of size n, y1,y2, ..., yn. We then sort it and comspare these empirical quantiles to the expected ones from F. For the plotting poition, there are several options: 1. i/(n+1) 2. (i-.375)/(n+.25) 3. (i- .3175)/ (n + .36...

Does anyone know how to extract residuals in lmer? Here's the error I get: > crop.lme=lmer(response~variety*irrigation*pesticide+(1|rep)+(1|rep: pesticide)+(1|rep:pesticide:irrigation), crop.data) > qqnorm(crop.lme) Error in qqnorm.default(crop.lme) : y is empty or has only NAs > resid(crop.lme) NULL Thanks! --Jake

...nput data are scaled so that they are insignificant compared to 1, say, then you get essentially the "theoretical quantiles" ie the "x" component of the list back from l$x - l$y. l$x is basically a sample from a normal(0,1) distribution so they do line up perfectly in the second qqnorm(). Is that what's happening? Reid Huntsinger -----Original Message----- From: r-help-bounces at stat.math.ethz.ch [mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of WeiWei Shi Sent: Thursday, April 28, 2005 1:38 PM To: Vincent ZOONEKYND Cc: R-help at stat.math.ethz.ch Subject: [R] ha...

...340631367 0.997186769599305 0.991632565640424 0.984047197106486 0.99225943762649 1.00555642128421 0.993725402926564 .... the data is saved in a file called f392.txt. I used the following codes to play around :) k<-read.table("f392.txt", header=F) # read into k kk<-k[[1]] l<-qqnorm(kk) diff=c() lenk<-length(kk) i=1 while (i<=lenk){ diff[i]=l$y[i]-l$x[i] # save the difference of therotical quantile and sample quantile # remember, my sample mean is around 1 while the therotical one, 0 i<-i+1 } hist(diff, breaks=300) # analyze the distr o...

Just found because an old (written for S+) function of mine did label plots wrongly with R. Example --- inspired from example(qqnorm) --- data(precip) qqnorm(precip, ylab = "Precipitation [in/yr] ...") qqnorm(precip, ylab = "Precipitation [in/yr] ...", datax = TRUE) this is all "fine" -- `datax' is a switch that just switches the axes but keeps the meaning of "ylab labels the data a...

Hi, When qqnorm on a vector of length 10M+ I get a huge pdf file which cannot be loaded by acroread or evince. Any suggestions? (apart from sampling the data). Thanks. -- Sam Steingold (http://sds.podval.org/) on Ubuntu 11.10 (oneiric) X 11.0.11004000 http://mideasttruth.com http://honestreporting.com http://came...

Dear all, I have two variables called c2 and c3 and want to plot these variables in the same qqnorm-plot with two different symbols or colors to distinguish them so I can easily compare the variables aginst each other. How can I do in R? I only manage to do two separated qqnorn-plots. Thanks for your help, All the best, Jenny c2= -0.1545775144 -0.0601161235 -0.145471...

Do you know of a reference that discusses alternative choices for plotting positions for a normal probability plot? The documentation for qqnorm says it calls ppoints, which returns qnorm((1:m-a)/(m+1-2*a)) with "a" = ifelse(n<=10, 3/8, 1/2)? The help pages for qqnorm and ppoints just refer to Becker, Chambers and Wilks (1988) The New S Language (Wadsworth & Brooks/Cole), and I couldn't find any discussion of this....

...ses if I keep bugging everyone in this list, but I am a newbie, and tend to find some behaviour that looks unexpected to me; and I would really appreciate to be pointed to some location that allows me to understand more about this software. Here is my next question: > par(mfrow=c(1,1)) > qqnorm(rnorm(20)) > qqmath(rnorm(20)) > par(mfrow=c(3,4)) > for(i in 1:12)qqnorm(rnorm(20)) Until here everything works as expected, and the last line prints 12 samples of qqnorm. However, > for(i in 1:12)qqmath(rnorm(20)) is doing nothing at all. As I wrote, I'd really appreciate th...