similar to: Performance note: Preallocating helps? and two questions

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]

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

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

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 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 want to pick 3 "integer" random numbers from 1 to 10. How to do this? Thanks, Yu-Ling Wu __________________________________________________ Do You Yahoo!? Get email at your own domain with Yahoo! Mail. http://personal.mail.yahoo.com/ -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hello all, to obtain a list of all objects in all search paths, I've found the following to work: > biglist <- sapply(1:length(search()), objects) This more obvious one, however, does not work: > biglist <- sapply(search(), objects) Error in pos.to.env(pos) : invalid argument Still, search() gives [1] ".GlobalEnv" "package:ctest" "Autoloads"

Currently (R-1.4.1 as well as R-devel, according to http://stat.ethz.ch/R-alpha/R-devel/library/base/html/par.html), the html version of help(par) shows [...] lab A numerical vector of the form c(x, y, len) which modifies the way that axes are annotated. The values of x and y give the (approximate) number of tickmarks on the x and y axes and len specifies the label size. The default

Hello, has anybody written a function to plot density ellipses (95%, 99% or anything) in a scatterplot? I found nothing in any package, nor in the list archives. There does seem to be a contributed package "ellipse" for S-Plus (on S-Archive), but it does a lot more than what I would need. Still, if anybody ported it to R, I'd be grateful for a link. I'm a bit afraid to try the

Hello, I'd like to propose an extension to the function summary.aov. In Splus (2000, I don't know about other versions), summary.aov allows a parameter ssType to be set to 1 or 3 (defaults to 1) to choose the type of Sums of Squares. I know I can get Type III SS in R with drop1(model), but including the functionality into summary.aov would, in my opinion, - yield a more usable table

Hello, first, I'd like to congratulate the core team to the new R version 1.1. I think it's a great update, with glimpses into an even greater future (tcltk!). While playing around with the new functions (on Win 95), I found the following: As the documentation states, when I generate a plot and save it with recordPlot, I can regenerate it by printing the variable: > plot(1:10) >

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, does anyone know of a quick way to create a data frame "like" another, but with more rows? What I'd like to do is this: if mydata is a data.frame like a b c 1 TRUE yes 2 FALSE no 3 TRUE yes I'd like to get mydata2 with the same column names and column types, but without the values and with more rows. All I could think of was to manually do

How can I access a function created in one data directory in another? I have been sourcing the text file that contains the function each time I go to another data dir. This is cumbersome, and I suspect there is a way to make functions global. Thanks in advance, S. David White Dept of Speech and Hearing Science Ohio State University Columbus, Ohio

On Thu, Jun 22, 2000 at 05:01:32PM +0200, p.dalgaard@biostat.ku.dk wrote: > I think this is pretty certainly a bug, so I'm cc'ing this to r-bugs. > > -p > > Kaspar Pflugshaupt <pflugshaupt@geobot.umnw.ethz.ch> writes: > > As the documentation states, when I generate a plot and save it with > > recordPlot, I can regenerate it by printing the variable:

On Thu, Jun 22, 2000 at 05:01:32PM +0200, p.dalgaard@biostat.ku.dk wrote: > I think this is pretty certainly a bug, so I'm cc'ing this to r-bugs. > > -p > > Kaspar Pflugshaupt <pflugshaupt@geobot.umnw.ethz.ch> writes: > > As the documentation states, when I generate a plot and save it with > > recordPlot, I can regenerate it by printing the variable: