[R] symm parameter of heatmap

Hi, What exactly is the difference between using symm=TRUE and symm=FALSE in a heatmap. More elaborately, when symm=FALSE, what does the heatmap plots because it seems like it changes the scale and makes it asymmetric, even if x is symmetric and square. Additionally, is it possible to view

Re: [R] Parallel Computing in R

Hi, You already asked this here (http://stackoverflow.com/questions/28357210/parallel-computing-in-r) and got some comments. Regards, Pascal On Fri, Feb 6, 2015 at 11:58 AM, Zhenchuan Wang wrote: > I usually need to compute P-values as following: > 1. generate one sample (usually it is a matri

[R] Parallel Computing in R

I usually need to compute P-values as following: 1. generate one sample (usually it is a matrix) 2. apply several methods (I already wrote a subfunction for each method, and they are independent) to the generated sample to get pvalues. 3. compare the pvalues. Since each method mentioned above ta

Re: [R] Rotate array by 90°

> lapply(1:2, function(x) t(A[rev(1:3),,x])) [[1]] [,1] [,2] [,3] [1,] "g" "d" "a" [2,] "h" "e" "b" [3,] "i" "f" "c" [[2]] [,1] [,2] [,3] [1,] "p" "m" "j" [2,] "q" "n" "k" [3,] "r" "o" "l" > Is this what you are looking for? I hope this helps. Chel Hee Lee On 02/05/2015

Re: [R] Scraping HTML using R

You want to take a look at rvest: https://github.com/hadley/rvest On Thu, Feb 5, 2015 at 2:36 PM, Madhuri Maddipatla wrote: > Dear R experts, > > My requirement for web scraping in R goes like this. > > *Step 1* - All the medical condition from from A-Z are listed in the link > below. > > http:/

[R] Scraping HTML using R

Dear R experts, My requirement for web scraping in R goes like this. *Step 1* - All the medical condition from from A-Z are listed in the link below. http://www.webmd.com/drugs/index-drugs.aspx?show=conditions Choose the first condition say Acid Reflux(GERD-...) *Step 2 *- It lands on the this

Re: [R] Rotate array by 90°

Hi, Here is an implementation. > rot90 function(a,times=1) { row <- dim(a)[1] col <- dim(a)[2] dep <- dim(a)[3] if (times %% 2 == 1){t <- row; row <- col; col <- t} tempA <- array(NA, c(row,col,dep)) for (i in 1:dep) { temp <- a[,,i] for (j in 1:times) { temp <-

Re: [R] Fixed Width EBCDIC Files in R

First off, thank you very much for taking a look at this. I didn't know "raw=TRUE" would be necessary here. Unfortunately, I'm stuck with the embedded nulls in the source data at this point. If worst comes to worst, does R have a way to do something like -- 1. Read the entire file in as raw bi

Re: [R] Rotate array by 90°

Thanks Jeff and Rui. On Thu, Feb 5, 2015 at 11:12 PM, Rui Barradas wrote: > Hello, > > Try the following. > > > rot90 <- function(x, n = 1){ > r90 <- function(x){ > y <- matrix(rep(NA, prod(dim(x))), nrow = nrow(x)) > for(i in seq_len(nrow(x))) y[, i] <

Re: [R] Rotate array by 90°

Hello, Try the following. rot90 <- function(x, n = 1){ r90 <- function(x){ y <- matrix(rep(NA, prod(dim(x))), nrow = nrow(x)) for(i in seq_len(nrow(x))) y[, i] <- rev(x[i, ]) y } for(i in seq_len(n)) x <- r90(x) x }

Re: [R] Fixed Width EBCDIC Files in R

On Thu, Feb 5, 2015 at 2:08 PM, Brian Trautman wrote: > I'm trying to read some mainframe data encoded as EBCDIC into R, and am at > a loss. I'd like to avoid using an external program to convert the files, > since I'm operating in a corporate environment. > > You can find the example files at at

Re: [R] Rotate array by 90°

?aperm aperm(A, c(2,1,3) )[,3:1,] --- Jeff NewmillerThe . . Go Live... DCN:Basics: ##.#. ##.#. Live Go... Live: OO#.. Dead: OO#

Re: [R] R2BayesX

Many thanks Z and Nikolaus! Much appreciated! Have a great weekend! Kind regards, Shanley -Original Message- From: Achim Zeileis [mailto:achim.zeil...@uibk.ac.at] Sent: Thursday, 5 February 2015 8:01 PM To: Shanley Chong Cc: r-help@r-project.org; nikolaus.uml...@uibk.ac.at Subject: Re: [R

[R] Rotate array by 90°

Dear aal, Is there a way to rotate array or a cube of matrices by Y axis? MatLab example: A = cat(3,{'a' 'b' 'c';'d' 'e' 'f';'g' 'h' 'i'},{'j' 'k' 'l';'m' 'n' 'o';'p' 'q' 'r'}) A(:,:,1) = 'a''b''c' 'd''e''f' 'g''h''i' A(:,:,2) = 'j''k''l'

Re: [R] How to download and unzip data in a loop

On Feb 5, 2015, at 10:03 AM, Alexandra Catena wrote: > Thank you guys for the response. > > I'm trying to download the last ten years of meteorology data from a > weather station in Livermore from the URL: > ftp://ftp.ncdc.noaa.gov/pub/data/noaa/2015/724927-23285-2015.gz > The Livermore station

[R] Fixed Width EBCDIC Files in R

I'm trying to read some mainframe data encoded as EBCDIC into R, and am at a loss. I'd like to avoid using an external program to convert the files, since I'm operating in a corporate environment. You can find the example files at at the link below, with both ASCII and EBCDIC versions. Note that t

Re: [R] How to download and unzip data in a loop

Dunno. Try posting your current code that fixes the previously mentioned problems, but this time use plain text so the HTML doesn't corrupt it. Usually you can solve this kind of issue by executing one line at a time and looking at each result to make sure it is what you think it is. You can als

Re: [R] Fastest way to calculate quantile in large data.table

Not sure if there is a question in here somewhere? But if I can point out an observation: if you are doing summary calculations across the rows like this, my guess is that using a data.table (data.frame) structure for that will really bite you, because this operation on a data.table/data.frame is

[R] Fastest way to calculate quantile in large data.table

In total I found 8 different way to calculate quantile in very a large data.table. I share below their performances for future reference. Tests 1, 7 and 8 were the fastest I found. Best, Camilo library(data.table) v <- data.table(x=runif(1),x2 = runif(1), x3=runif(1),x4=runif(100

Re: [R] How to download and unzip data in a loop

Thank you guys for the response. I'm trying to download the last ten years of meteorology data from a weather station in Livermore from the URL: ftp://ftp.ncdc.noaa.gov/pub/data/noaa/2015/724927-23285-2015.gz The Livermore station code is 724927-23285. If I wanted to download data from 2005, the

[R] Creating two autocorrelated and correlated time series

I can create both correlated series and autocorrelated time series with well defined correlation factors. I am however having trouble creating two series given fixed autocorrelation for some number of lags in the series (autocorrelation might differ between the two as well), that also has a specifi

Re: [R] one sample test

Very sorry that I ignore the question type.Yest,it's not an R question. Many thanks for your answer. -- QQ: 1733768559 At 2015-02-05 18:26:01,"peter dalgaard" wrote: > >On 05 Feb 2015, at 09:20 , meng wrote: > >> Hi all: >> If I want to test whether the mean of a set of normal distrib

[R] Needing help about recosystem package

Dear expeRt,   I am using recosystem package which is an wrapper of the LIBMF library,  it implements matrix factorize which needed in latent factor model, but the package only provides  the predict method, not include matrix P,Q which estimated. I just need the P and Q matrix, so ,if anyone h

Re: [R] FFT Normalization Documentation

> Hiroyuki Kawakatsu > on Thu, 5 Feb 2015 12:40:28 + writes: > On 2015-02-05 Martin Maechler wrote: > [...] >> Both a definition of the DFT (in LaTeX or directly as \deqn{}{}) which >> matches *our* fft(), and an extended commented example (ca 20 lines, >> rath

Re: [R] save structure to be accesible later

Hi, I generally prefer using the functions saveObject() and loadObject() from the R.utils package. I like that you load directly to an object in the R workspace. HTH, Ivan -- Ivan Calandra, ATER University of Reims Champagne-Ardenne GEGENAA - EA 3795 CREA - 2 esplanade Roland Garros 51100 Re

Re: [R] save structure to be accesible later

Hi, have a look at save() and load(). You can save part of your workspace to a file (save()) and load it from this file (load()). Of course the file format is R specific (but its specifications are free) You may easily share the file with other R users. Olivier. -- Olivier Crouzet LLING - Labo

[R] save structure to be accesible later

Dear all,I have a function that returns the following list. At the end I will call my function 1000 times and I want to keep for each of these 1000 "results" (the structure as given below)in an order to be accessible later (Load the 1000 results and access them within a for loop for example) How

Re: [R] FFT Normalization Documentation

On 2015-02-05 Martin Maechler wrote: [...] > Both a definition of the DFT (in LaTeX or directly as \deqn{}{}) which > matches *our* fft(), and an extended commented example (ca 20 lines, > rather than ca 100 lines!) would be welcome additions. [...] Here is a previous attempt: https://stat.ethz.

Re: [R] How to download and unzip data in a loop

In addition to following Jim's suggestion, you should probably also use full.names = TRUE, otherwise you will try to open a connection to files in your current directory, not in tmpdir. Another thing is that the unzipped files appear irregular with respect to columns, so read.table might not wor

Re: [R] How to download and unzip data in a loop

try taking the quotes off of 'files' Jim Holtman Data Munger Guru What is the problem that you are trying to solve? Tell me what you want to do, not how you want to do it. On Wed, Feb 4, 2015 at 5:24 PM, Alexandra Catena wrote: > Hi All, > > I need to loop through and download the past 10 yea

Re: [R] one sample test

On 05 Feb 2015, at 09:20 , meng wrote: > Hi all: > If I want to test whether the mean of a set of normal distributed data is > different from a value(e.g. 0), I can use one sample t test . But if the data > is not normal distributed,what kind of method should be used? > Not really an R quest

Re: [R] Rhive installation problem in Windows RStudio

On 05/02/2015 09:30, Pascal Oettli wrote: Hello, This page says that RHive is not available as Windows binaries: http://cran.r-project.org/web/packages/RHive/index.html and OS_type:unix explains why it cannot even be installed from the sources. Regards, Pascal On Thu, Feb 5, 2015

Re: [R] Rhive installation problem in Windows RStudio

Hello, This page says that RHive is not available as Windows binaries: http://cran.r-project.org/web/packages/RHive/index.html Regards, Pascal On Thu, Feb 5, 2015 at 3:35 PM, deepu wrote: > Hi, > > I am facing problem in RHive installation in Rstudio Version 0.98.1091 for R > version 3.1.2 on W

[R] How to unload R.dll successfully by FreeLibrary

hi, I've some C++ code which try to load R.dll and execute some R command, then unload R.dll, yet it seems fail to unload it, so second time to call the function hangs since it is trying to Initialize R again(this is not allowed as documented). Anyone can tell me why, thanks. //the following

[R] Rhive installation problem in Windows RStudio

Hi, I am facing problem in RHive installation in Rstudio Version 0.98.1091 for R version 3.1.2 on Windows machine. Warning message : > install.packages("RHive") Warning in install.packages : package ‘RHive’ is not available (for R version 3.1.2) Please help me its urgent Thanks Deepika --

[R] one sample test

Hi all: If I want to test whether the mean of a set of normal distributed data is different from a value(e.g. 0), I can use one sample t test . But if the data is not normal distributed,what kind of method should be used? Many thanks! My best. [[alternative HTML version delete

[R] Rolling window VAR (vector autoregressive) estimation

Hi, I am completely new to R. I have a dataset of several thousand, weekly observations with ten variables. I would like to create a program (algorithm) that does a rolling window VAR regression of window width "w" (for example, a window width of 200 weeks) for all of the variables within the m

[R] ASA John M. Chambers Statistical Software Award - 2015

John M. Chambers Statistical Software Award - 2015 Statistical Computing Section American Statistical Association The Statistical Computing Section of the American Statistical Association announces the competition for the John M. Chambers Statistical Software Award. In 1998 the Association fo

Re: [R] R2BayesX

Shanley, thanks for the example. Thank you for your help. Here is my R codes: data("MunichBnd") N <- length(MunichBnd); n <- N*5 dat <- data.frame(x1 = runif(n, -3, 3),id = as.factor(rep(names(MunichBnd), length.out = n))) dat$sp <- with(dat, sort(runif(N, -2, 2), decreasing = TRUE)[id])

Re: [R] how to draw a legend outside of the plot

Hi Kathryn, Try this at the end of your example: plot(y3,type="n",axes=FALSE,xlab="",ylab="") legend(35,1,c("y1","y2","y3"),pch=1:3) Jim On Thu, Feb 5, 2015 at 6:42 PM, Kathryn Lord wrote: > Dear R users, > > I have three plots, so I tried, for exmple, > > par(mfrow=c(2,2)) > > y1 <- rnorm(100)

Re: [R] FFT Normalization Documentation

Dear Jeff, Franklin, Eike and other interested parties, I agree that the help page for fft() deserves a bit more detail; its source is https://svn.r-project.org/R/trunk/src/library/stats/man/fft.Rd For instance, it has the comment %% %% Here, we should really have a nice \deqn{}{} giving t