similar to: sub in boxplot doesn't do subtiles

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

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

Hi I was trying to install a package (sgeostat) under windows but the r-faq only tells how to do it under unix. Question: 1) is there documentation on how install a package under windows? Thanx bob sandefur -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

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

Hi list: I work with a lot of laboratory analytical data and I often have inconsistent names of files and variables within those files so I wrote this sourcearg function to facilitate handling file and variable names as both character and R names. The source of the function is given below sourcearg=function(arg){ su=sprintf("%s\n",arg) sink("TMPTMPRarg")

Hello! I have a problem with integrate() in my function nctspa(). Integrate produces an error message "evaluation of function gave a result of wrong length". I don't know what that means. Could anyone suggest me what is wrong with my function? These are the examples of function calls that work OK: nctspa(a=1:10,n=5) nctspa(a=1:10, n=5, mu=2, theta=3, renorm=0) This does not work:

HI >i have a problem please help me to solve it: http://r.789695.n4.nabble.com/file/n4639434/aj.pdf aj.pdf >i want to calculate the vecteur a[j] where j: 1...8 >this is the code in R: >aj.fun <- function(j, i, X, z, E, beta0, beta1){ + n <- length(X) + iX <- order(X) + iz <- order(z) + e1 <- -(beta)*z[ iz[1:(i - 1)] ] + numer <- E[j] - sum( X[ iX[1:(i - 1)] ]

Hi- There is a public doman geostatistics package used in the minieral industry called GSLIB and I would like to use it from R. However it is fortran (but portable and easy to compile on windows and linux)and has its own file format and control file format. It seems easiest to write an r function to dump r data to an ascii file, build a control file in ascii; call the GSLIB program with

Hi- I wanted percentage labels on a qqnorm x axis. I used the following: #make up some data jt<-rnorm(100) #qqnorm with percentile x labels qqnorm(jt,axes=F,xlab="Percentiles") box() l<-c(1,5,10,30,50,70,90,95,99) axis(1,at=qnorm(l/100),label=l) which was ok for this data set. Does anyone have a more general solution which would allow one to specify the number of x labels

Hello! I have a problem with integrate() in my function nctspa(). Integrate produces an error message "evaluation of function gave a result of wrong length". I don't know what that means. Could anyone suggest me what is wrong with my function? These are the examples of function calls that work OK: nctspa(a=1:10,n=5) nctspa(a=1:10, n=5, mu=2, theta=3, renorm=0) This does not work:

I have connected asterisk 1.0.7 with Avaya Definity via E1 with a TE100P Digium Card. Inbound calls are working perfectly and I dont have any problem. But when I try to make an outgoing call with my softphone (xlite) I am getting the following messages. Hungup 'Zap/13-1' Executing Dial("SIP/IZ-bc0a", "Zap/g1/3118") in new stack Called g1/3118 Channel 0/1, span 1 got

I want to add an outer subtitle to my 2x3-plot, but it's distance to the lowest plots is very high: 6 lines or more?! If I choose oma=c(5,0,2,0), it lies out of the plotting region and disapprears. Obviously I make something wrong here. Help appreciated, Thomas x=seq(from=1, to=3.5, length=100) par(mfrow = c(2, 3), oma=c(6,0,2,0)) ## oma=c(5,0,2,0) will not work?! plot(x,x^2) plot(x,sin(x))

Hi, I've read up on readBin() and chapter 6 in the R Data Import/Export manual, but I still can't read a binary file. Here is how the creator of the file described the code that would be needed in Fortran: "Every record has a return in fortran. The length of each record is nx*ny*4. To read you would use the following: nlayx = nx*ny*4 do iz=1,nz,4 read(binary file) var(1:nlayx)