search for: 2e3

Dear all, I think this is an easy task, but I don't know how to do it. Specifically, I have 69 columns with 300.000 rows. In each cell there is a code like "2E3", "4RR", etc. I now have a list that replaces this with values, e.g., old new 2E3 5 4RR 3 etc. The list ist about 1600 rows long, so also to extensive for normal solutions Do you how to do that? It would be great if you could help me, best, tobebryant

....data.i" not found' + for(k in 1:ncol(auto.data.i)) { + auto.data.i[j,k] <- 'some value' + }}} Thanks a bunch for all your help. Best, Michael Michael Drescher Ontario Forest Research Institute Ontario Ministry of Natural Resources 1235 Queen St East Sault Ste Marie, ON, P6A 2E3 Tel: (705) 946-7406 Fax: (705) 946-2030

...how I could use, e.g., the 'charmatch' or 'parse' commands to achieve this. Thanks a lot in advance for your help. Best, Michael Michael Drescher Ontario Forest Research Institute Ontario Ministry of Natural Resources 1235 Queen St East Sault Ste Marie, ON, P6A 2E3 Tel: (705) 946-7406 Fax: (705) 946-2030 [[alternative HTML version deleted]]

I have an R function defined as: f<-function(x){ return(dchisq(x,9,77)*((13.5/x)^5)*exp(-13.5/x)) } Numerically integrating this function, I observed a couple of things: A) Integrating the function f over the entire positive real line gives an error: > integrate(f,0,Inf) Error in integrate(f, 0, Inf) : the integral is probably divergent B) Increasing the interval of integration

...give you (after a long while) a bunch of libatlas*.deb files that you can easily install with sudo dpkg -i | | Now if you simply configure R --with-blas --with-lapack it will work, and you'll get a nice and fast BLAS, but unfortunately it will be single-threaded: | | > m=matrix(rnorm(4e6),2e3) | > system.time(tcrossprod(m)) | user system elapsed | 0.684 0.012 0.697 | | The problem is that all those nice libraries that are handled by update-alternatives are just the single-threaded part of ATLAS. If you want to use the parallelized part, you'll need to add the libpt*...

.... However, I'm not a "CG" expert and can live with a suboptimal result. --- optim.c (Revision 37878) +++ optim.c (Arbeitskopie) @@ -946,6 +946,8 @@ G3 = G1 / G2; else G3 = 1.0; + if (G3 > 2e3) + G3 = 2e3; gradproj = 0.0; for (i = 0; i < n; i++) { t[i] = t[i] * G3 - g[i]; Andreas -- Andreas Westfeld, 0432 01CC F511 9E2B 0B57 5993 0B22 98F8 4AD8 EEEA <westfeld at inf.tu-dresden.de> http://www...

...nge on my Linux (Ubuntu 7.10) machine, with occasional segfaults. Today something else and reproducible happened: If I type the code below (meant for calibrating data), I get the error message that "the C stack usage is too close to the limit". calcurve <- cbind(1:2e4, 1:2e4, 1:2e3); #dummy curve, real one is more complex caldist <- function(cage=Cage, error=Error, sdev=Sdev, times=Times, By=By) { theta <- seq(min(calcurve[,1]), max(calcurve[,1]), by=By); interpolate <- function(th, col) { if(th==calcurve[1,1]) {calcurve[1,col]}else if(th==cal...

