similar to: implementing a foreign system with system?

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 R-Users Is there any package in R that work Indicator Kriging? Or somebody is working in a link between GSLIB library in FORTRAN for R? Thank you very much for your help. --

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

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

The NOAA National Centers for Coastal Ocean Science is hiring a Computational Ecologist, a statistical/computational ecologist with experience fitting advanced spatial models to marine wildlife survey data (e.g., seabirds and marine mammal transects, fisheries trawl surveys) in R and other statistical languages. This is a full-time, long-term stable contract position. We are looking for an

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

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, all, Have anyone used kringing included in R? How is it? Does it handle anisotropy data well? How does it compare with Kriging in Arc/Info? or other geostatistics software customized to do kriging or other geostatistics functions? I tried Easykriging, a geostatistics tool developed for Matlab. It has very nice GUI, but it does not provide library which i can call in my programs. so it is

hello I need basic explanation of variogram and how to use in R. I will appreciate if you supply them. regards Ahmet Temiz ______________________________________ XamimeLT - installed on mailserver for domain @deprem.gov.tr Queries to: postmaster at deprem.gov.tr ______________________________________ The views and opinions expressed in this e-mail message are ...{{dropped}}

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

Dear Users, I use Sgeostat to calculate co variogram (function est.variogram with two variables). If i inverse the 2 variables in the functions, the co variogram change! Normaly it is symmetric, isn't it? Someone could me explain why, or give me an idea to calculate co variogram in a easy way... Thanks in advance Laurent ps : sorry for my english Laurent William Dubroca Centre de

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

List: I have 2 related questions: (1) first I have x-y-z data, where x & y are the geographic locations of point values, z. I need to perform a normal score transform on the z-values and maintain their geographic location. So, how do I go from columns x-y-z to x-y-z-t (or x-y-t), where the t-values are the normal score transforms of the z-values? Can I use qnorm(ppoints(data)) to do

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 to everybody ! I don't have a question directly connected with R, but it takes starting point from the use of the R function Krig for kriging interpolation and from the idea of making a different implemetation of the function for my purposes. Doas anyone know where can I find free C or C++ codes of kriging ? As always, thanks a lot in advance, Paola.