similar to: difficult R-problem

Dear R-community The following loops produce the error message: Error in round(x, digits) : Non-numeric argument to mathematical function after performing the outer loop 6 times Splus can execute these loops but is very much slower than R . Interestingly, if the scale-argument in the step-function is omitted R performes these loops a few times more. What could be the reason for that error? Any

> install.packages("gam") Warning in install.packages("gam") : argument 'lib' is missing: using /usr/local/lib/R/site-library --- Please select a CRAN mirror for use in this session --- Loading Tcl/Tk interface ... done trying URL 'http://cran.cnr.Berkeley.edu/src/contrib/gam_0.97.tar.gz' Content type 'application/x-gzip' length 89613 bytes opened URL

Hola buenas noches, tengo un problema bastante gordo. ¿A alguno le ha dejado de funcionar las funciones pamer.fnc y mcp.fnc con la nueva versión de R? La semana pasada formatee el ordenador y ahora scripts antiguos no funcionan. La cuestión es que me precupa que no funcione el ejemplo de tutorial del autor. Os dejo un script que debería de funcionar y no lo hace

Hi all, My question is not really urgent. I can write a loop and solve the problem. But I know that I'll be in a similar situation many more times so it would be useful to find out the answer Is there a fast way to perform linear fit to all the columns of a matrix? (or in the one dimension of a multi-dimensional array.) I'm talking about many single linear fits, not about a multiple fit.

Is there a way to add residuals to plots produced by plot.gam in the mgcv package? I'm looking for something like what you get using resid=T in Splus plot.gam. Thanks in advance Toby -----Original Message----- From: Simon Wood [mailto:snw at mcs.st-and.ac.uk] Sent: 23 January, 2002 8:14 PM To: Toby.Patterson at csiro.au Cc: r-help at stat.math.ethz.ch Subject: Re: [R] multiple surfaces in

Sorry for re-posting the question, I did not get any reply. Kindly reply please if any one can. ### Hello R users, I have a montly time series over a several year period. It's easy to compute a monthly climatology (12 values), Now, I would like to calculate the corresponding standard deviation, ie the 12 values calculated from the january values, february values, etc. What's the best way

I have an irregular time series, which cannot easily be converted to a ts object (long stretches with NA). I have calculated the climatological mean and st.dev, but need help on how to calculate the anomalies from the climatology and my original data set. Below is an example, where I have indicated the basic idea of what I want to do under Question 1 and 2. I assume I will have to create a custom

University of California Merced Merced, CA Programmer Analyst II/III (Research Assistant) Job Code SSNRI723A Open until filled. ? In the Sierra Nevada Research Institute at UC Merced, act in support of research in applied climatology and statistical modeling for wildfire, energy and water resource management applications and assist the Principle Investigator with the development of software,

Hello R users, I have a montly time series over a several year period. It's easy to compute a monthly climatology (12 values), Now, I would like to calculate the corresponding standard deviation, ie the 12 values calculated from the january values, february values, etc. What's the best way for such a calculation ? Regards, Yogesh -- Yogesh K. Tiwari (Dr.rer.nat), Scientist, Indian

Dear R list members, I am seeking for an elegant way for matching NA elements of a matrix. Everthing I tried, the result was a vector. Here ist an example with vectors that works: a and b are vectors of same lenght representing pairs of data (a[i] b[i]) with some NA elements: ## R Code a <- (1:5) a[2] <- NA b <- (6:10) b[3] <- NA a 1 NA 3 4 5 b 6 7 NA 9 10 ## /R

R-users E-mail: r-help@r-project.org I have a quenstion on "gam()" in "gam" package. The help of gam() says: 'gam' uses the _backfitting algorithm_ to combine different smoothing or fitting methods. On the other hand, lm.wfit(), which is a routine of gam.fit() contains: z <- .Fortran("dqrls", qr = x * wts, n = n, p = p, y = y *

Dear all I use the library(netCDF) to read in NCEP data. The file I want to read has size 113 Mb. When i try to read it I get the following message: Error: cannot allocate vector of size 221080 Kb In addition: Warning message: Reached total allocation of 255Mb: see help(memory.size) I get a similar message when I try to read a file with 256 Mb in a PC with 2 GigaByte RAM. Is there something

Dear R-users I would like to draw a barplot with a special y-axis. Most of my data points are in the range from -50 to 50, apart from one value that is 550. I would now like to have a y-axis from -50 to 50, then a "break" and the "rest" of the y-axis, let's say from 500 to 600. In order to make clear what I'm trying to do, I attached a figure. Is it possible to do

Dear all, I have in my possession a netcdf from which I want to extract some data files. I have used the "ncdf" package to read the netcdf file and used the "get.var.ncdf" function to identify the variable i wish to use. The data is in the form of a time-series of geographical data points that relate to climatology variables. As such I have a large data frame,

I have a data set something like this: "YYYY", "Value" 1972 , 117 1984 , 73 1969 , 92 1976 , 113 1999 , 80 1996 , 78 1976 , 98 1984 , 106 1976 , 99 it could be created with: > dafSamp <- data.frame(cbind(c(1972,1984,1969,1976,1999,1996,1976,1984,1976),c(117,73,92,113,80,78,98,106,99))) The real dataset is of cause much larger, app. 100.000 samples

Good morning, suppose the following: m <- round(matrix(rnorm(16), ncol=4), 3) a <- rev(c(0.01, 0.025, 0.05, 0.1)) rownames(m) <- a colnames(m) <- c("0.25,0.75", "0.4,0.6", "0.1,0.9", "0.4,0.9") m 0.25,0.75 0.4,0.6 0.1,0.9 0.4,0.9 0.1 1.034 -0.119 -1.213 0.619 0.05 0.035 1.074 0.525 1.671 0.025 -1.687 0.960

Hi, I'm estimating a loglogistic aft (accelerated failure time) model, just a simple plain vanilla one (without time dependent covariates), I'm comparing the results that I obtain between aftreg (eha package) and survreg(surv package). If I don't use any covariate the results are identical , if I add covariates all the coefficients are the same until a precision of 10^4 or 10^-5 except

Hello, are there people outside who use R for analysis of wind-measurement data (meteorological or for planning of wind power stations)? Are there already scripts/modules available for analysing and displaying/plotting wind data in the way it is done in projection/planning of wind power stations? If not, would it be of interest to use R for this, and therefore adapt data-logger output to R (by

The Climate Corporation is a Google Ventures funded big data startup with a mission to help people adapt to climate change. We make sense of huge amounts of complicated data and run simulations <http://goo.gl/9Fv0J> to determine the risk of adverse weather, which influences 25% of GDP. We are looking for a Web Application Developer who are passionate about innovation and share the value

Hi all I have the next data frame year STOD SLAGNR TAL TALT TALVEKT 1 2002 2120006 57 1 NA 1 2 1997 97030032 57 NA NA NA 3 1997 97030071 27 1 NA NA 4 1997 97030005 57 1 NA NA 5 1997 97020127 60 NA 1 NA 6 2001 1160025 27 1 NA 1 7 1998 98020069 60 1 NA NA 8