similar to: Tab Separated File Reading Error

Hi All, I wonder if there should be one character for quote= in read.table, i.e., > args(read.table) function (file, header = FALSE, sep = "", quote = "\"'", dec = ".", ... I have a file containing the following lines, 08248-GOTERM 3'-phosphoadenosine 5'-phosphosulfate biosynthetic process 08279-GOTERM 3'-phosphoadenosine

Hi there, I wonder if there is a way of efficiently generating a correlation matrix of two expression matrices. I want to correlate miRNA and mRNA expression and used the following code: ##dat.mi miRNA expression matrix, dat.m mRNA expression matrix nc <- nrow(dat.mi) cor.mat <- data.frame(rep(NA,nrow(dat.m))) pval.mat <- data.frame(rep(NA,nrow(dat.m))) for(i in 1:nc) { cr <- vector()

Dear All, I have two matrices A (40 x 732) and B (40 x 1230) and would like to calculate correlation between them.  I can use: cor(A,B, method="pearson") to calculate correlation between all possible pairs. But the issue is that there is one-many specific mappings between A and B and I just need to calculate correlations for those pairs (not all). Some variables in A (proteins, say p1)

Hi, After manipulate my data I have ended up with 5 different data frames with different number of observations but the same number of variables (columns) An example, if I write str(object1), I see this, data.frame': 47 obs. of 3 variables: $ ORF : Factor w/ 245 levels "YAL038W","YAL054C",..: 10 19 38 39 44 45 50 51 59 60 ... $ mRNA : num 0.891 1.148 1.202

Full_Name: Ross Darnell Version: 2.4.1 OS: Linux Submission from: (NULL) (130.102.133.33) rdarnell at pc:~/Data$ ls -al .RData -rwxrwxrwx 2 root root 13551 2006-12-06 08:58 .RData rdarnell at pc:~/Data$ R R version 2.4.1 (2006-12-18) Copyright (C) 2006 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to

Hi all Observe: x <- c(1,2) y <- c(1,-1) phi <- 1 p <- 2 par( mfrow=c(1,2)) plot(x , y, main=bquote( paste( p==.(p)," and ",phi==.(phi)) ) ) plot(y ~ x, main=bquote( paste( p==.(p)," and ",phi==.(phi)) ) ) par( mfrow=c(1,2)) On my system (details below), the first plot is correct (in my understanding), and produces a title reading "p=2 and phi=1"

Buenas tardes a tod@s, Estoy incursionando en el analisis de tweets utilizando el paquete twitteR y siguiendo http://www.webmining.cl/2012/07/text-mining-de-twitter-usando-r/ Desafortunadamente cuando ejecuto # cargar librerias library(twitteR) library(tm) library(wordcloud) # recolecta tweets de @camila_vallejo tweets = userTimeline("camila_vallejo", 2000) obtengo Error in

I have an S4 object class defined in a Bioconductor package which contains multiple slots, some of which are tibbles, whilst others are vectors. If I call slot(object, name) where 'name' is an slot that contains a vector, everything works as expected. However, when I call slot(object, name) where 'name' is an slot that contains a tibble I get the following warning: Warning

Hola a todos, Me gustaria leer el texto que se encuentra en http://dl.dropboxusercontent.com/u/9601860/txt.txt He intentado txt <- 'http://dl.dropboxusercontent.com/u/9601860/txt.txt' r <- scan(txt) #Error in scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, : # invalid multibyte string at '<ff><fe>M' r <- read.table(txt, header = FALSE)

Bug summary: glm() causes a segfault if the argument 'data' is a data frame with more than 16384 rows. Bug demonstration: -------input --------------- N <- 16400 df <- data.frame(x=runif(N, min=1,max=2),y=rpois(N, 2)) glm(y ~ x, family=poisson, data=df) ------ output --------------- *** caught segfault *** address (nil),

Hello, I have a minimal example that shows a problem I'm having with parallel processing. library(parallel) mcmapply(function(x, y) { print("Running") Sys.sleep(10) }, as.list(1:10), as.list(10:1), mc.cores = 16, SIMPLIFY = FALSE) I see "Running" printed once every ten seconds. I read the documentation for mcmapply, but I don't understand why it wouldn't

Hi there, using cbind with a numeric and raw argument produces an incorrect result. I've posted some details below, kind regards, Mike. e.g. > cbind(0, as.raw(0)) [,1] [,2] [1,] 0 6.950136e-310 A longer example shows that the result is not a rounding error, is not consistent, and repeated applications get different results. > cbind(0, as.raw(1:10))

The documentation states : An alternative parametrization (often used in ecology) is by the mean ?mu?, and ?size?, the dispersion parameter. However, this fails : > rnbinom(10, mu = 100, size = 0) [1] NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN Warning message: In rnbinom(10, mu = 100, size = 0) : NAs produced For dispersion set to 0, it should work like drawing from a Poisson distribution.

Hi, I'm trying to use Frank Harrell's Design:::ols function to do regression of y (numeric) on the interaction of two factors (x1 and x2), but Lapack throws an error: > library(Design) ... > load(url("http://www.csse.unimelb.edu.au/~gabraham/x")) > ols(y ~ x1 * x2, data=x) Error in chol2inv(fit$qr$qr) : 'size' cannot exceed nrow(x) = 20 > traceback()

Hi all I have just downloaded the sma package from CRAN. On installing on my linux machine, I get the message > library(sma) > data(MouseArray) Warning message: file 'MouseArray.RData' has magic number 'RDX1' Use of save versions prior to 2 is deprecated Hereafter, MouseArray is not found: > MouseArray Error: object "MouseArray" not found We were

Dear list, I'm trying to suppress/redirect/squash the output from commands like install.packages, or download.file. The problem is that none of: sink(..., type="message"), sink(..., type="output"), capture.output, suppressMessages are quite doing the trick. Output gets written to the stderr stream, despite any combination of the above suppression commands. According to

Muchas gracias Marcelino. Desafortunadamente al ejecutar el ejemplo que aparece en ?registerTwitterOAuth obtengo twitCred$handshake() Error: Unauthorized registerTwitterOAuth(twitCred) Error in registerTwitterOAuth(twitCred) : oauth has not completed its handshake Alguna otra sugerencia? Muchas gracias! Jorge.- 2013/7/25 Marcelino de la Cruz <marcelino.delacruz@upm.es> > Mira

Buenos dias a todos, Me gustaria genear n muestras de una N(0,1) a determinados intervalos de tiempo. Asi por ejemplo, si el vector de tiempos es tiempos <- c(2, 4, 6) la primera columan debe generarse a los 2 segundos de ejecutado el script, la 2a columna a los 4 segundos y la ultima a los 6 segundos. Al final el tiempo total de ejecucion sera aproximadamente 12 segundos para completar

I'm observing the following on different platforms: > parse(text='"\\x7F"') expression("\177") > parse(text='"\\x80"') Error: invalid multibyte string ... > parse(text='"\\xFF"') Error: invalid multibyte string However, cat("\x7F\n\x80\n...\xFF\n") works. Using R --vanilla. SYSTEMS GIVING THE ERROR: >

Wanted to raise two questions: 1. Is bugs.r-project.org down? I haven't been able to reach it for two or three days: ``` ping bugs.r-project.org PING rbugs.research.att.com (207.140.168.137): 56 data bytes Request timeout for icmp_seq 0 Request timeout for icmp_seq 1 Request timeout for icmp_seq 2 Request timeout for icmp_seq 3 Request timeout for icmp_seq 4 Request timeout for icmp_seq 5