Re: [R] FFT Normalization Documentation

2015-02-05 Thread Martin Maechler
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

[R] Parallel Computing in R

2015-02-05 Thread Zhenchuan Wang
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

Re: [R] Parallel Computing in R

2015-02-05 Thread Pascal Oettli
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 zwan...@mtu.edu wrote: I usually need to compute P-values as following: 1. generate one sample (usually

[R] symm parameter of heatmap

2015-02-05 Thread Aziz, Muhammad Fayez
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

Re: [R] save structure to be accesible later

2015-02-05 Thread Olivier Crouzet
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 -

Re: [R] save structure to be accesible later

2015-02-05 Thread Ivan Calandra
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: [R] FFT Normalization Documentation

2015-02-05 Thread Hiroyuki Kawakatsu
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:

[R] save structure to be accesible later

2015-02-05 Thread Alaios via R-help
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-es] Interpretación de coeficientes en un cox proportional hazards con variable strata

2015-02-05 Thread José Luis Cañadas Reche
Gracias por los comentarios Fran. El tema era porque si no estratificaba por país sino que consideraba la interacción con movilidad no se me cumplían las hipótesis de hazard proportionals al estratificar por Q6 pero no por país. Me sale que hay diferentes hazards basales por país y por Q6,

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

2015-02-05 Thread David Winsemius
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 code

Re: [R] R2BayesX

2015-02-05 Thread Shanley Chong
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] Fixed Width EBCDIC Files in R

2015-02-05 Thread Brian Trautman
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

[R] Rotate array by 90°

2015-02-05 Thread Karim Mezhoud
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 draw a legend outside of the plot

2015-02-05 Thread Jim Lemon
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 kathryn.lord2...@gmail.com wrote: Dear R users, I have three plots, so I tried, for exmple, par(mfrow=c(2,2)) y1 -

Re: [R] R2BayesX

2015-02-05 Thread Achim Zeileis
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 download and unzip data in a loop

2015-02-05 Thread Alexandra Catena
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

Re: [R] Rotate array by 90°

2015-02-05 Thread Jeff Newmiller
?aperm aperm(A, c(2,1,3) )[,3:1,] --- Jeff NewmillerThe . . Go Live... DCN:jdnew...@dcn.davis.ca.usBasics: ##.#. ##.#. Live Go...

Re: [R] Rotate array by 90°

2015-02-05 Thread Rui Barradas
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] Rotate array by 90°

2015-02-05 Thread Karim Mezhoud
Thanks Jeff and Rui. On Thu, Feb 5, 2015 at 11:12 PM, Rui Barradas ruipbarra...@sapo.pt 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)))

Re: [R] Fixed Width EBCDIC Files in R

2015-02-05 Thread Brian Trautman
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

Re: [R] Fixed Width EBCDIC Files in R

2015-02-05 Thread John McKown
On Thu, Feb 5, 2015 at 2:08 PM, Brian Trautman btrautma...@gmail.com 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

Re: [R] Rotate array by 90°

2015-02-05 Thread JS Huang
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] Scraping HTML using R

2015-02-05 Thread Steve Lianoglou
You want to take a look at rvest: https://github.com/hadley/rvest On Thu, Feb 5, 2015 at 2:36 PM, Madhuri Maddipatla madhuri@gmail.com 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

[R] Scraping HTML using R

2015-02-05 Thread Madhuri Maddipatla
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

[R] Creating two autocorrelated and correlated time series

2015-02-05 Thread Daniel Marmander
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

Re: [R] one sample test

2015-02-05 Thread meng
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 pda...@gmail.com wrote: On 05 Feb 2015, at 09:20 , meng laomen...@163.com wrote: Hi all: If I want to test whether the mean

Re: [R] FFT Normalization Documentation

2015-02-05 Thread Martin Maechler
Hiroyuki Kawakatsu hkawa...@gmail.com 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,

[R] Needing help about recosystem package

2015-02-05 Thread PO SU
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

Re: [R] Rotate array by 90°

2015-02-05 Thread Chel Hee Lee
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 03:17 PM, Karim Mezhoud wrote: Dear

[R] one sample test

2015-02-05 Thread meng
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

[R] Rolling window VAR (vector autoregressive) estimation

2015-02-05 Thread Wesley Burnett
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

Re: [R] one sample test

2015-02-05 Thread peter dalgaard
On 05 Feb 2015, at 09:20 , meng laomen...@163.com 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

Re: [R] Rhive installation problem in Windows RStudio

2015-02-05 Thread Prof Brian Ripley
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,

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

2015-02-05 Thread Munjal, Aarti
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

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

2015-02-05 Thread 暮如雪
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

Re: [R] Rhive installation problem in Windows RStudio

2015-02-05 Thread Pascal Oettli
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 praka...@hcl.com wrote: Hi, I am facing problem in RHive installation in Rstudio Version 0.98.1091 for R

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

2015-02-05 Thread Jon Skoien
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

Re: [R-es] Interpretación de coeficientes en un cox proportional hazards con variable strata

2015-02-05 Thread Francisco Viciana
Me corrijo a mi mismo: Si puede tener sentido sentido estratificar por la variable País, si no tenemos interés en estimar el efecto de especifico asociado al país, si no en su iteración con otra variable de exposición. Para aclararme las ideas he construido unos datos ficticios sobre los que poder

[R] Rhive installation problem in Windows RStudio

2015-02-05 Thread deepu
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 --

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

2015-02-05 Thread jim holtman
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 amc5...@gmail.com wrote: Hi All, I need to loop through and download

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

2015-02-05 Thread Camilo Mora
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),

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

2015-02-05 Thread Steve Lianoglou
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

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

2015-02-05 Thread Jeff Newmiller
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