...ound this same error message mentioned in another context (levelplot) as indicating a bug (original bug report PR# 6005 on Mon, 22 Dec 2003) ------------------------ Michael Drescher Ontario Forest Research Institute Ontario Ministry of Natural Resources 1235 Queen St East Sault Ste Marie, ON, P6A 2E3 Tel: (705) 946-7406 Fax: (705) 946-2030 ------------------------ Code: > sfalls.plot.comp <- matrix(nrow=plot.count, ncol=spec.count, byrow=T) > colnames(sfalls.plot.comp) <- levels(SPECIES) ### SPECIES, SPP_VOL, & PLOT are columns/variables in a previously read data file > sfa...

...I am puzzled by this error message and can't find anything in the archives etc. Your input would be greatly appreciated. Best regards, Michael Michael Drescher Ontario Forest Research Institute Ontario Ministry of Natural Resources 1235 Queen St East Sault Ste Marie, ON, P6A 2E3 Tel: (705) 946-7406 Fax: (705) 946-2030 [[alternative HTML version deleted]]

...r example: q<-sample(1e7,2e6); r<-sample(1e7,1e4) system.time(intersect(q,r)) => 70 mS qb<-as.bit.which(q,1e7); qr<- as.bit.which(r,1e7); system.time(qb&qr) => 3 mS system.time(set_intersection(as.set(q),as.set(r))) => many minutes (with |q|=2e3, time is 3 S) I was wondering if there were other set packages which use different representations of sets (e.g. sorted lists, heaps, etc.) to make some set operations more efficient. In particular, I'd be interested in an intersection operator which is fast and space-efficient when the sets a...

...#39;t R use KB, MB and GB when reporting on sizes kilobytes, megabytes and gigabytes? More specifically, format() for object_size objects (returned by object.size()) uses Kb, Mb and Gb, which AFAIK is more often used to indicate sizes in kilobits, megabits and gigabits, e.g. > format(structure(2e3, class="object_size"), units="auto") [1] "2 Kb" > format(structure(2e6, class="object_size"), units="auto") [1] "1.9 Mb" > format(structure(2e9, class="object_size"), units="auto") [1] "1.9 Gb" > form...

Hi, I was trying to evaluate the pearson correlation and the p-values for an nxm matrix, where each row represents a vector. One way to do it would be to iterate through each row, and find its correlation value( and the p-value) with respect to the other rows. Is there some function by which I can use the matrix as input? Ideally, the output would be an nxn matrix, containing the p-values

...ere: fisher exact test p-value) iteratively using 4 matrices with identical dimensions. And as a result I aim for an array with identical dimensions containing the corresponding p-values. Please consider some code using a trivial example with 3x4 arrays below. Eventually I would like to run code on 2e3 x 7e6 arrays, for which someone suggested Amazon EC2 already... Q1: would you agree that fisher.test() is not vectorizable? e.g. fisher.test( matrix(c(Ax,Ay,Bx,By),ncol=2) ) does not work Q2: direct use of Ax, Ay, Bx, By as input instead of a (list) transform for the input would seem beneficial fo...

I am trying to share a CD-RW from the samba server (SuSE Linux 7.3, kernel 2.4.10) using Samba that came with the distribution. I am using user level sharing and have successfully shared out disk shares with both r and rw access. I can also access the CD-ROM and CD-RW in read mode and even do an installation from the CD-ROM to a client. When I try to create a CD using CD-RW, I get the following

Hi, More than 1000''s of classes and filters gets created/deleted in run time. Are we doing it correctly( as our requirement is such). Since few days we are facing lots of problems, like the server gets hanged, we get some junk messages in our logs. we have to reboot our server every now and then. What is the max limit of classes/filters can be created. we are using redhat 7.3.

...36+16u78Yf?=k)S1O5Ch-Oeq#{^fzoqt?8&e5%WUvG>k}UHkX#fxmVJ;$*8; zS0}158}u&2N|Cg<<lC(VQ-y1a*=!_rCRZprWo9esTc}AjS6R#^(qvB~ZDm%21xmWN zW{{Drg|cfg>guTaZm<Gb154R_i*!bV9_wWuX{$Ec(-GJm=Ig5M<|>@wrAcPfRDGer zWS5s4b+$^n88NOV8B}d5L#2E*=~bnqUQg-=nN9XQvch1qTkG7Z8hQ#7KPQY*VL6?_ zL|R?4t6^2E3}R5FGwCWwt4*Pyv_Fe_cGnCtmsQ*7Cpjz#jflD0T1MJDHF<QB&Q)mg z>uh$iihh(*7*vqvDz`k?tIkwmBq{o2nf291a*(bJQmhWq$7ZiHk~W{z{QP{BN+ZjR zI?8~#<Xfx;JLv{Hvyf0LF|uWF09LaRqJUaM!y#+!$aJ_njnjjzlxj3UlF=BU!X5?i z4>S~H774jqW^>8X!QWbOju}$hO9LCULhKARB-N5Ey}TU$V8T=Q!d7kKeDotl{KeTe zR}EE8UZ_Fs4>DUJBSqD)N...