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 >

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 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 that contents w

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] Gender balance in R

2014-11-25 Thread Scott Kostyshak
On Tue, Nov 25, 2014 at 1:15 PM, Martin Morgan wrote: > On 11/25/2014 04:11 AM, Scott Kostyshak wrote: >> >> On Mon, Nov 24, 2014 at 12:34 PM, Sarah Goslee >> wrote: >>> >>> I took a look at apparent gender among list participants a few years ago: >>> https://stat.ethz.ch/pipermail/r-help/2011-Ju

Re: [R] Gender balance in R

2014-11-25 Thread Scott Kostyshak
On Tue, Nov 25, 2014 at 8:24 AM, Maarten Blaauw 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, black=unknown). After a c. 5% f

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

2014-11-25 Thread Michael Mason
Thanks! That worked From: William Dunlap mailto:wdun...@tibco.com>> Date: Tuesday, November 25, 2014 12:53 PM To: Rolf Turner mailto:r.tur...@auckland.ac.nz>> Cc: Michael Mason mailto:mma...@benaroyaresearch.org>>, R help mailto:R-help@r-project.org>> Subject: Re: [R] plot.hclust point to older

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 persis

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 th

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,

[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), Stre

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 peo

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 wrote: > On 25/11/2014 3:12 PM, Tom Wright wrote: > >> Hi, >> This problem has me stumped so I thought I'd a

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
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] 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 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 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 .Intern

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

[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. require(tables)

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 o

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

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

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 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 li

Re: [R] RStuido seg faults

2014-11-25 Thread Berend Hasselman
> On 25-11-2014, at 00:58, Stack Kororā 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 help me fix RS

[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 .RD

[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 b

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" wrote: > >On 24 Nov 2014, at 18:34 , Sarah Go

[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 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] 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
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] 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? https://support.rstudio.com/hc/communities/public/ques

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

Re: [R] Gender balance in R

2014-11-25 Thread peter dalgaard
On 24 Nov 2014, at 18:34 , Sarah Goslee 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 any sign that that

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

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 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", varying=list(3:6), >>

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 Ntirengany

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==6

[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] Gender balance in R

2014-11-25 Thread Scott Kostyshak
On Mon, Nov 24, 2014 at 12:34 PM, Sarah Goslee 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 any sign that

[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: water_blnce=function(data,capaci

Re: [R] Gender balance in R

2014-11-25 Thread Rainer M Krug
Sarah Goslee 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 changed in the last thr

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 thr