similar to: Error in vector("double", length) : vector size specified is too large....VLDs

Dear R users, I'm trying to write a small function in order to compute empirical cumulative density function.All seems to work but when I try to plot the function, I always get error messages. This is the function I use mecdf<-function(u,v,z) { u=sort(u) v=sort(v) n=length(u) nb=0 for (i in seq(1,n)) { if (u[i]<z & v[i]<z) { nb<-nb+1 } } nb=nb/n

>From the following: basin.map <- readAsciiGrid("c:/temp/area.asc", colname="area") I have a SpatialGridDataFrame which has the x and y cordinate of a cell, and the drainage area of that cell. There are many cells with a low drainage area (in my case, 33000 with an area of 37.16) and one cell with the highest drainage area (again, in my case, a drainage area of of

Hi, I've collected quite a bit of elevation data (LIDAR elevation points) and am looking for a suitable platform to do analysis and modeling on it. The data is sitting in an Oracle database, one table, 200 million rows of x,y, and z. I'm trying to figure out what hardware resources we need to reserve in order to run 64 BIT R for this size data. Here's my question: Is the 64 BIT

ComRades: I am geeting the error message Error in ldots[[which(MET)]] : attempt to select less than one element when I try to fit the geostatistical model with the likfit() function of geoR. I have tried with old data for which likfit() successfully maximised the likelihood in previous versions of geoR, and yet the current version fails. I have tried in Windows Vista and Windows XP (I haven't

Hello, I am attempting to get the standard deviation in multiple distance bins in my spatial data. It appears as though the 'variog' command in the geoR package will do the trick, as one of the outputs from 'variog' is 'variog$sd', which, according to the manual, is the "standard deviation of the values in each bin". However, when I run this command, the

Hi, I'm a new R user. My goal is to do a variogram using geoR. I started by trying to do the example in the geoR Illustrative Session using my own data. I am able to read in my Ascii data using: D <- matrix(scan("file.dat", n=530*3), 530,3, byrow=TRUE). Then I use: as.geodata(D, coords.col=1:2, data.col=3) to make the object D geodata. I check the descriptive statistics,

I'm using R for Mac OSX version 1.14 (2129) and the geoR package version 1.6-5 (the current version in the R repository). I'm running R in OS 10.4.6 on a Mac G4 iBook (933MHz, 640 MB DDR SDRAM). I searched the R archive and did not find a posting on this issue. I want to use the variog and variog4 functions of geoR to characterize the pattern of spatial autocorrelation of tree

dear list I am student M.S. statistics in department statistics . I am working in the function "nls" in the [R 2.3.1] with 246 data and want to fit the "exp" model to vectors( v and u ) but I have a problem to use it u 5.000000e-13 2.179057e+03 6.537171e+03 1.089529e+04 1.525340e+04 1.961151e+04 2.396963e+04 2.832774e+04 3.268586e+04 3.704397e+04 4.140209e+04

Dear GeoR-er, If I use the variog function in the latest release of geoR, the first lag is always ignored. For instance, if you read in geodata, calculates the variogram using the variog function and give in a uvec like uvec=seq(0,max,by=2.44), it only starts giving results from distance=4.88 and ignores 2.44! This wasn't the case in former versions of geoR. Is this done on purpose? Why?

I wonder what is the default value for the argument 'cutoff' when not specified in the variogram.formula function of gstat. Computing variogram envelops within gstat, I am comparing the results obtained with variog in geoR and variogram in gstat, and it took me a while before understanding that the cutoff default value is not the maximum distance. Can Edzer tell us about it? All the

Um texto embutido e sem conjunto de caracteres especificado associado... Nome: n?o dispon?vel Url: https://stat.ethz.ch/pipermail/r-help/attachments/20060628/2ff7714d/attachment.pl

Hey compiler peeps, Cray is ramping up a number of exciting projects and we're looking for new hires. Obviously parallelism has been our core focus, but the challenges of manycore and accelerator technology are presenting new twists requiring us to bring new solutions to bear. After the successful launch of the Jaguar machine (the fastest for open science,

I believe that actually means configure was able to use the BLAS you specified (MKL in this case). You can scan through config.log to make sure. Andy > From: > > Fedora Core 3 > > I installed the Intel MLk, > > and tried : --with-blas="-lmkl -lguide -lpthread" > > and got : External libraries: readline, BLAS(generic) > > thus I'm

Hello Everyone, We are using the LLVM infrastructure inside the RTSC, a operating system aware compiler used for research in the field of real-time systems. For further information please check out our website: https://www4.cs.fau.de/Research/RTSC/ In the future we would like to use this new logo for our project: https://www4.cs.fau.de/~klaus/rtsc-wyvern.png Since we directly modified the

Hi While performing migration of a guest virtual machine from a physical node to another, it would be nice if we can figure out the free pages in the guest VM and opt not to transfer them. 1.I came across in the paper "Xen and the Art of Virtualization", that guest VM can be made to give the free pages back to Xen before doing migration. 2.What is the way to figure out free pages in

Hi all, I tried to produce some kriged surfaces with geoR (latest version). The size of the grid should be around 900 x 650 cells (what I find is not a very big grid), and the number of points is around 2500. The command krige.conv stopped after arround 5 min saying it can not allocate a vector with around 1.5 billion units. Sounds reasonable. Is there a workaround? How would I partition the

Hi all, I use kriging to interpolate the precipitation from stations, but the map of this results show lots of stripes. (please see the attachment)I think there's something wrong with the setting of the dimension of this matrix, however, I have no idea how to know or test to see if this setting is correct or not.I've tried to switch the latitude and longitude, but still got the same

Hello again, since I am back from my holidays and there is still no reply to my post, I was wondering if there is no further interest in new projects using the LLVM framework or if it was just overlooked. I know that the copyright question is probably a tough one (therefore we couldn't decide it for ourselves), but we somehow expected to get some feedback or being listed on

I'm not sure I'm barking up the right tree here, but would I need to make use of groups to plot two separate datasets within ONE panel in xyplot? The desired end result is a single xy plot of two separate (but similar in values and ranges). Full code follows, xyplot code at bottom #########Determine Frequencies ##########coastal_slope #needs the maptools package to read ESRI grid

On 9/24/13 5:40 AM, Tobias Klaus wrote: > Hello again, > > since I am back from my holidays and there is still no reply to my > post, I was wondering if there is no further interest in new projects > using the LLVM framework or if it was just overlooked. It was just overlooked. If you don't get a response right away, feel free to ping as you've done here. > > I