similar to: multiple plot devices

Displaying 20 results from an estimated 2000 matches similar to: "multiple plot devices"

2000 Dec 11
1
qqline (PR#764)
I think qqline does not do exactly what it is advertised to do ("`qqline' adds a line to a normal quantile-quantile plot which passes through the first and third quartiles."). Consider the graph: tmp <- qnorm(ppoints(10)) qqnorm(tmp) qqline(tmp) The line (which I expected go through all the points), has a slightly shallower slope than does the points plotted by qqnorm. I think
2010 Aug 06
1
qqline error: Error in sort.list(x, partial = unique(c(lo, hi))) : 'x' must be atomic for 'sort.list' Have you called 'sort' on a list?
Hi all, I'm modeling using lme in the nlme package. qqnorm makes plots just find, but when I try to add a line with qqline, I get the following error: "Error in sort.list(x, partial = unique(c(lo, hi))) : 'x' must be atomic for 'sort.list' Have you called 'sort' on a list?" For example > modelincr10<-lme(lruin~can.pos.num, data=rwushi30, random=(~1|
2011 Apr 30
4
QQ plot for normality testing
Hi all, I am trying to test wheater the distribution of my samples is normal with QQ plot. I have a values of water content in clays in around few hundred samples. Is the code : qqnorm(w) #w being water content qqline(w) sufficient? How do I know when I get the plots which distribution is normal and which is not? Thanks, m [[alternative HTML version
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 quantiles. Yet I expect such platykurtic
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.
2009 Nov 02
7
qqplot
Hi, We could use qqplot to see how two distributions are different from each other. To show better how they are different (departs from the straight line), how is it possible to plot the straight line that goes through them? I am looking for some thing like qqline for qqnorm. I thought of abline but how to determine the slope and intercept? Best wishes, Carol
2010 May 16
2
Box-Cox Transformation: Drastic differences when varying added constants
Dear experts, I tried to learn about Box-Cox-transformation but found the following thing: When I had to add a constant to make all values of the original variable positive, I found that the lambda estimates (box.cox.powers-function) differed dramatically depending on the specific constant chosen. In addition, the correlation between the transformed variable and the original were not 1 (as I
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
2006 Feb 01
1
Difficulty with qqline in logarithmic context
Hi, R friends. I had some difficulty with the following code: qqnorm(freq, log='y') qqline(freq) as the line drawn was seemingly random. The exact data I used appears below. After wandering a bit within the source code for "abline", I figured out I should rather write: qqnorm(freq, log='y') par(ylog=FALSE) qqline(log10(freq)) par(ylog=TRUE)
2010 Dec 18
3
use of 'apply' for 'hist'
Hi all, ########################################## dof=c(1,2,4,8,16,32) Q5=matrix(rt(100,dof),100,6,T,dimnames=list(NULL,dof)) par(mfrow=c(2,6)) apply(Q5,2,hist) myf=function(x){ qqnorm(x);qqline(x) } apply(Q5,2,myf) ########################################## These looks ok. However, I would like to achieve more. Apart from using a loop, is there are fast way to 'add' the titles to be
2007 Oct 03
2
Ferret and Mongrel and redirect_to
Hello, I try and run mongrel 1.0.1 and ferret 0.11.4 on Windows XP. To break mongrel all I have to do is invoke a class with Ferret in it. Here''s my class: require ''ferret'' include Ferret class SearchIndex def initialize print "SearchIndex:Initialize\n"; end end When I put this in my controller, I break mongrel: searchIndex =
2010 Jun 24
4
Simple qqplot question
I am a beginner in R, so please don't step on me if this is too simple. I have two data sets datax and datay for which I created a qqplot qqplot(datax,datay) but now I want a line that indicates the perfect match so that I can see how much the plot diverts from the ideal. This ideal however is not normal, so I think qqnorm and qqline cannot be applied. Perhaps you can help? Ralf
2005 Oct 29
1
how to get colnames of a dataframe within a function called by 'apply'
Hello alltogether, how is it possible to assign the colnames of a data.frame to a function called by apply, e.g. for labeling a plot? Example: I want to plot several qqnorm-plots side by side and there should be a maintitle for each qqnorm-plot which is identical to the respective colname. I checked, but the column which is processed by the function called by apply does not contain a colname
2005 Apr 28
3
have to point it out again: a distribution question
Stock returns and other financial data have often found to be heavy-tailed. Even Cauchy distributions (without even a first absolute moment) have been entertained as models. Your qq function subtracts numbers on the scale of a normal (0,1) distribution from the input data. When the input data are scaled so that they are insignificant compared to 1, say, then you get essentially the
2008 Jul 06
1
Backgrounds in Multiple Plots made with "fig"
The following code was adapted from an example Vincent Zoonekynd gave on his web site http://zoonek2.free.fr/UNIX/48_R/03.html: n <- 1000 x <- rnorm(n) qqnorm(x) qqline(x, col="red") op <- par(fig=c(.02,.5,.5,.98), new=TRUE) hist(x, probability=T, col="light blue", xlab="", ylab="", main="", axes=F) lines(density(x),
2008 Aug 05
2
qqline function doesn't plot
I have a data vector x. When I try qqline(x) I get the following error: Error in int_abline(a = a, b = b, h = h, v = v, untf = untf, ...) : plot.new has not been called yet And a blank plot appears. Can anybody help? What am I doing wrong? Thanks, Scotty _________________________________________________________________ Contest [[alternative HTML version deleted]]
2012 Feb 10
1
Formatting Y axis.
I've looked around and I just can't find anything that will work for my needs. This is a bit of a 2 part question but pertaining to the same topic so bare with me. The first is with my qq plot. On the Y axis of my qq plot it'll have my sample quantities but because my data is log-normal it'll show numbers between 0 - 5 (depending on the data). I'd like to know how to get it,
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. -jiong The email message (and any attachments) is for the sole use of the intended
2000 Feb 24
2
(-1 as index) OR (envelope for QQ)
I'm new to R (and to S) and am wondering about code from pages 72 and 83 of MASS (Venables+Ripley, 3rd edition), to draw an envelope on a QQ plot. Copying from the book, I've got: #... code whose gist is "a.fit <- nls(..." num.points <- length(resid(a.fit)) qqnorm(residuals(a.fit)) # illustrate data-model residuals qqline(residuals(a.fit)) samp <-
2014 Jul 11
2
outliers (Marta valdes lopez)
Tu fichero tiene los decimales como puntos y no como comas como tu le indicas. Te dejo un ejemplo #--------------------------------------------------------------------------------------------------------------------- setwd(dir="c:/Users/usuario/Desktop/") library(outliers) filename<-"timediff.csv" time<-read.csv(filename, sep=";",header=TRUE,dec=".")