search for: qqnorm

Displaying 20 results from an estimated 206 matches for "qqnorm".

Did you mean: qnorm
2007 Jun 29
2
why this doesn't work 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...
2001 Nov 27
2
overlaying qqnorm plots...
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&...
2009 Apr 26
3
Flipping axes of qqnorm
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...
2010 Apr 15
1
classes and functions for qqnorm and stem
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...
2005 Apr 03
2
how to draw a 45 degree line on qqnorm() plot?
# 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.
2000 Sep 21
2
qqnorm(), is it "backwards"?
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...
2009 Apr 23
1
qqnorm.lme & pairs.lme
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...
1999 Mar 16
1
qqnorm in R-0.63.3
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...
2011 Jul 25
1
scripting qqplot and qqnorm
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(...)...
2008 Oct 22
1
R 2.8.0 qqnorm produces error with object of class zoo?
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...
2009 Sep 08
1
Plotting two qqnorm plots:
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...
2007 Feb 08
1
the plotting position of theoretical quantile for qqnorm
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...
2005 Apr 25
2
residuals in lmer
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
2005 Apr 28
3
have to point it out again: a distribution question
...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...
2005 Apr 22
1
an interesting qqnorm question
...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...
2003 Oct 15
0
qqnorm(*, datax=TRUE, xlab,ylab) -- S+ compatibility problem
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...
2011 Dec 21
4
qqnorm & huge datasets
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...
2006 Dec 08
1
How to plot two variables in the same qqnorm-plot?
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...
2006 Apr 13
2
Plotting positions in qqnorm?
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....
2010 Apr 15
2
Consistent behaviour of for-loop
...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...