Re: [R] C stack error in as.vector() starting in R 3.3.0

Luke, Thanks for this point out. I started removing packages one-by-one in R v.3.3.1. Oddly enough, it went away after I removed the Matrix package. When I restarted R, the same version of Matrix (1.2-6) had been "reinstalled" (?), but now no C stack error. Lets see if the others that were having

Re: [R] GAMS, std errors and confidence intervals

> On Jul 4, 2016, at 11:44 AM, Dan Jaffe wrote: > > Can anyone help me calculating CIs from a GAM analysis? > > I have calculated a GAM fit (m3) and the associated std errors using > predict.gam > I assume that the 95% CI around each fit value would be 1.96 > times the se.But when I do this

[R] Checking for modality (R - circular)

R Users Is there a way to check for modality using the "circular" package in R or any other package ? Circular time data. Jeff [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see

Re: [R] C stack error in as.vector() starting in R 3.3.0

I cannot reproduce this. My best guess is that there is a problem, maybe a version incompatibility, with one of the packages loaded when load("problem.method.rdata"). This includes adegenet apex copula Rcpp DBI sp colorspace and any dependencies these are bringing in. Best, luke On Fri, 1 Ju

Re: [R] BCa Bootstrapped regression coefficients from lmrob function not working

Dear Duncan, Many thanks for your reply. What you propose is a very good idea yes. However, I am interested in the coefficients... Best, S De : Duncan Murdoch Cc : R-help Mailing List Envoyé le : Mercredi 6 juillet 2016 14h13 Objet : Re: [R] BCa Bootstr

[R] Performing Principal Cluster Analysis, k-means clustering etc. with PDB/DCD trajectory

Hello If someone knows how to convert an MD trajectory into a format acceptable to R or has done any of the above analyses, please let me know. Thanks and regards Sourav [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] dplyr : row total for all groups in dplyr summarise

Thanks muchly. I hate the smart quotes! -Roy > On Jul 6, 2016, at 11:42 AM, David Winsemius wrote: > > >> On Jul 6, 2016, at 9:45 AM, rmendelss gmail wrote: >> >> >>> On Jul 6, 2016, at 9:36 AM, David Winsemius wrote: >>> >>> n this case the text was cut from the R session console text

Re: [R] dplyr : row total for all groups in dplyr summarise

> On Jul 6, 2016, at 9:45 AM, rmendelss gmail wrote: > > >> On Jul 6, 2016, at 9:36 AM, David Winsemius wrote: >> >> n this case the text was cut from the R session console text and pasted >> without modification into Mail.app version 8.2. In replicating this action, >> I see now that hitti

Re: [R] dplyr : row total for all groups in dplyr summarise

My point is that this is highly software-dependent. Certain email programs and editors are worse than others in inclusion of configuration settings that allow you to avoid this problem. In general you need to look for "plain text" options, and some software has "Auto-Correct" options turned on b

Re: [R] dplyr : row total for all groups in dplyr summarise

> On Jul 6, 2016, at 9:36 AM, David Winsemius wrote: > > n this case the text was cut from the R session console text and pasted > without modification into Mail.app version 8.2. In replicating this action, I > see now that hitting "return" then unfortunately converts the final > double-quote

Re: [R] How to extract "specific"/"last" intercept value from segmented package.

Hello, I'm glad that it helped. Note, however, that you don't need inter.m <- as.matrix(...) because inter$x already is a matrix. You can simply do inter.row <- nrow(inter$x) answer <- inter$x[inter.row,1] Rui Barradas   Citando Narendra Modi : > Thanks! that worked. > > I also tested with t

Re: [R] How to extract "specific"/"last" intercept value from segmented package.

Thanks! that worked. I also tested with the below method although your solution is faster and done in fewer steps. inter <- intercept(segmented.mod) inter.m <- as.matrix(inter$x) inter.row <- nrow(inter.m) answer <- inter.m[inter.row,1] PD On Tue, Jul 5, 2016 at 2:28 PM, wrote: > Hell

Re: [R] dplyr : row total for all groups in dplyr summarise

> On Jul 6, 2016, at 8:24 AM, Jeff Newmiller wrote: > > Cut and paste is not to blame... it is the use of word processing software > rather than text editors for manipulating code that is the problem. In this case the text was cut from the R session console text and pasted without modificati

Re: [R] tclArray() --- "active"?!

Thanks, Peter. When it comes time to read the array back to R, I can explicitly delete the "active" cell first... On Wed, Jul 6, 2016 at 8:35 AM, peter dalgaard wrote: > I am not up to speed on Tcl/Tk these days, but I would suspect that the > issue lies in tkTable and not in tclArray. The tabl

Re: [R] tclArray() --- "active"?!

I am not up to speed on Tcl/Tk these days, but I would suspect that the issue lies in tkTable and not in tclArray. The table widget has the notion of an active cell, which can be used in indexing operations. According to the widget documentation -- Command-Line Name:-variable Database Name: var

Re: [R] Formatting ggplot2 graph

Hi Georg, Your problem with the geom_text was that you tried to count Var1 in some bins which cannot be defined on characters - the counts you really wanted to use are in Anzahl. This should do what you want: library(ggplot2) freq_ls <- structure(list(Var1 = c("zldkkd", "aakdkdk","aaakdkd", "aai

Re: [R] WG: Fw: Re: dplyr : row total for all groups in dplyr summarise

Cut and paste is not to blame... it is the use of word processing software rather than text editors for manipulating code that is the problem. Georg: note that plyr does not mix very well with dplyr... try to pick one and stick with it. -- Sent from my phone. Please excuse my brevity. On Jul

Re: [R] WG: Fw: Re: dplyr : row total for all groups in dplyr summarise

It's the cut and paste monster. Somewhere along the way, the final " got converted to ” which R does not see. > "100%” + > "100%" [1] "100%" - David L Carlson Department of Anthropology Texas A&M University College Station, TX 77840-4352 -Original Messa

[R] tclArray() --- "active"?!

Sometimes when working with tclArray's, an empty element called "active" is appended to the array. Does anyone know when (and why) this happens? And how to prevent it (or at least predict it so that it can be removed by hand)? E.g., library(tcltk); library(tcltk2) tt<-tktoplevel() dat<-tclArray()

[R] Formatting ggplot2 graph

Hi All, my current code looks lke this: freq_ls <- structure(list(Var1 = c("zldkkd", "aakdkdk", "aaakdkd", "aaieiwo", "vöalsl", "ssddkdk", "glowowp", "laoiw", "ruklow", "rolsl", "delk

Re: [R] BCa Bootstrapped regression coefficients from lmrob function not working

On 06/07/2016 7:56 AM, varin sacha via R-help wrote: Dear Professor Dalgaard, Okay, this is what I was afraid of. Many thanks for your response. I know that my next question is off-topic here (on this website) but is it nevertheless possible to calculate confidence intervals for MARS regressio

Re: [R] BCa Bootstrapped regression coefficients from lmrob function not working

Dear Professor Dalgaard, Okay, this is what I was afraid of. Many thanks for your response. I know that my next question is off-topic here (on this website) but is it nevertheless possible to calculate confidence intervals for MARS regression or CIs for MARS will just be impossible to calculate

Re: [R] BCa Bootstrapped regression coefficients from lmrob function not working

Offhand, I would suspect that the cause is that earth() does not return a coefficient vector of the same length in every bootstrap iteration (_adaptive_ regression splines). This makes the bootstrap rather tricky to even define. -pd > On 06 Jul 2016, at 10:15 , varin sacha wrote: > > Dear Ber

[R] WG: Fw: Re: dplyr : row total for all groups in dplyr summarise

Hi All, if I run the suggested code mtcars %>% group_by (am, gear) %>% summarise (n = n()) %>% mutate(rel.freq = paste0(round(100 * n / sum(n), 0), "%")) %>% ungroup() %>% plyr::rbind.fill(data.frame(n = nrow(mtcars), rel.freq = "100%”)) I

Re: [R] BCa Bootstrapped regression coefficients from lmrob function not working

Dear Bert, You are right. > results=boot(data=newdata, statistic=boot.MARS, R=1000,formula=PIBparHab ~ > QUALITESANSREDONDANCE + competitivite + innovation) Erreur dans boot(data = newdata, statistic = boot.MARS, R = 1000, formula = PIBparHab ~ : le nombre d'objets à remplacer n'est pas multi