similar to: matrix matching NA

Dear Everyone, I am using R 2.6.2 on my Redhat AS4. I installed the RPM offered on the cran website. I would like to use R from my bash console with the following command: /path/to/R --vanilla --slave --args < /data/myscript.R this script invokes a range of functions and tools, among others to process a bunch of raster images in a loop using GDAL. The thing is that releasing this command in

Some attribute x from 17 individuals was recorded repeatedly on 6 time points using a Likert scale with 7 distractors. Which statistical test(s) can I apply to check whether the changes along the 6 time points were significant? set.seed( 123 ) x <- matrix( sample( 1:7, 17*6, repl=T ), nrow = 17, byrow = TRUE, dimnames = list(1:17, paste( 'T', 1:6, sep='' )) ) I found

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.

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

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

Greetings, I have two questions that I could not answer from the documentation. A - ecdf and confidence intervals : Is there a (simple) way to generate confidence intervals (95%) for a ecdf? B - cross-validation of rpart trees : a colleague is using S to generate decision tree and mentioned to me the use of cross-validation. Is this function enabled in R ? if so, how should one proceed to

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,

Hi there In the course of my diploma thesis in climatology I have encountered a difficult R-Problem that I cannot solve. I want to fill R-Objects (whose names should depend on j) with numbers at the i-th position. The resulting Objects should be something like: RQuadratStep1, RQuadratStep2, RQuadratStep3 ... filled with Elements like c(0.324, 0.456, 0.657 ...) Below is a short version of

Dear List, after updating the exactRanksumTests package I receive a warning that the package is not developed any further and that one should consider the coin package. I don't find the signed rank test in the coin package, only the Wilcoxon Mann Whitney U-Test. I only found a signed rank test in the stats package (wilcox.test) which is able to calculate the exact pvalues but unfortunately

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

Hi, I was trying to draw an empirical distribution function with uniform confidence bands. So I tried to find a way to calculate values of the Kolmogorov-Smirnov Distribution but failed. I guess it must be hidden somewhere (since the ks-test is implemented), but I was unable to find it. Is there any way to do this? Thanks Leif Boysen

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,

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

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

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

Dear R users, I'm coming back to R after while. I have a data frame with 200 columns, each column has a name. How to extract all columns to a new dataset, but the specified (by names) ones? I was playing with that for a little bit using the vector syntax but got several syntax errors. Thanks, Rado

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, I'd like to draw a polygon graph. I used the package lattice extra which includes a function for that (see: http://latticeextra.r-forge.r-project.org/#panel.xyarea&theme=default). But i don't want the polygon ending with its filled border at the x-axis (like in my code). Instead it should be rotated about 90 degrees and end with its filled border at the y-axis. I'm pleased