Re: [R] More RGoogleDocs - R1C1 reference style

2014-11-25 Thread Jeff Newmiller
You appear to be using a package that are not even on CRAN, much less base R, so you need to tell us that in your (missing) reproducible example. It is self-described as primitive, so don't be surprised when it acts odd. Note that your creation of an exposed Google account in your previous

Re: [R] Gender balance in R

2014-11-25 Thread Rainer M Krug
Sarah Goslee sarah.gos...@gmail.com writes: I took a look at apparent gender among list participants a few years ago: https://stat.ethz.ch/pipermail/r-help/2011-June/280272.html Same general thing: very few regular participants on the list were women. I don't see any sign that that has

[R] Creating a new column inside a function.

2014-11-25 Thread Frederic Ntirenganya
Dear All, i need a help on how I can create a new column on my dataset and use it as argument inside the following function. The column i want to create and vary is Evapolation. It varies that'S why I need it as argument. When I make it like this is not working:

Re: [R] Gender balance in R

2014-11-25 Thread Scott Kostyshak
On Mon, Nov 24, 2014 at 12:34 PM, Sarah Goslee sarah.gos...@gmail.com wrote: I took a look at apparent gender among list participants a few years ago: https://stat.ethz.ch/pipermail/r-help/2011-June/280272.html Same general thing: very few regular participants on the list were women. I don't

[R] Error Missing values where true/false needed

2014-11-25 Thread Frederic Ntirenganya
Dear All, I am getting this error and don't know why it comes. can you please help ? Error in if (data$Rain[i_NA] == 60) { : missing value where TRUE/FALSE needed The loop is : indicNAs - which(data$Rain %in% NA) ind_nonleap = c() # NAs due to non leap years ind_nonrecord = c() # NAs due

Re: [R] Error Missing values where true/false needed

