similar to: overlaying qqnorm plots...

I'm trying to build a flexible OLS/SUR/2SLS/3SLS package and I'm having trouble getting some information out of a formula. For example.... # set up the system of equations demand <- q ~ p + d supply <- q ~ p + f + a inst <- ~ d + f + a systemeq <- list( demand, supply ) ...blah, blah, blah... # get the number if instruments... # or build the "bigX"... #I

Hi all, I want to create a rather standard overlaid qqnorm plot on a single variable, with different subgroups of the same dataset plotted using different colors/symbols/etc. (I don't want side-by-side, rather different-colored curves on the same graph) I managed to do it rather tediously using "lapply" and "split", and wondered whether there is any single-command

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: >

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

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

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 similar and qqnorm(x,type="s") # in 0.63.2 and 63.3 Well, the

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,

I've been trying to figure out how to get the names of the parameter vector variables when inside the function that nlm calls to return the objective function value: knls <- function( theta, eqns, data, fitmethod="OLS", instr=NULL, S=NULL ) { ## print( names( theta ) ) # returns NULL ## get the values of the parameters for( i in 1:length( theta ) )

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

I've figured out how to use optim (barely) to estimate 2 parameter = weibull distributions. I can't get over how easy this is. What I need to = do is use a weight in the observations..... For example,=20 the tree diameters and weights are are=20 4.70 , 100 6.00, 98 7.10, 75.0 8.10, 86.3 8.60, 80.456 8.90, 20.5 9.50, 16.6 11.40, 12.657 11.80, 12.47 14.50,

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-tp25352893p25352893.html Sent from the R help

I thought I had plotted a map of oregon and washington and overlayed some points using R. Is that possibe since I can't seem to find the source file I used or commands for plotting maps and points? If not, sorry for the bother... Thanks, Jeff. Jeff D. Hamann Hamann, Donald & Associates, Inc. PO Box 1421 Corvallis, Oregon USA 97339-1421 Bus. 541-753-7333 Cell. 541-740-5988 jeff_hamann

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 mailing list archive at Nabble.com.

I'm having some trouble installing the add on packages. I would like to install the RMySQL_0.4-3.tar.gz add on package. I've downloaded it and unzipped the package. The INSTALL file says R CMD INSTALL [-l<alternate-dir>] RMySQL_<version>.tar.gz does this mean in the R system enter the following? > INSTALL -l /root/RMySQL_0.4-3.tar.gz (?) when I do this I get... Error:

----- Original Message ----- From: "Richard A. O'Keefe" <ok at cs.otago.ac.nz> To: <paul at datavore.com> Sent: Thursday, September 04, 2003 2:56 AM Subject: Re: [R] Overlaying graphs > I do not know how to overlay the curve graphic on top of hist graphic. > > Do you know about the "add=TRUE" option for plot()? > > I am hoping to show visually

# 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.

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

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 + .365) etc.

Hi, r-gurus: I happened to have a question in my work: I have a dataset, which has only one dimention, like 0.99037297527605 0.991179836732708 0.995635340631367 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 :)

I've developed a dll (win32) of a bunch of functions prototypes as: /* for building on *nix */ #ifndef WIN32 #define __stdcall /*nothing*/ #endif void __stdcall dbh_2_height( const unsigned long *func_idx, const unsigned long *metric, const unsigned long *species, const double *dbh, double *pred_height ); and have created a R interface