similar to: reading tables into R. .

I have been using R and Perl. When I read in a text file using the read.table option, and I try to mathematically manipulate the individual elements in the table, I keep getting an "object is not subsettable" error. If I try to use a different method, it works, but takes too much time(basically, I then need to read in values individually into R instead of as a 2D array, so the number of

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

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

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

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

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]] :

Recently, I ran a series of Kruskal-Wallace tests [kruskal.test()] using by() to group by site Output is a list: >Herb.KW Herb.df$ID: 10-1 Kruskal-Wallis rank sum test data: Indicator_Rating by Year Kruskal-Wallis chi-squared = 15.24, df = 7, p-value = 0.03302 ----------------------------------------------------------------------------------------------------- Herb.df$ID: 18-1

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

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

hello everybody, I have a directory with over 3000 files with different names. I would like to make some vectors with the file names which are belong together. I'm trying to do it with a for loop in R: SF <- c("ad", "cd", "cer", "stress", "salty", "PC", "high", "transfer", "cold",

Hi all, when i run this script: >read.table("Angelika.txt",header=T,sep="\t") >mytable=read.table("Angelika.txt",header=T,sep="\t") >for ( dye in c("A","B","C","F","G","K","L","M")) + { + for (cond in 1:8) + { + measurement =

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"))

Dear Sir/Madam, I have just installed R for Windows. I am trying to open a stata file and I have a problem. I have used the following commands: > install.packages("foreign") > library(foreign) > data.name <- read.dta(file.choose()) # then I choose a kelleya.dta file and click on open and then I use the following command: > summary(kelleya) Error in object[[i]]

Hello, I am a newbie. I can run the following code stored in "test.txt" without error using my XP machine: x <- scan("C:\\Rwork\\A.txt") x10 = filter(x, rep(1/10,10), sides=1) x x10 for(i in 10:length(x)){ if (x[i] > x10[i]) diff[i]="b" else diff[i]="s" } However, if I run it in another PC that uses Windows 7, the following error occurs:

I'm attempting to write a general function to implement Faraway's bootstrapping algorithm for mixed models with lmer, but have run into a curious problem. I'm comparing two models model.1<-lmer(Response ~ Treatment + (1|Trial), data=exp.data, method="ML") model.2<-lmer(Response ~ 1 + (1|Trial), data=exp.data, method="ML") When I attempt to update

Hello folks: I have a data frame which includes two numeric vectors: "Pa1" and "Ingr" (which represents income).  I run the following code on the R console: ingreso <- factor(ifelse(Ingr<=3000,"bajo","alto")) tapply(Pa1,ingreso,summary) boxplot(Pa1~ingreso) and it works fine.  However, when I "Sweave" through an rnw file containing that

hey there! The object 'cit' contains: > cit ##################################################### # Johansen-Procedure Unit Root / Cointegration Test # ##################################################### The value of the test statistic is: 5.3484 9.0681 10.6433 --------------- I want R to save the value 5.3484 in a new object. I am used to use the command x=cit[a] where a

Hello, I am trying to use the p.adjust function for multiple testing. here is what i have 9997 201674_s_at 0.327547396 9998 221013_s_at 0.834211067 9999 221685_s_at 0.185099475 I import them from excel have have the gene symbol as well as the pvalue here is the issue > pa<-p.adjust(pt,method="BH") Error in p[nna] : object is not