similar to: To load Add-ons

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

The solve function in r-patched has been changed so that it applies a tolerance when using Lapack routines to calculate the inverse of a matrix or to solve a system of linear equations. A tolerance has always been used with the Linpack routines but not with the Lapack routines in versions 1.7.x and 1.8.0. (You can use the optional argument tol = 0 to override this check for computational

Folks, I asked a question on this mailing list about the subset support of lm(). In a flash, I got three helpful responses from Rajarshi Guha <rxg218 at psu.edu> <http://jijo.cjb.net> Erin Hodgess <hodgess at gator.uhd.edu> and Peter Dalgaard <p.dalgaard at biostat.ku.dk> :-) and it was just great. The mistake I was making was in not understanding the notion of a

Hallo, I'm using the R Version 2.3.0 (2006-04-24) on Suse Linux 10.1. With an older R and Linux version I could write a R-function into a file and execute it with the command: R --vanilla < script.r In the file script.r was code like this: postscript(file="results.eps") x<-2 y<-3 plot(x,y) dev.off() It worked fine, but now nothing happens. In the manual page this option

Hi, I am trying to do a bayesian prediction for soil pollution data above a certain threshold, using geoR. Everything is working fine until i am doing the krig.bayes. I tried to do the prediction on a grid 67 by 113 cells and my computer is freezing to death. At larger numbers of cells it tells me after a while that it reaches the max. memory of 511 Mb. My computer has only 512 Mb of RAM.

I've got a bunch of similar datasets, all of which I've fit to linear models. I'd like to easily create arrays of a specific parameter from each linear model (e.g., all of the intercepts in one array). I figured I'd put the model objects into an array, and then (somehow) I could easily create corresponding arrays of intercepts or residuals or whatever, but I can't the

Dear All, I learned about R during my studies at Cologne University of Applied Science. Now I work at Allianz Dresdner Bauspar AG and I would like to install R here too. Is there any license issues that need to be taken in consideration, any fees to pay by the company? I know there are answers to this in the FAQs but I didn't really understand the legal language. Could you therefore just give

Dear subscribers of R-devel I am experiencing that R crashes (further details are given below) in some complex matrix calculations when EISPACK=TRUE has been specified in eigen(). I discovered the behaviour some months ago just after the release of R-2.2.0, and it has been lying on my desk since. I apologise for not having nailed the problem down to a simple function call, but I thought I

Dear all, how can I plot a line graph and a bar graph in one single figure? I tried to combine "barplot" and "plot". Even though they both have the same x-values (1 to 55),  it just doesnt look as if they match in their scale (the barplot is much wider than the "plot"....even though I tried to put limits on the x-axis). Here is an example of what I did:

Dear R-helpers, I would like to do a Spearman rank order test, and used the cor() function with the method "spearman". It gives me a number (correlation coefficient?) , but how can I get the p-value? Thank you for the help in advance! Regards, Anne-Katrin -- [[alternative HTML version deleted]]

Dear All, I have data in HDF file format and would like to read it into R. I have tried the package hdf5 without success. Any ideas and suggestions?? Kind regards, Katrin -- Katrin Fleischer Vrije Universiteit Amsterdam Faculty of Earth and Life Sciences Subdepartment Hydrolgy and Geo-Environmental Sciences Room E-360 De Boelelaan 1085 1081 HV AMSTERDAM Tel: +31 20 59 87391

Dear Madams and Sirs, because the DOI foundation recommends a new, secure resolver [1], I wanted to suggest the attached patch. It a) updates a static DOI link in the docu, b) the code chunks that generate new ones, and c) a reg-ex that detects DOI links. Hopefully, my first venture into SVN & patch files was done correctly. My apologies, if there is anything amiss. In that case, I'll

Full_Name: Ole Christensen Version: 1.6.1 OS: linux-gnu Submission from: (NULL) (130.225.18.176) In package boot : > inv.logit(800) [1] NaN where it should have been 1. The problem is caused by exp(x) returning Inf when x is large. One way of fixing the problem [there may be better ways] would be to include the line out[x > 709] <- 1 in inv.logit() Cheers Ole >

Hello, I would like to label outliers (or all dots) in a plot plot(as.matrix(ValAddInd_byYear), as.matrix(Cons_Elec_Ind_byYear)) ideally by the 30 different row.names, but anything else (x or y values, for example) would already be helpful. I've tried various things with pch, but they didn't work and I can't find anything else. Any hint appreciated! Thanks, Katrin

Hello, I've downloaded the package "mass" and I've tried to install it, but it didn't work. Can you help me (I use Windows)? Katrin -------------- next part -------------- An HTML attachment was scrubbed... URL: https://stat.ethz.ch/pipermail/r-help/attachments/20020707/e560f42b/attachment.html

Hi From the documentation I can''t see if I can define the following schedule: I want a certain command to be run once a day, monday-friday on a certain time of the day. Can I define that somehow...? -- Med venlig hilsen Juri Rischel Jensen Fab:IT ApS Vesterbrogade 50 DK-1620 København Tlf: 70 202 407 / Fax: 33 313 640 www.fab-it.dk / juri@fab-it.dk

Dear all, a program that worked well for weeks today gave me consistently the error message no lines available in input referring to the lines for (i in (0:(timeintervals-1))) { j=subjectquantity+6+i*(subjectquantity+7) print (j) results<-read.table(file, header=F, skip=j, nrows=subjectquantity) timeintervals have been specified as argument as 12 subjectquantity as 9 the file it should

Dear All, I fear that this is a really easy question but I do seem to go around in circles.. I have 2 points on a plot and would like to calculate the slope of the line drawn through these 2 points. that cant be so hard?! Thank you in advance, Katrin -- Katrin Fleischer Vrije Universiteit Amsterdam Faculty of Earth and Life Sciences Subdepartment Hydrolgy and Geo-Environmental Sciences Room

Hi, I have run into this problem a couple of times now and hope you can help! If I want to plot mulitiple series with differing x-axis values (however, all in the same range) in 1 plot with one common axis R obstruses the plots. E.g. if I plot water content against time and I start with the sampler at 5 cm depth, it plots that one right. But the next depth, 15 cm, is measured at slightly

Dear List, How do I pick the p-value out of the urppTest result? For adfTest the p-value can be extracted by A2 at test$p.value following A2 <- adfTest(myData[,i], lags=2, type=c("c")) What do I do for urppTest? The above doesn't seem to work. There is a slot @test with $output, which is a list of various test results that didn't want to give away only the p-value