similar to: reading tables into R

Hi, The file I am reading is a text file, whose contents are a matrix that has 15 rows and 58 columns. The first row has column names, and the first column has row names, so the format is correct as far as using read.table is concerned. The other values in the table are all float values (numeric). So when I read in the file using data1 <- read.table("HAL001_HAL0015_Signals.txt"), it

Hi, I have been using the prcomp function to perform PCA on my example microarray data, (stored in metric text files) which looks like this: 1a 1b 1c 1d 1e 1f ...................................................4r 4s 4t g1 1.2705 1.2766 ...........................................................2.0298 g2 0.1631

Hi, I've just started using R and RSPerl. I have some code as follows: &R::initR("--no-save"); &R::call("t.test", (\@array1, \@array2)); where @array1 and @array2 are both 1-dimensional arrays in Perl having 54675 elements each. On execution the output is as follows: Calling R function name `t.test', # arguments: 3 1) Arg type 3 Got a reference to a

I have passwords turned off, and require keys to match. The zombie armies swarming outside are trying brute force attacks that in part involve guessing login NAMES. If they guess the wrong NAME, this is logged in syslog. If they guess a working user name, then the attack has PARTIALLY SUCCEEDED, but this information is IGNORED. That is, it is not logged. If the zombie army has tell when it

What is the status of i7 architecture support for CentOS-5? Do the latest updates support it? -- This e-mail may contain technical information which is controlled by the United States Government, Department of State, International Traffic & Arms Regulations (ITAR) (22 CFR 120-130) which requires an export license prior to sharing with foreign persons. Lacking such license, ITAR technical

I've got a Solaris 2.6 SPARC system running Samba 2.0.6 and an NT-4.0sp5 PDC system. If Samba is running on the Sun system when the NT PDC is rebooted it (the NT system) complains that there is already another PDC in the domain. This should not be as far as I can see as I've tried specifically to configure the Samba systems not to do this and in fact prior to upgrading to Samba version

I am running Red Hat 7.2 and using smbmount to mount a PC share from a NT4 SP5 Workstation. I have a folder, in this share, of around 300 files and when I do ls -l I get inconsistencies. Sometimes files are missing, and if I keep repeating the command the files are then appear again. Other times the size information is wrong. No one else is using, or changing, any of the files on the source

Does anyone know of an ACL/extended attributes enabled version of the Red Hat 9 kernel exist in RPM form? -- Michael St. Laurent Hartwell Corporation

Hi all, How can I overcome the error "object of type 'closure' not subsettable" I ran the following script seq <- paste(seq(1914, 1916, by=1), "*.y", sep=".") # make sequence c <- 3 # total number of files d2 <- file # creates dummy file # Input sequence in loop for (i in 1:3){ list <- list.files("~/ukcp09/txt/x.djf", seq[[i]]) file

I came across this in connection with an unrelated issue > beta[2] Error in beta[2] : object is not subsettable > beta[2] <- 5 Error in "[<-"(`*tmp*`, 2, value = 5) : object is not subsetable One of the messages must be wrong, but I need a native English speaker to tell me which one. -- O__ ---- Peter Dalgaard Blegdamsvej 3 c/ /'_ --- Dept. of

Dear all, Is there a less cumbersome way to rename a column by name (as opposed to index) than -- names( X)[ names[ X] == "bob"]<-"sue" ? A semi-related question: how does one get the index of a column by name, something along the lines of col.index( X, "sue") ? Chip Barnaby --------------------------------------------------------- Chip Barnaby

I'm seeing random disk full errors when saving large Excel files to a Samba share. The users will retry and the second or third time (usually the second time) it will save correctly. This is with version 3.0.21a on i386 Red Hat 9, kernel 2.4.20-31 patched with ACL. -- Michael St. Laurent Hartwell Corporation "This is the captain speaking. We may experience some turbulence... and

Dear R-users I need to calibrate kappa, rho, eta, theta, v0 in the following code, see below. However when I run it, I get: y <- function(kappahat, rhohat, etahat, thetahat, v0hat) {sum(difference(k, t, S0, X, r, implvol, q, kappahat, rhohat, etahat, thetahat, v0hat)^2)} > nlminb(start=list(kappa, rho, eta, theta, v0), objective = y, lower =lb, > upper =ub) Error in dots[[1L]][[1L]] :

Hello, I would like to create a subscriptable collection (presumably a list) of lm() models. I have a data frame DX containing 6 groups of data. The general idea is (NOT RUN) ... for (i in 1:6) { DXS = subset( DX, <whatever>); LMX[ i] = lm( <formula>, data = DXS); } Now access model results by subscript ... e.g. coefficients( LMX[ 2]). Or would it be [[ 2]]? I

Hi everyone, Would somebody please explain (or point me to a reference that explains) the following error: "Error: object of type 'closure' is not subsettable" I was trying to use rep() to replicate a function: > example_function <- function() { return(TRUE) } > rep(example_function, 3) Error: object of type 'closure' is not subsettable But I just cannot

Hola, Alguien podría ayudarme en la siguiente cuestión: Trato de abrir múltiples archivos ncDF para dejarlos disponibles para trabajar luego con ellos. Intento hacerlo con un blucle pero no consigo lo que quiero. El código es el siguiente: setwd("/Users/jaime/Desktop/2008") a<-list.files() a library(ncdf) for (i in 1:length(a)) {

Hello R list, I have a little problem with splom. I'd like to wrap it in a function, for example: multi.scatterplot <- function(data,groups,cols,colors) { splom(~data[,cols], groups = as.symbol(groups), data = data, panel = panel.superpose, col=colors) } and then call it like in multi.scatterplot(iris,"Species",1:4,c("green","blue","red"))

Hi all, I want to write a wrapper for an analysis of variance and I face a curious problem. Here are two different wrappers: fun.1 <- function(formula) { summary(aov(formula)) } fun.2 <- function(formula) { oneway.test(formula) } values <- c(15, 8, 17, 7, 26, 12, 8, 11, 16, 9, 16, 24, 20, 19, 9, 17, 11, 8, 15, 6, 14) group <- rep(1:3, each=7) # While the first

Hi All, I have a data frame where nearly 10K columns of data, where most of them have standard deviation( of all rows) as zero. I want to exclude all the columns from the data frame and proceed to further processing. I tried like blow. *data <- read.csv("data.CSV", header=T) for(i in 2:ncol(data)) if(sd(data[,i])==0){ df[,i] <-NULL } * where I have the data columns from

Hi, I am attempting to write code which will read in my data which is of this form: X1 Y1 X2 Y2 .... Xn Yn 0 0 0 0 0 0 1 0 1 255 1 0 2 255 2 0 2 255 3 0 3 0