similar to: references for 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

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

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

# 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 R People: Here is an interesting question(I think) Suppose I want to calculate Centered Moving Averages; i.e. x[1] <- ( sum(y[1:12]) )/12 x[2] <- ( sum(y[2:13]) )/12 and so on. Of course, this is easily done through loops. However, I have been trying to do this more elegantly, but have failed. I have tried things like j1 <- 1:109 j2 <- 12:120 x[1:109] <- ( sum( y[j1:j2])

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|

Dear R People: Is there a function for exponential smoothing, please? R version 1.4.1 for Windows. Thanks in advance! Have a great weekend! Sincerely, Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown 1 Main Street Houston, TX 77002 mailto: hodgess at uhddx01.dt.uh.edu

Dear R People: I have the source for R version 1.1.1 on a UNIX system alpha-dec-osf4. The ./configure works fine. When I do the make, it works fine until it tries to build the package tcklc. I get an error of Tk_Saf_Init not identified Then it kicks out. Any suggestions, please? Thanks! Sincerely, Erin M. Hodgess, Ph.D. Associate Professor Department of Computer and Mathematical Sciences

Dear R people is there a module for quality control, please? X bar charts, and such? perhaps a contributed package? Thanks in advance! Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown 1 Main Street Houston, TX 77002 mailto: hodgess at uhddx01.dt.uh.edu

Dear R People: Do any of the linear model or regression function calculate the Durbin-Watson test statistic for serial correlation, please? I found the hat matrix, studentized residuals, and so on, but no D-W. Thanks so much! Sincerely, Erin M. Hodgess, Ph.D. Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown One Main Street Houston, TX 77002

Dear R People: Here is a very simple question, please: I have a data frame with the state name, state division, and population. Note: state division is something like Atlantic, Mountain, etc. I want to produce a table with the state division and the total of the populations for the states within that division. I've been messing with "table" but so far, no luck. This is R1.3.1,

Hello R People: Suppose I create a function, x.e, via the fix command: > x.e <- fix(x.e) When I type in x.e, the function code appears, along with the name/number of the environment. However, if I create a function, y.e, and source it in, when I type in y.e, the function code only displays(no environment). Why is this, please? (It's not a bug....I just wondered). I am using

Dear R People: Does anyone have any code for Fractals, chaos, or anything like that, please? This is strictly for demo purposes...decorative only. This is R version 1.5.1 for Windows. Thank you in advance! Have a great weekend! Sincerely, Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown 1 Main Street Houston, TX 77002

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

Hello, Although your message is not very explicit about the structure of your data frame, I think you can get what you want by using a combination of the function "split" and the "apply" functions family: apply, lapply and sapply. If you provide a more detailed organization of your data frame, perhaps we can give to you a more wiser advice. Thanks, Carlos. P.S: By the way,

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)

Hello, R People: How do I change the size of a graph, please? I tried par(cin=c(8.5,11)) and it said that Warning message: parameter "cin" can't be set in: par(args) > par(fin=c(8.5,11)) > plot(xa,ya) Error in plot.new() : Figure region too large Also, can I set for portrait or landscape, please? Thanks in advance! Erin M. Hodgess, Ph.D. Associate Professor Department

Dear R People: I'm trying to use the predict command on an arima0 object. I do the following: xm.arma <- arima0(xm2,order=c(1,0,1)) predict(xm.arma,n.ahead=2) and I get the message: Error in round(x, digits) : Non-numeric argument to mathematical function Any ideas what the problem might be, please? R version 1 4 1 on Windows. Thanks in advance! Sincerely, Erin Hodgess Associate

Dear R People: Here is a non-R probability question, please: There is a problem, originally posed from an American game show, called "Let's Make a Deal." A contestant makes a selection of one door, from Doors 1 through 3. He or she knows that there are 2 goats and 1 car behind those doors. Once the contestant makes the selection, say Door 2, then the host, Monty Hall, tells him

I recently wrote a function to compute functions over moving windows. (See below). It varies from your request by using the last N points rather than the last N/2 and next N/2 points, so you may have to modify it to get what you want. It will be in the next gregmisc package, which will go out in a month or so when the publication review process here finishes. (I hate red tape!) -Greg