2014-11-25 Thread PIKAL Petr
Hi Error message seems to be clear Error in if (data$Rain[i_NA] == 60) { : missing value where TRUE/FALSE needed data$Rain[i_NA] produces probably NA x-NA if(x==60) print(1+1) else print(Errrorrr) Error in if (x == 60) 1 + 1 : missing value where TRUE/FALSE needed x-10 if(x==60)

Re: [R] Error Missing values where true/false needed

2014-11-25 Thread Michael Dewey
You do not tell us what you are trying to do but I think there is something wrong in the logic of your thinking as on the one hand you are selecting just precisely those elements of data$Rain which are NA and then testing whether any of them equals 60. On 25/11/2014 12:19, Frederic

Re: [R] More elegant way of stacking the data

2014-11-25 Thread Lee, Chel Hee
If you do not want to use the loop, a function called 'reshape' may be useful: df - data.frame(a=1:5,b=letters[1:5],c1=1:5,c2=2:6,c3=3:7,c4=4:8) out2 - reshape(data=df, direction=long, varying=list(3:6), times=paste(c,1:4,sep=)) out2 a b time c1 id 1.c1 1 a c1 1 1 2.c1 2 b c1 2

Re: [R] More elegant way of stacking the data

2014-11-25 Thread Dimitri Liakhovitski
Thanks a lot, guys! On Tue, Nov 25, 2014 at 9:42 AM, Lee, Chel Hee chl...@mail.usask.ca wrote: If you do not want to use the loop, a function called 'reshape' may be useful: df - data.frame(a=1:5,b=letters[1:5],c1=1:5,c2=2:6,c3=3:7,c4=4:8) out2 - reshape(data=df, direction=long,

Re: [R] Gender balance in R

2014-11-25 Thread peter dalgaard
On 24 Nov 2014, at 18:34 , Sarah Goslee sarah.gos...@gmail.com wrote: I took a look at apparent gender among list participants a few years ago: https://stat.ethz.ch/pipermail/r-help/2011-June/280272.html Same general thing: very few regular participants on the list were women. I don't see

[R] question regarding rmvDAG in package pcalg

2014-11-25 Thread Jack Luo
Hi, I am trying to use rmvDAG in pcalg package to generate data from DAG structure. One thing I found is that when the number of variables gets large, there can be really large numbers in the data matrix. I played around with different parameters and it looks like the same case. library(pcalg)

[R] RStuido seg faults

2014-11-25 Thread Stack Kororā
Greetings, I am having a big issue with RStudio segfaulting recently. It is becoming a very big problem for me. I have attached most of the information to the support site but no one has responded there. Can someone please help me fix RStudio?

Re: [R] Gender balance in R

2014-11-25 Thread Maarten Blaauw
Thanks for the responses so far. The gender ratio in R should reflect the gender ratio of the potential users, as this is the pool the R users / developers are coming from. I agree with this, but then again I don't think R really has 0% female users/developers as the R member list suggests.

[R] Converting list to character

2014-11-25 Thread Massimiliano Tripoli
Dear all, I can't convert the result of aggregate function in a dataframe. My data looks like: mydata - structure(list(ID = c(11, 11, 460, 460, 986, 986, 986, 986, 1251, 1251, 1251, 1251, 1251, 1251, 1251, 1251, 1801, 1801, 1801, 1801 ), YEAR = c(2009, 2010, 2010, 2011, 2008, 2009, 2010, 2011,

[R] Packages for Handling Large Data Set

2014-11-25 Thread jeeth ghambole
Hello All, I am working on BackTesting Strategies on stocks using daily prices. Initially the size of data was very limited and can be easily handled using R and SQL, but now my analysis has been extending on large set of data. Can anyone suggest me the best packages available for handling large

Re: [R] Gender balance in R

2014-11-25 Thread Maarten Blaauw
Nice graph, Scott, thanks! Based on your code I plotted not the absolute numbers but the ratios, which show slowly increasing relative participation of female Rhelpers over time (red = women, blue=men, black=unknown). After a c. 5% female contribution in 1998, this has grown to about 15% now.

[R] Checking the proportional odds assumption holds in an ordinal logistic regression using polr function

2014-11-25 Thread Charlotte Whitham
Dear list, I have used the ‘polr’ function in the MASS package to run an ordinal logistic regression for an ordinal categorical response variable with 15 continuous explanatory variables. I have used the code (shown below) to check that my model meets the proportional odds assumption following

Re: [R] Gender balance in R

2014-11-25 Thread Hatcher, Catherine
I just saw this comment and I agree with Peter. I have occasion to ask questions and get help on the R forum but I am not a programmer and use programs as I need them and I suppose I must comment more often. :) On 11/25/14, 11:28 AM, peter dalgaard pda...@gmail.com wrote: On 24 Nov 2014, at

[R] plot.hclust point to older version

2014-11-25 Thread Michael Mason
Hello fellow R users, I have recently updated to R 3.1.2. When trying to plot an hclust object to generate the dendrogram I get the following error: Error in .Internal(dend.window(n, merge, height2, hang, labels, ...)) : there is no .Internal function 'dend.window' I am indeed using R3.1.2

[R] porting an access database to sqlite

2014-11-25 Thread ravi
Hi,All my data is presently locked in a Microsoft access database. This has huge data in a number of large tables. Using RODBC and connecting to it takes too long a time, sometimes making the system to hang up. To make things more manageable, I have tried to transfer the data to manageable

Re: [R] RStuido seg faults

2014-11-25 Thread Berend Hasselman
On 25-11-2014, at 00:58, Stack Kororā i.am.st...@gmail.com wrote: Greetings, I am having a big issue with RStudio segfaulting recently. It is becoming a very big problem for me. I have attached most of the information to the support site but no one has responded there. Can someone please

Re: [R] Gender balance in R

2014-11-25 Thread Martin Morgan
On 11/25/2014 04:11 AM, Scott Kostyshak wrote: On Mon, Nov 24, 2014 at 12:34 PM, Sarah Goslee sarah.gos...@gmail.com wrote: I took a look at apparent gender among list participants a few years ago: https://stat.ethz.ch/pipermail/r-help/2011-June/280272.html Same general thing: very few regular

Re: [R] Converting list to character

2014-11-25 Thread Lee, Chel Hee
do.call(rbind, TAB$x) [,1] [,2] 1 GR.3.8 GR.3.8 2 GR.3.1 GR.3.8 4 GR.3.8 GR.3.8 5 GR.3.7 GR.3.7 6 GR.3.8 GR.3.8 7 GR.3.1 GR.3.8 9 GR.3.8 GR.3.8 10 GR.3.7 GR.3.7 11 GR.3.1 GR.3.1 12 GR.3.8 GR.3.8 13 GR.3.1 GR.3.8 15 GR.3.8 GR.3.8 16 GR.3.1 GR.3.1 17 GR.3.8 GR.3.8 18 GR.3.1 GR.3.8 20

Re: [R] Converting list to character

2014-11-25 Thread David L Carlson
Or just modify your aggregate() command: TAB - aggregate(mydata$CODE, by=list(ID=mydata$ID, +YEAR=mydata$YEAR), FUN=paste0, collapse=, ) TAB ID YEAR x 1 986 2008 GR.3.8 2 1251 2008 GR.3.1, GR.3.8 3 1801 2008 GR.3.8 411 2009 GR.3.7 5

Re: [R] Packages for Handling Large Data Set

2014-11-25 Thread Greg Snow
Look at the task view for High Performance Computing ( http://cran.r-project.org/web/views/HighPerformanceComputing.html) there is a section on packages for large memory and out-of-memory analyses. There are also sections on parallel computing which is one way to deal with large data if you have

Re: [R] porting an access database to sqlite

2014-11-25 Thread Jeff Newmiller
Counting rows is not something RODBC is supposed to do. That is a very basic SQL query that you can use RODBC to execute: SELECT COUNT(*) FROM yourtablename --- Jeff NewmillerThe .

Re: [R] plot.hclust point to older version

2014-11-25 Thread Rolf Turner
Reproducible example??? (I know from noddink about hclust, but I tried the example from the help page and it plotted without any problem.) cheers, Rolf Turner On 26/11/14 06:13, Michael Mason wrote: Hello fellow R users, I have recently updated to R 3.1.2. When trying to plot an hclust

[R] Presentation tables in R (knitr)

2014-11-25 Thread Tom Wright
Hi, This problem has me stumped so I thought I'd ask the experts. I'm trying to create a pretty summary table of some data (which patients have had what tests at what times). Ideally I'd like to knitr this into a pretty PDF for presentation. If anyone has pointers I'll be grateful.

Re: [R] RStuido seg faults

2014-11-25 Thread ~Stack~
On 11/25/2014 12:09 PM, Berend Hasselman wrote: If SL is Snow Leopard (no system mentioned by you) this belongs on R-SIG-Mac. Actually, I did mention it: My OS: Scientific Linux 6.6 But this is really for RStudio support. If RStudio needs fixing the RStudio people will have to do that. I

Re: [R] RStuido seg faults

2014-11-25 Thread ~Stack~
On 11/25/2014 12:01 PM, Mark Sharp wrote: Have you tried the current version of R, 3.1.2? I have not. I haven't had many issues in the past using what was in the EPEL repos. Let me take one of my dev boxes and give it a try. I will post back what I find. Thanks! signature.asc Description:

Re: [R] plot.hclust point to older version

2014-11-25 Thread Rolf Turner
On 26/11/14 08:53, Michael Mason wrote: Here you are. I expect most folks won't get the error. N = 100; M = 1000 mat = matrix(1:(N*M) + rnorm(N*M,0,.5),N,M) h = hclust(as.dist(1-cor(mat))) plot(h) Error in .Internal(dend.window(n, merge, height2, hang, labels, ...)) : there is no

Re: [R] Presentation tables in R (knitr)

2014-11-25 Thread Mark Sharp
Tom, If you are wanting PDF as your output, are you wanting to use LaTeX or Markdown with knitr. LaTeX will give you more options. You have not shown an attempt to use either for your table construction. Can you define what you mean by pretty? Is it the underscores in the column names that are

Re: [R] plot.hclust point to older version

2014-11-25 Thread William Dunlap
You probably have a local copy of an old version of plot.hclust or plot.dendrogram in your global environmenet or another package that masks the one in package:stats. E.g., I fired up R-2.14.2 and copied those 2 plot methods to .GlobalEnv and then saved by workspace when quitting R. I then fired

Re: [R] Presentation tables in R (knitr)

2014-11-25 Thread Tom Wright
Hi Mark, It is the underscores that are my issue, I'd prefer multiple level row titles: |ID1 |ID2 |Time1 |Time2 |Time1 |Time2 |OD |OS |OD |OS|OD |OS |OD |OS Height | 1| 1 | 1 | 1 | 1 | 1 | 1 | 1 Weight| 1| 1 | 1 | 1 | 1

Re: [R] Presentation tables in R (knitr)

2014-11-25 Thread Duncan Murdoch
On 25/11/2014 3:12 PM, Tom Wright wrote: Hi, This problem has me stumped so I thought I'd ask the experts. I'm trying to create a pretty summary table of some data (which patients have had what tests at what times). Ideally I'd like to knitr this into a pretty PDF for presentation. If anyone has

Re: [R] Presentation tables in R (knitr)

2014-11-25 Thread Tom Wright
Thanks Duncan, Dropping the extra columns might be the way forward. I'm sure I can work out how to embed latex into a markdown document ;-) On Tue, Nov 25, 2014 at 4:37 PM, Duncan Murdoch murdoch.dun...@gmail.com wrote: On 25/11/2014 3:12 PM, Tom Wright wrote: Hi, This problem has me stumped

Re: [R] Presentation tables in R (knitr)

2014-11-25 Thread Jeff Newmiller
Pebbles in sand are also an alternative to Excel, but that doesn't mean people want to switch to abaci. Horse, meet water. If you use LaTeX (or Rmarkdown-to-pdf in RStudio), then you can obtain much better looking tables using latex.tabular(). Unfortunately, LaTeX is just too scary for some

[R] covariate or predictor

2014-11-25 Thread Kristi Glover
Hi, I am wondering how I can separate whether it is covariate or predictor in the ANOVA analysis. For example A-structure(list(Machine = c(1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L), Diameter = c(20L, 25L, 24L, 25L, 32L, 22L, 28L, 22L, 30L, 28L, 21L, 23L, 26L, 21L, 15L),

Re: [R] covariate or predictor

2014-11-25 Thread Rolf Turner
On 26/11/14 13:57, Kristi Glover wrote: Hi, I am wondering how I can separate whether it is covariate or predictor in the ANOVA analysis. For example A-structure(list(Machine = c(1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L), Diameter = c(20L, 25L, 24L, 25L, 32L, 22L, 28L, 22L,

Re: [R] covariate or predictor

2014-11-25 Thread Bert Gunter
Yes, Rolf -- she seems to think that covariates must be categorical and predictors categorical -- or maybe it's vice-versa. Anyway, she apparently has not done any homework (e.g. by reading an Intro to R) and so doesn't understand the use of modeling formulas in lm() and thus does not understand

Re: [R] covariate or predictor

2014-11-25 Thread Rolf Turner
On 26/11/14 15:49, Bert Gunter wrote: Yes, Rolf -- she seems to think that covariates must be categorical and predictors categorical -- or maybe it's vice-versa. You of course meant ... covariates must be categorical and predictors numerical -- or maybe it's vice-versa. (I can't help

Re: [R] plot.hclust point to older version

2014-11-25 Thread Michael Mason
Thanks! That worked From: William Dunlap wdun...@tibco.commailto:wdun...@tibco.com Date: Tuesday, November 25, 2014 12:53 PM To: Rolf Turner r.tur...@auckland.ac.nzmailto:r.tur...@auckland.ac.nz Cc: Michael Mason mma...@benaroyaresearch.orgmailto:mma...@benaroyaresearch.org, R help

Re: [R] Gender balance in R

2014-11-25 Thread Scott Kostyshak
On Tue, Nov 25, 2014 at 8:24 AM, Maarten Blaauw maarten.bla...@qub.ac.uk wrote: Nice graph, Scott, thanks! Based on your code I plotted not the absolute numbers but the ratios, which show slowly increasing relative participation of female Rhelpers over time (red = women, blue=men,

Re: [R] Gender balance in R

2014-11-25 Thread Scott Kostyshak
On Tue, Nov 25, 2014 at 1:15 PM, Martin Morgan mtmor...@fredhutch.org wrote: On 11/25/2014 04:11 AM, Scott Kostyshak wrote: On Mon, Nov 24, 2014 at 12:34 PM, Sarah Goslee sarah.gos...@gmail.com wrote: I took a look at apparent gender among list participants a few years ago:

Re: [R] Error Missing values where true/false needed

2014-11-25 Thread Frederic Ntirenganya
Hi PIKAL, The error seems to be starnge to me because i access the indices of NAs. Indices can't be non-applicable. This is the output of indecs having the NA in my dataset. my dataset is very big that's why I did not provide it. indicNAs - which(data$Rain %in% NA) indicNAs [1] 426 792

Re: [R] Error Missing values where true/false needed

2014-11-25 Thread Boris Steipe
On Nov 26, 2014, at 1:27 AM, Frederic Ntirenganya ntfr...@gmail.com wrote: The error seems to be starnge to me because i access the indices of NAs. No you don't. You access the contents of the cell via an index for which you have previously determined that the contents is NA. Then you compare

Re: [R] Error Missing values where true/false needed

2014-11-25 Thread PIKAL Petr
Huh and what is the result of e.g. data$Rain[indicNAs[1]] I bet that you will see [1] NA and your if question askes if (NA eqals 60) which results in NA and if is telling you that it expects TRUE or FALSE but not NA. I do not see how much clearer the error message shall be. Cheers Petr

[R-es] help: Shiny y denisty plots en R maps

2014-11-25 Thread Javier Villacampa González
Hola buenas, estoy en una época convulsa de de dudas, ya que ando pegandome con shiny y los mapas. Quería pegar un heatmap dinámico (mapa de paises) en una aplicación shiny. El problema es que no termino de hayar la manera. Os paso las apps. A ver si alguien sabe algo. La idea es poder mostrar el

[R-es] Duda sobre cómo analizar un experimento factorial con algoritmos de extracción de características, clustering y clasificación como factores

2014-11-25 Thread Daniel Carrillo Zapata
Hola compañeros Soy Daniel Carrillo, y os escribo porque me ha surgido una duda sobre si puedo tratar algoritmos de clustering como un factor en un experimento. Concretamente, tengo un conjunto de datos sin etiquetar, y quiero probar los siguientes algoritmos sobre él: 1) Extracción de