similar to: percentile labels in qqnorm

hi- I have been using postscript()and contour() to do quick plots of geological sufaces and it would be handy to be able to dump a plot as an autocad dxf file for vector import into autocad. A search of rhelp of 99 thru 01 gave no matches for dxf. Any suggestions for a dxf or dump of a plot as vectors? Thanx Thanx Robert (Bob) L Sandefur Principal Geostatistician Pincock Allen & Holt

hi Suppose I have a mixture of 2 distributions generated by rtwonormals <- function(npnt,m1,s1,m2,s2,p2){ rv<-vector(npnt,mode="numeric") for( i in seq(1:npnt)){ if(runif(1,0,1)<=p2){ rv[i]<-rnorm(1,m2,s2) } else{ rv[i]<-rnorm(1,m1,s1) } } return(rv) } x <- rtwonormals(50000,0,100,500,500,0.05) #and I try to fit these with (based on thread: [R]

hi listees- Running windows 2000 and r 0.99 or 1.00; I get different lines on screen and print out: > j<-read.table("jnk.rin",head=T,sep="\t") > attach(j) > #gives lines I want > qqnorm(Dog,pch="o",main="Cumulative Frequency Log Dog Size",log="y",sub="Figure 2",x =c(-3,3),ylim=c(2000,1000000)) >

hi- in r 1.1 on windows 2000 with length(AU) of 35833 AUcap30<-0 for(i in 1:length(AU))AUcap30[i]<-min(30,AU[i]) took over an hour on pentium II 300 mhertz (I esc'ed before it finished) but AUcap30<-AU for(i in 1:length(AU))AUcap30[i]<-min(30,AU[i]) is very quick (a few seconds) Is this performance difference common in r (ie is linux the same way)? Are there other tricks

hi- in r 1.1 on windows 2000 with length(AU) of 35833 AUcap30<-0 for(i in 1:length(AU))AUcap30[i]<-min(30,AU[i]) took over an hour on pentium II 300 mhertz (I esc'ed before it finished) but AUcap30<-AU for(i in 1:length(AU))AUcap30[i]<-min(30,AU[i]) is very quick (a few seconds) Is this performance difference common in r (ie is linux the same way)? Are there other tricks

hi- If I source the following for(k in seq(1:20)){ x<-runif(20000,min=-500,max=2000) y<-runif(20000,min=-500,max=2500) z<-runif(20000,min=-10,max=10) cat(k,"file",memory.size()) cc<-rainbow(11) plot(x,y,asp=1i,xlim=c(-500,2000),ylim=c(-500,2500),main=k,cex=1.0) for(i in seq(-10,10,2)){ points(x[z > i],y[z > i],col=cc[(12+i)/2],cex=1.0) } rm(x,y,z)

hi- I downloaded r 0.90.0 for windows and ran rwinst.exe and the base package would not install for me. I had to unzip rw0900b1.zip and r20900b2.zip manually and then do the rest of install with rwinst. Is this due to b being replaced by b1 and b2? Robert (Bob) L Sandefur Principal Geostatistician PAH rls at pincock.com (primary) rlsandefur at aol.com (secondary) 303 914 4467 v 393 987-8907 f

Hi- I have about 20 groups for which I know the mean, variance, and number of points per group. Is here an R function where I can plot (3 group example) something like | |----- 2 | |-----------| | | |----- 1 | | |------| | |----------------- 3 | | 0 25 50 75 100 ie 1 and 2 are different at 75% level of confidence 1 2 combined are

Dear all, Will R run under Win2000 or is it safer to stay with NT4.0 for a while? Regards -- Dr. Jens Oehlschl?gel Analyse BBDO InterOne Gr?nstr. 15 40212 D?sseldorf Tel.: +49 (0)211 1379-187 Fax.: +49 (0)211 1379-461 http://www.bbdo-interone.de -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

R list- Under Windows 2000 and r 99: I wanted a subtitle in boxplot and tried: > boxplot((AU[ROCK==2 & AU>-1]+0.1~iz),log="y",main="Boxplot of AU by 25 M elevations",sub="aa") And I got: Error in sort(x) : only vectors can be sorted In addition: Warning messages: 1: is.na() applied to non-(list or vector) in: is.na(x) 2: is.na() applied to non-(list

hi- I would like to get a value set with the dos set command (e.g. set dog=cow) in a rterm session on rterm 1011 under windows 2000. if I try in a dos window: set dog=cow echo %dog% i get cow (what I want) then i start rterm system("echo %dog%",intern=T) gives "%dog%" I can of course do something like set dog=cow echo %dog% >jnk.rin and then start rterm and

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

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

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

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

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

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 )

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,