similar to: Retain parts of a matrix

Hello, I'm trying to develop a box plot of time series data to look at the range in the data values over the entire period of record. My data initially starts out as a list of hourly data, and then I've been using this code to make this data into the final ts array. # Read in the station list stn.list <- read.csv("/home/kbennett/fews/stnlist3", as.is=T, header=F) # Read in

Hi, I'm installing rgdal but I keep having failures because I have not been able to find a good source of information for the correct configuration settings when installing GDAL. My error from the R install.packages("rgdal") is below. Can someone point me to a good source to tell me how to set after ./configure when installing GDAL? I'd like to be able to work with raster

I'd like to sum a matrix only where the matrix meets a specific condition. The matrix has 365 rows and about 50,000 columns. str(cdem.mat.yr) num [1:365, 1:41772] -43.5 -48.4 -45.9 -38.4 -32 ... I'm having trouble replicating this because it is not working out for me, so I'm unsure I can provide an solid working example (apologies). I would like to subset my matrix where the values

Hello, I am wondering if someone can help me. I have the following function that I derived using nls() SSlogis. I would like to find its derivative. I thought I had done this using deriv(), but for some reason this isn't working out for me. Here is the function: asym <- 84.951 xmid <- 66.90742 scal <- -6.3 x.seq <- seq(1, 153,, 153) nls.fn <- asym/((1+exp((xmid-x.seq)/scal)))

Hi R-help, I am working with US COOP network station data and the files are concatenated in single rows for all years, but I need to pull these apart into rows for each day. To do this, I need to extract part of each row such as station id, year, mo, and repeat this against other variables in the row (days). My problem is that there are repeated values for each day, and the files are fixed width

I have a zooreg object and I want to be able to generate a value for seasons and 11-day composites paste it onto my zoo data frame, along with year, month and days. Right now I have the following to work from: eg. dat.zoo.mdy <- with(month.day.year(time(dat.zoo)), cbind(dat.zoo, year, month, day, quarter = (month - 1) %/% 3 + 1, dow = as.numeric(format(time(dat.zoo), "%w")))) For

Hello everyone, I have a dataset in the following format: col1 col2 # # # # # # # # # # # # What I want to do is: loop a random sample 10 times, and for each time it is sampled I want to run a correlation between both columns. What I have so far is this: >feb <- read.csv("corr.csv") # where the dataset is for February >attach(feb) >for(i in

Dear all, I need to do the same procedure on several files. But I don't know how to refer to the file name. Here is an example of what I am trying to do. List of files: file1(A,B,C, D1...Dn), file2(A,B,C,E1,...,En), file3(A,B,C,F1,...,Fn) Procedure I want to apply on each file: dft <- melt(df,id=c('A','B','C')) dft$X <- substr(dft$variable,1,3) dft$Y <-

Could any of you provide me information on a good Polycom phone dealers to utilize. One who provides firmwares ..etc Thank you! Kenny ______________________________________________________ Click here to donate to the Hurricane Katrina relief effort. http://store.yahoo.com/redcross-donate3/

Dear mailinglist members, I was wondering if there was a way to re-evaluate the instances of a tree (in the forest) again after I have manually changed a splitpoint (or split variable) of a decision node. Here's an illustration: library("randomForest") forest.rf <- randomForest(formula = Species ~ ., data = iris, do.trace = TRUE, ntree = 3, mtry = 2, norm.votes = FALSE) # I am

Hi, I am hoping some one can help me. I am learning to use C and would like to learn how to call c code in R. I have look at Writing R Extensions and I tried to copy the example on page 38 void convolve(double *a, int *na, double *b, int *nb, double *ab) { int i, j, nab = *na + *nb - 1; for(i = 0; i < nab; i++) ab[i] = 0.0; for(i = 0; i < *na; i++) for(j = 0; j <

Well, I'm a recent convert from WBEL. My biggest concern with CentOS is that the community here seems to want to be more than a recompile of RHEL. But WBEL is floundering, what with Katrina and Rita, and there really being only 1 developer behind it, etc. I offer an automated shell script to switch from WBEL4 to CentOS4 (easy, it's hosted on my home DSL line!) It assumes that

I am trying to find a definition for the two-tailed exact binomial test but have been unsuccessful. Can you help? --------------------------------- [[alternative HTML version deleted]]

Greetings R Users! I am posting to inquire about the proximity matrix in the randomForest R-package. I am having difficulty pushing very large data through the algorithm and it appears to hang on the building of the prox matrix. I have read on Dr. Breiman's website that in the original code a choice can be made between using an N x N matrix OR to increase the ability to compute large

Okay, here is the background. I have a PRI with 15 active channels on it. I originally setup all of them in group=1 and all outgoing and incoming calls used this group. The phone number that I have associated with these channels ends with 750 and that is how I direct the calls. i.e. In my extensions.conf I have: exten => 750,1,Dial(SIP/120,20) All this works fine. Now I have the need

Hi, I installed Wine07252005, my printer is EPSON LX 300X, it's CUPS, the name is EPSON9Pin. I tested it with a simple Delphi program and have two problems: 1) If i print using a report builder like RBuilder, it prints quite well, but the font does look good in the print out. 2) If I use spooler to print, like following: if StartDocPrinter(PrnHandle, 1, @DocInfo) = 0 then begin

I apologize for the double post. I am curious as to what the usefullness is of the multiple line appearance feature on Polycom phones. I setup our phones to register one line per extension but I hear the IP501's can do three line appearances. Why and how could this feature be applied? Thanks again all. Kenny ______________________________________________________ Click here to

I have a lot (more than one hundred) of files with tables of the same kind (quadratic, same size) and I want to obtain some statistics for every position on them. Therefore, as a result I want another table. I import every table, and create an object “read.table” for it, then I have try to create a “list” or a “data frame“ and directly utilize some functions like sd( ) without success, because it

Hello R-helpers. We have a question about AIC values. We ran several different GLMs (quadratic, interactions) and our lowest AIC values were the same for several different models. We don't know how to interpret this. Any thoughts? Thanks in advance. Katrina Shelding _________________________________________________________________ e. [[alternative HTML version deleted]]

hi how does one control traffic of ftp using tc.I tried few things ,dont know how to do it.Could you please guide me in it. Thanks in advance Nampreet _________________________________________________________________ All that you wanted to know about Ms Beautiful Lips http://server1.msn.co.in/Profile/katrina.asp