[R] reticulate via callr

2025-04-09 Thread Benjamin Tyner
Greetings, I would like to programmatically interact with multiple versions of a given Python module via the reticulate package, without having to restart the R process over and over. (My understanding is reticulate does not provide such a capability by itself, as once reticulate binds to a P

Re: [R] installation: while running make, unable to run pdflatex on 'NEWS.tex'

2024-03-02 Thread Benjamin Tyner
r to that date. On 3/2/24 02:42, Ivan Krylov wrote: В Fri, 1 Mar 2024 10:46:53 -0500 Benjamin Tyner пишет: my platform info: Platform: x86_64-pc-linux-gnu (64-bit) Running under: Ubuntu 22.04.4 LTS Quick things first: have you installed all the build-dependencies? apt build-d

[R] dput(..., file = stderr())

2024-03-01 Thread Benjamin Tyner
Curious to know if this warning is expected behavior, and if so, what is the recommended way instead: > dput(letters, file = stderr()) c("a", "b", "c", "d", "e", "f", "g", "h", "i", "j", "k", "l", "m", "n", "o", "p", "q", "r", "s", "t", "u", "v", "w", "x", "y", "z") Warning messa

[R] installation: while running make, unable to run pdflatex on 'NEWS.tex'

2024-03-01 Thread Benjamin Tyner
A kind member of R-core suggested this is due to a misconfiguration on my system, and to post it to the mailing list for troubleshooting. When trying to build R version 4.3.3, in at least two places during the process it gives LaTeX errors of the form: (example 1) you should 'make docs

[R] .BaseNamespaceEnv locking in non-interactive mode

2022-11-23 Thread Benjamin Tyner
Greetings, Noticed that starting with R version 4.1.0, the base environment is locked by default. So for example in interactive mode, we have > .BaseNamespaceEnv$foo <- "bar" Error in .BaseNamespaceEnv$foo <- "bar" :   cannot add binding of 'foo' to the base environment However it appears tha

Re: [R] chaining closure arguments on-the-fly

2020-06-21 Thread Benjamin Tyner
On 6/20/20 5:04 PM, Duncan Murdoch wrote: I think you effectively did that in your original post (all but encapsulating the expression in a function), so yes, it's possible. However, it's a really bad idea.  Why use non-standard evaluation when standard evaluation is fine?  Standard evaluation

Re: [R] chaining closure arguments on-the-fly

2020-06-20 Thread Benjamin Tyner
On 6/20/20 5:49 PM, Bert Gunter wrote: Gents: (with trepidation) f(x = 3, y = g(expr)) **already** evaluates g in the environment of f, **not** in the environment of the caller. (This does not contradict Duncan's example -- 3 is a constant, not a variable). e.g. > f <- function(x = 3, y =

Re: [R] chaining closure arguments on-the-fly

2020-06-20 Thread Benjamin Tyner
On 6/20/20 9:00 AM, Duncan Murdoch wrote: How about g <- function(x, y = x) {   f(x, y) } g(x = 3) or even yEqualsX <- function(f) function(x, y = x) f(x, y) yEqualsX(f)(x = 3) These are a lot like currying, but aren't currying, so they may be acceptable to you.  Personally I'd choose the f

[R] chaining closure arguments on-the-fly

2020-06-20 Thread Benjamin Tyner
Greetings, Occasionally, I desire to call a function with one argument set to equal to another. Here is a toy example:    f <- function(x, y) {        x + y    }    f(x = 3, y = x) # Error in f(x = 3, y = x) : object 'x' not found So far, the most concise way I found to accomplish this is:

[R] [R-pkgs] forestError 0.1.0: Random forest conditional prediction intervals, squared errors, and biases

2020-01-21 Thread Benjamin Lu
Hi all, I’m writing to introduce a new package, forestError. This package estimates conditional response quantiles (prediction intervals), conditional mean squared prediction errors, conditional biases, and conditional prediction error distribution functions for random forests using new methods pr

[R] safest way to subset (and replace) a data.table

2019-12-07 Thread Benjamin Tyner
Hi I would like to replace a data.table 'DT' with a subset of itself, where the subset is determined by an expression 'expr' which evaluates to logical. Thus far I've been using:    DT <- DT[expr, ] however this frequently results in a SIGABRT from glibc of the form:    "double free or corr

Re: [R] [SPAM] Re: The "--slave" option

2019-09-23 Thread Benjamin Lang
Hi Abby, I don’t really understand why you’re upset with me, but a) they’re cultured cell lines, not animals, b) they might cure people, c) I don’t do experiments, d) modern slavery, dated today: https://www.theguardian.com/world/2019/sep/21/such-brutality-tricked-into-slavery-in-the-thai-fishi

Re: [R] [SPAM] Re: The "--slave" option

2019-09-23 Thread Benjamin Lang
Hah, fair. I do hope somebody does see it and gives it a thought. Thanks, Ben On Sun, 22 Sep 2019 at 01:29, Roy Mendelssohn - NOAA Federal < roy.mendelss...@noaa.gov> wrote: > Please All: > > While as I said in my first post I am still not convinced that the OP was > in good faith to improve R a

Re: [R] [SPAM] Re: The "--slave" option

2019-09-20 Thread Benjamin Lang
> <https://en.wikipedia.org/wiki/Political_correctness> in 2004, and named > it the most politically > incorrect term of that year." > > The one thing "slave" does not mean in technology is any kind of human > being. > > On Thu, 19 Sep 2019 at 21:51, Benjami

Re: [R] [SPAM] Re: The "--slave" option

2019-09-20 Thread Benjamin Lang
s are > totally aside from what I assume is your main point, and because of this my > first reaction was don't feed the trolls. > > My $0.02. > > -Roy > >> On Sep 19, 2019, at 2:51 AM, Benjamin Lang wrote: >> >> Dear Richard, >> &g

Re: [R] [SPAM] Re: The "--slave" option

2019-09-19 Thread Benjamin Lang
Dear Richard, Thank you, that’s interesting. There is also something called an “etymological fallacy”. I think current usage is more useful here than the “science of truth”, i.e. the Ancient Greek idea that the (sometimes inferred) derivation of a word allows us to grasp “the truth of it”. In

[R] The "--slave" option

2019-09-18 Thread Benjamin Lang
he historically inclined, it does not make much of a difference whether the term evokes the Roman, Greek, American or modern kind of slavery for you: it is as disgusting as it gets. Thank you, Ben -- Benjamin Lang, PhD http://orcid.org/-0001-6358-8380 Marie Sklodowska-Curie Postdoctoral Fell

[R] Matrix::bdiag doesn't like being given a single named argument

2019-05-31 Thread Benjamin Tyner
Hello, Perhaps not a bug, but interesting because the error only happens when there is a single named argument.    > m <- matrix(1, 1, 1)    > library(Matrix)    > bdiag(m)    1 x 1 sparse Matrix of class "dgCMatrix"    [1,] 1    > bdiag(a = m)    Error in is.list(...) : supplied argument nam

Re: [R] undo compile? (or: remove bytecode from closure)

2018-07-16 Thread Benjamin Tyner
pe this helps, Rui Barradas Às 03:25 de 16-07-2018, Benjamin Tyner escreveu: Hi Given a closure which has been compiled, what's the recommended way to recover the original? For example,       > f <- function(x) x+1       > fc <- cmpfun(f)       > rm(f)       > fc       funct

[R] undo compile? (or: remove bytecode from closure)

2018-07-15 Thread Benjamin Tyner
Hi Given a closure which has been compiled, what's the recommended way to recover the original? For example,     > f <- function(x) x+1     > fc <- cmpfun(f)     > rm(f)     > fc     function(x) x+1     what's the best way to recover f from fc ? Regards Ben ___

Re: [R] Efficient manipulation with list object

2018-06-10 Thread Benjamin Christoffersen
You may be able to speed it up further by using `data.table`'s `rbindlist` or a similar function as shown here https://stackoverflow.com/a/49772719/5861244. 2018-06-10 21:20 GMT+02:00 Christofer Bogaso : > Using do.call() reduces my calculation time significantly. > > On Sun, Jun 10, 2018 at 10:45

[R] system() or pipe(..., open = "r") without child process?

2018-05-22 Thread Benjamin Tyner
Greetings On linux, is it possible to invoke an OS command from within R without spawning a child process? If not, is it possible to avoid copying the "parts of the caller's context" that are mentioned on the clone manpage?  ENOMEM Cannot  allocate  sufficient memory to allocate a task struct

Re: [R] rJava garbage collect

2018-02-06 Thread Benjamin Tyner
esolve the issue. Regards Ben On 6 February 2018 at 04:34, Benjamin Tyner https://stat.ethz.ch/mailman/listinfo/r-help>> wrote: >/Hi />//>/Does rJava offer a way to instruct the JVM to perform a garbage col

Re: [R] rJava garbage collect

2018-02-06 Thread Benjamin Tyner
Thanks Jeff; indeed it works:    .jcall("java/lang/System", method = "gc") On 02/05/2018 11:53 PM, Jeff Newmiller wrote: rJava offers a mechanism to call arbitrary methods in Java. Wouldn't you use that mechanism to call whatever you would call if you were programming in Java (e.g. System.gc

[R] rJava garbage collect

2018-02-05 Thread Benjamin Tyner
Hi Does rJava offer a way to instruct the JVM to perform a garbage collection? Regards Ben __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

[R] setSessionLimit

2018-01-19 Thread Benjamin Tyner
Hello The doc/NEWS.2 file mentions a setSessionLimit function, added with version 2.8.0     o   setTimeLimit() function to set limits on the CPU     and/or elapsed time for each top-level computation, and     setSessionLimit() to set limits for the rest of the session. However, I no l

[R] trying to find the multiple combinations...

2017-12-07 Thread Benjamin Sabatini
Hi, I'm trying to find a way to determine what multiples of the combination of three or more numbers equals a forth number. So, if I had a number set like: c(13.4689, 12.85212, 17.05071) What combination and multiples of these numbers would average to 15.0078? (so, something that would tell m

Re: [R] read.table(..., header == FALSE, colClasses = )

2017-10-24 Thread Benjamin Tyner
, colClasses = c(y = "character", x = "numeric"), col.names = c("foo", "bar"))      foo  bar    1   a 3.14 so I agree it is good that we are checking for that now. Regards Ben On 10/24/2017 08:55 AM, Martin Maechler wrote: Benjamin Tyner

Re: [R] read.table(..., header == FALSE, colClasses = )

2017-10-24 Thread Benjamin Tyner
Jeff, Thank you for your reply. The intent was to construct a minimum reproducible example. The same warning occurs when the 'file' argument points to a file on disk with a million lines. But you are correct, my example was slightly malformed and in fact gives an error under R version 3.2.2.

[R] read.table(..., header == FALSE, colClasses = )

2017-10-23 Thread Benjamin Tyner
Hello I noticed that starting with R version 3.3.0 onward, this generates a warning:    > txt <- c("a", "3.14")    > read.table(file = textConnection(txt), header = FALSE, colClasses = c(x = "character", y = "numeric")) the warning is "not all columns named in 'colClasses' exist" and I gues

Re: [R] Creating New Variable Using Ifelse

2017-08-10 Thread Courtney Benjamin
Thanks very much; with your tips, I was able to get the nested ifelse statement to work properly! Courtney Benjamin From: PIKAL Petr Sent: Thursday, August 10, 2017 5:39 AM To: Courtney Benjamin; r-help@r-project.org Subject: RE: Creating New Variable

[R] Creating New Variable Using Ifelse

2017-08-09 Thread Courtney Benjamin
oach at trying to handle the NAs properly-returns an error elsq2wbl$NONWBLRELJOB <- ifelse(elsq2wbl$PSWBL=="No" & elsq2wbl$EVERRELJOB=="Yes",1,ifelse(is.na(elsq2wbl$PSWBL)&is.na(elsq2wbl$EVERRELJOB),NA, ifelse(elsq2wbl$PSWBL!="No" & elsq2wbl$EVER

Re: [R] Svyglm Error

2017-07-05 Thread Courtney Benjamin
?Problem solved; I did not have the most updated version of R. When I updated it and ran the code, all worked well again. Thanks very much for helping me solve the mistake of an R novice! Courtney Benjamin Broome-Tioga BOCES Automotive Technology II Teacher Located at Gault Toyota

[R] Svyglm Error

2017-07-05 Thread Courtney Benjamin
a.action=na.omit) summary(Othpared)? Any help in resolving this concern would be greatly appreciated. Sincerely, Courtney Courtney Benjamin Broome-Tioga BOCES Automotive Technology II Teacher Located at Gault Toyota Doctoral Candidate-Educational Theory & Practice State University of Ne

[R] plot command error message

2017-06-04 Thread Benjamin Gföhler
I tried to plot a clustered linear regression model with the cplot command in R (code below). Leaflet is a binary variable (I know logit would be better), partisan is nummeric variable (0-4) and partisan_mis a dummy (0,1). As you can see it is clustered around two variables: around individuals

[R] Unknown anomaly

2017-04-18 Thread Benjamin Robira
ues had an answer. I hope that you would be able to help me fix that as it must be a pretty straightforward error that I do not realise. I would be thankful for any help, With my very Best Regards, Benjamin. __ R-help@r-project.org mailing list --

[R] flushing on.exit prior to q()

2017-02-21 Thread Benjamin Tyner
Hi, When using a custom error function that calls q(), what is the recommended way to "flush" the calling function's on.exit ? For example, say I have a script: #!/usr/bin/Rscript --no-init-file options(error = function() { cat("on error message\n", file = stderr()) q(sav

Re: [R] [FORGED] file.exists() on device files

2017-01-13 Thread Benjamin Tyner
Thank you for the insights, Rolf and Henrik. To give another example, this time in non-interactive mode, Rscript -e "file.exists(commandArgs(TRUE))" <(echo "Hi") [1] TRUE versus Rscript -e "normalizePath(commandArgs(TRUE))" <(echo "Hi") [1] "/dev/fd/63" Warning message: In n

[R] file.exists() on device files

2017-01-11 Thread Benjamin Tyner
Hi, On my linux machine (Ubuntu, and also tested on RHEL), I am curious to know what might be causing file.exists (and also normalizePath) to not see the final device file here: > list.files("/dev/fd", full.names = TRUE) [1] "/dev/fd/0" "/dev/fd/1" "/dev/fd/2" "/dev/fd/3" > file.exis

Re: [R] modify the imported version of a function

2016-12-19 Thread Benjamin Tyner
= FALSE, dotall = FALSE, ...) { options <- stri_opts_regex(case_insensitive = ignore_case, multiline = multiline, comments = comments, dotall = dotall, ...) structure(pattern, options = options, class = c("regex", "patt

[R] modify the imported version of a function

2016-12-16 Thread Benjamin Tyner
Hi I saw on the assignInNamespace help page, that it changes "the copy in the namespace, but not any copies already exported from the namespace, in particular an object of that name in the package (if already attached) and any copies already imported into other namespaces." So now I'm wonder

[R] Svyolr-Properly Specifying Start Values

2016-11-29 Thread Courtney Benjamin
[,1:16], repweights = elsq1ch[,18:217], weights = elsq1ch[,17], combined.weights = TRUE, type = "BRR") elsq1ch_brr allCColr <- svyolr(F3ATTAINMENT~F1PARED+BYINCOME+F1RACE+F1SEX+F1RGPP2+F1HIMATH+F1RTRCC,design=subset(elsq1ch_brr,BYSCTRL==1&G10COHRT==1),na.action=na.omit) ?

[R] Further Subsetting of Data for Log. Reg. Results in qr.default Error

2016-11-23 Thread Courtney Benjamin
is resulting in the qr.default error CTE <- svyglm(formula=F3ATTAINB~F1PARED+BYINCOME+F1RACE+F1SEX+F1RGPP2+F1HIMATH,family="binomial",design=elsq1ch_brr,subset=BYSCTRL==1&G10COHRT==1&F1RTRCC=="Academic",na.action=na.omit) summary(CTE) Courtney Benjamin Broome-Tioga

Re: [R] Archer-Lemeshow Goodness of Fit Test for Survey Data with Log. Regression

2016-11-18 Thread Courtney Benjamin
-svrepdesign(variables = elsq1ch, repweights = elsq1ch[,18:217], weights = elsq1ch[,17], combined.weights = TRUE, type = "BRR") decilemodel<- svyglm(r~your_g, design=newdesign,subset=BYSCTRL==1&G10COHRT==1) regTermTest(decilemodel, ~your_g) ​ ​ Courtney Benjamin Broome-

Re: [R] Archer-Lemeshow Goodness of Fit Test for Survey Data with Log. Regression

2016-11-18 Thread Courtney Benjamin
​Thank you; I appreciate your advisement. Courtney Benjamin Broome-Tioga BOCES Automotive Technology II Teacher Located at Gault Toyota Doctoral Candidate-Educational Theory & Practice State University of New York at Binghamton cbenj...@btboces.org<mailto:cbenj...@btboces.org>

Re: [R] Archer-Lemeshow Goodness of Fit Test for Survey Data with Log. Regression

2016-11-18 Thread Courtney Benjamin
indicate a poor model fit. Sincerely, Courtney Courtney Benjamin Broome-Tioga BOCES Automotive Technology II Teacher Located at Gault Toyota Doctoral Candidate-Educational Theory & Practice State University of New York at Binghamton cbenj...@btboces.org<mailto:cbenj...@btboces.or

[R] Archer-Lemeshow Goodness of Fit Test for Survey Data with Log. Regression

2016-11-16 Thread Courtney Benjamin
naming/specifying the new variables of r and g properly transform(elsq1ch,r=r,g=g) elsq1ch_brr <- update(elsq1ch_brr,tag=g,tag=r) #then: decilemodel<- svyglm(r~g, design=newdesign) regTermTest(decilemodel, ~g) #is the F-adjusted mean residual test from the Archer Lemeshow paper Thank you, Courtney

Re: [R] Resetting Baseline Level of Predictor in svyglm Function

2016-11-03 Thread Courtney Benjamin
Thank you, Anthony; it worked flawlessly.? Courtney Benjamin Broome-Tioga BOCES Automotive Technology II Teacher Located at Gault Toyota Doctoral Candidate-Educational Theory & Practice State University of New York at Binghamton cbenj...@btboces.org<mailto:cbenj...@btboces.org>

[R] Resetting Baseline Level of Predictor in svyglm Function

2016-10-31 Thread Courtney Benjamin
for the predictors allCC <- svyglm(formula=F3ATTAINB~F1PARED+BYINCOME+F1RACE+F1SEX+F1RGPP2+F1HIMATH+F1RTRCC,family="binomial",design=elsq1ch_brr,subset=BYSCTRL==1&G10COHRT==1,na.action=na.omit) summary(allCC) Any guidance is greatly appreciated.? Sincerely, Courtney? Courtney B

Re: [R] Significance of Svyrepdesign Object Warning

2016-10-27 Thread Courtney Benjamin
​Thank you; I will do so. Courtney Benjamin Broome-Tioga BOCES Automotive Technology II Teacher Located at Gault Toyota Doctoral Candidate-Educational Theory & Practice State University of New York at Binghamton cbenj...@btboces.org<mailto:cbenj...@btboces.org> 60

Re: [R] Significance of Svyrepdesign Object Warning

2016-10-27 Thread Courtney Benjamin
") elsq1ch_brr #Logistic regression call which yields a warning regarding svyrepdesign object allCC <-svyglm(formula=F3ATTAINB~F1PARED+BYINCOME+F1RACE+F1SEX+F1RGPP2+F1HIMATH+F1RTRCC,family="binomial",design=elsq1ch_brr,subset=BYSCTRL==1&G10COHRT==1,na.action=na.exclude) summar

Re: [R] Significance of Svyrepdesign Object Warning

2016-10-23 Thread Courtney Benjamin
el of two of my predictor variables (BYINCOME & F1HIMATH) were changed. Courtney Benjamin Broome-Tioga BOCES Automotive Technology II Teacher Located at Gault Toyota Doctoral Candidate-Educational Theory & Practice State University of New York at Binghamton cbenj...@btboces.org<mailto

[R] Significance of Svyrepdesign Object Warning

2016-10-23 Thread Courtney Benjamin
es.1252 #attached base packages: # [1] grid stats graphics grDevices utils datasets methods base #other attached packages: #[1] survey_3.31-2 survival_2.39-4 Matrix_1.2-6RCurl_1.95-4.8 bitops_1.0-6 #loaded via a namespace (and not attached): #[1] tools_3.3.1 splines

Re: [R] Svyglm Error in Survey Package

2016-09-25 Thread Courtney Benjamin
To: Courtney Benjamin Cc: r-help@r-project.org Subject: RE: Svyglm Error in Survey Package Dear Courtney, You're confusing a function call, na.action(na.omit), with an argument specification, na.action=na.omit (and, in any event, there is no na.action() function). But you don't have

[R] Svyglm Error in Survey Package

2016-09-23 Thread Courtney Benjamin
TRL==1 & G10COHRT==1) dim(samp1) mc1 <- svyglm(F3ATTAINMENT~F1SES2QU+F1RGPP2,elsq1ch_brr,subset=samp1,na.action) summary(mc1)? Both attempts resulted in the same error stated above. Any advisement in how to resolve this error would be greatly appreciated. Sincerely, Courtney Benjamin ?

[R] commercial license

2016-09-12 Thread benjamin . stocker
Dear r-project Team How does It cost a commercial license for the R Console and the R-comander GUI without the Rstudio enviroment. Thanks for helping me. Freundliche Gr�sse/Kind regards Benjamin Stocker Reporting/Controlling MBC Mercedes-Benz Schweiz AG Bernstrasse 55 8952 Schlieren

[R] How to Display Value Labels in R Outputs?

2016-07-28 Thread Courtney Benjamin
(file="els-Q1-04-21-16.dat", header = TRUE, sep = "\t", quote = "\"", dec =".") or I am not specifying some detail that is required to manually assign labels to the values of the variables.? Courtney Benjamin Broome-Tioga BOCES Automotiv

[R] R_DirtyImage and Rprof

2016-05-11 Thread Benjamin Tyner
Hello, I have some code which was running in interactive mode while Rprof(..., line.profiling = TRUE). Near the end of my script, it opens up a pipe(..., open = "w") to a perl script, and at that point the execution gets stuck using 100% cpu. (The perl script itself never showed up in pstree

[R] party package: is cforest really bagging?

2015-12-15 Thread Benjamin
Hi all, I'm using the "party" package to create random forest of regression trees. I've created a ForestControl class in order to limit my number of trees (ntree), of nodes (maxdepth) and of variables I use to fit a tree (mtry). One thing I'm not sure of is if the cforest algo is using subsets of

Re: [R] choropleth packages (US)

2015-12-10 Thread Benjamin Tyner
ec 7, 2015 at 9:34 AM, Erich Neuwirth wrote: ggplot2 also can do this with fortify geom_polygon Von meinem iPad gesendet Am 06.12.2015 um 21:03 schrieb Benjamin Tyner : Hi I wish to draw a basic choropleth (US, by state) and am wondering if anyone has any recommendations? I've t

[R] choropleth packages (US)

2015-12-06 Thread Benjamin Tyner
Hi I wish to draw a basic choropleth (US, by state) and am wondering if anyone has any recommendations? I've tried the following thus far: 1. choroplethr: this works, but required installation of 30+ dependencies. I would prefer something with fewer dependencies. 2. tmap: this also seems prom

[R] efficiently multiply each column of a sparse Matrix by a sparse vector

2015-11-28 Thread Benjamin Tyner
Hi, Say I have a sparse Matrix X, and a sparse vector (stored as a 1-column sparse Matrix A), with X and A having the same number of rows, and I wish to multiply each column of X by A, but would like the operation to take full advantage of the sparseness of both X and A. In other words I want

Re: [R] backsolve, chol, Matrix, and SparseM

2015-09-25 Thread Benjamin Tyner
path in the first place. By the way, is R-forge still the correct place to report bugs in package:Matrix? Regards Ben On 09/25/2015 04:25 AM, Martin Maechler wrote: > Dear Ben, > >>>>>> Benjamin Tyner >>>>>> on Thu, 24 Sep 2015 13:47:58 -0400 writes:

[R] backsolve, chol, Matrix, and SparseM

2015-09-24 Thread Benjamin Tyner
Hi I have some code which does (on a symmetric matrix 'x') backsolve(chol(x), diag(nrow(x))) and I am wondering what is the recommended way to accomplish this when x is also sparse (from package:Matrix). I know that package:Matrix provides a chol method for such matrices, but not a backsolve

[R] Rprof and system

2015-08-11 Thread Benjamin Tyner
Hi I have an R script which invokes WriteXLS() (from the package of the same name) which as you may know, calls perl via system(). I've noticed that when I enable profiling using Rprof(), when the script gets to the part where perl is called, it gets "stuck": it just sits there using 99-100% CPU a

[R] Combining multiple probability weights for the sample() function.

2015-06-02 Thread Benjamin Ward (ENV)
Dear R-List, I have a set of possibilities I want to sample from: bases <- list(c('A', 'C'), c('A', 'G'), c('C', 'T')) possibilities <- as.matrix(expand.grid(bases)) >possibilities Var1 Var2 Var3 [1,] "A" "A" "C" [2,] "C" "A" "C" [3,] "A" "G" "C" [4,] "C" "G" "C" [5,] "A" "A" "T" [6,]

[R] Scatterplot : smoothing colors according to density of points

2015-06-02 Thread Benjamin Dubreuil
Hello everyone, I have a data frame D with 4 columns id,X,Y,C. I want to plot a simple scatter plot of D$X vs. D$Y and using D$C values as a color. (id is just a text string not used for the plot) But actually, I don't want to use the raw values of D$C, I would prefer to calculate the average v

[R] specific package to laply

2015-05-14 Thread Benjamin
fun(obj, substitute(ex), parent.frame(), e$data) : worker initialization failed: there is no package called ‘rvest’ Calls: %dopar% -> Presumably, I need a way to export my .libPaths() to the nodes. Any suggestions? Thanks, Benjamin [[alternative HTML version deleted]] _

Re: [R] package implementing continuous binomial?

2015-05-07 Thread Benjamin Tyner
Thanks David! I'll take a look at zipfR. Regards Ben On 05/07/2015 03:10 PM, David Winsemius wrote: > On May 6, 2015, at 7:00 PM, Benjamin Tyner wrote: > >> Hi >> >> I'm wondering if anyone is aware of an R package implementing (i.e., >> providing a pdf

[R] package implementing continuous binomial?

2015-05-06 Thread Benjamin Tyner
Hi I'm wondering if anyone is aware of an R package implementing (i.e., providing a pdf, cdf, and/or quantile function) for the continuous binomial distribution? Specifically the one characterized here: http://www2.math.uni-paderborn.de/fileadmin/Mathematik/AG-Indlekofer/Workshop/Satellite_meetin

Re: [R] non-terminal token lacking children from utils::getParseData

2015-05-02 Thread Benjamin Tyner
Thank you Yihui for also reporting the bug here: https://bugs.r-project.org/bugzilla/show_bug.cgi?id=16354 and thank you Duncan for finding the issue and fixing it! I definitely like your idea to report a summary message instead of the long text string. Regards Ben > I tried to reduce the offen

Re: [R] Having trouble with gdata read in

2015-03-27 Thread Benjamin Baker
y.ComAMI                     : num  0 NA NA NA 1600 NA NA NA NA NA ... 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, Mar 25, 2015 at 5:01 PM, Benjamin Baker wrote: Trying to rea

Re: [R] Having trouble with gdata read in

2015-03-27 Thread Benjamin Baker
: num 0 0 0 155 0 0 0 0 0 0 ... > $ Q8.Hourly.ResAMI : num 0 NA NA NA 16100 NA NA NA NA > NA ... > $ Q8.Hourly.ComAMI : num 0 NA NA NA 1600 NA NA NA NA > NA ... > > Jim Holtman > Data Munger Guru > What is the proble

Re: [R] Having trouble with gdata read in

2015-03-27 Thread Benjamin Baker
t;- tempfile() > ami <- " > http://www.ferc.gov/industries/electric/indus-act/demand-response/2008/survey/ami_survey_responses.xls > " > download.file( ami , tf , mode = 'wb' ) > ami.data2008 <- read.xlsx( tf , sheetIndex = 1 ) > On Wed, Mar 25, 2015 at 5:0

[R] Having trouble with gdata read in

2015-03-25 Thread Benjamin Baker
I can find some “solutions” for read.table, but nothing specific to read.xls Many thanks, Benjamin Baker — Sent from Mailbox [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, se

Re: [R] mboost: Proportional odds boosting model - how to specify the offset?

2015-03-20 Thread Benjamin Hofner
). More details are given in the example below. Best, Benjamin [1] http://r-forge.r-project.org/projects/mboost/ Example code library(MASS) library(mboost) data(iris) iris$Species <- factor(iris$Species, ordered = T) p.iris <- polr(Species ~ Sepal.Length, data = iris) p.ir

Re: [R] non-terminal token lacking children from utils::getParseData

2015-02-19 Thread Benjamin Tyner
I tried to reduce the offending portion as best I could to a more-or-less minimal example (1136 bytes), which can be downloaded via: wget https://www.dropbox.com/s/74rgxr5x2aalr99/badstring.R then once in R, > b <- parse(file = "~/badstring.R", keep.source = TRUE) > d <- getParseData

[R] quadprog::solve.QP sometimes returns NaNs

2014-09-26 Thread Benjamin Tyner
Hello, Here is an example; hopefully it is reproducible on others' platform: library(quadprog) n <- 66L set.seed(6860) X <- matrix(1e-20, n, n) diag(X) <- 1 Dmat <- crossprod(X) y <- seq_len(n) dvec <- crossprod(X, y) Amat <- diag(n) bvec <- y + runif(n)

[R] scaling of Amat still recommended for quadprog::solve.QP ?

2014-09-25 Thread Benjamin Tyner
Greetings, I ran across this recommendation, to keep the norms of the columns of the Amat on similar scales, https://stat.ethz.ch/pipermail/r-help/2007-September/141335.html However, when I looked at the code, I noticed that the norms are already being calculated: c c calculate the

Re: [R] from cut.Date

2014-09-17 Thread Benjamin Tyner
t; args(cut.Date) function (x, breaks, labels = NULL, start.on.monday = TRUE, right = FALSE, ...) I suppose the latter does make sense, given that days, months, years etc are right-continuous functions of time. Regards Ben > On 17/09/2014 12:04, Benjamin Tyner wrote: > >/ Hello, > />/

[R] from cut.Date

2014-09-17 Thread Benjamin Tyner
Hello, I'm wondering if this is expected? > cut(structure(1, class="Date"), structure(c(11100,1), class="Date")) [1] Levels: 2000-05-23 The help page says that "for ‘"Date"’ objects, only ‘"day"’, ‘"week"’, ‘"month"’, ‘"quarter"’ and ‘"year"’ are allowed" [for the 'breaks' a

Re: [R] lattice: packet.number() versus panel.number()

2014-08-27 Thread Benjamin Tyner
On 08/27/2014 02:10 AM, David Winsemius wrote: > On Aug 26, 2014, at 3:25 PM, Benjamin Tyner wrote: > >> Hi, >> >> According to >> https://svn.r-project.org/R-packages/trunk/lattice/R/print.trellis.R, >> >>"[panel.number] is usually the same as, bu

[R] lattice: packet.number() versus panel.number()

2014-08-26 Thread Benjamin Tyner
Hi, According to https://svn.r-project.org/R-packages/trunk/lattice/R/print.trellis.R, "[panel.number] is usually the same as, but can be different from packet.number" and I had been under the impression that as long as the user is not using a custom index.cond nor perm.cond, the panel.numbe

Re: [R] non-interactive subset, `$`, and quotes

2014-06-26 Thread Benjamin Tyner
Duncan, How embarrassing! Thanks. -Ben On 06/26/2014 10:25 PM, Duncan Murdoch wrote: > On 27/06/2014, 4:08 AM, Benjamin Tyner wrote: >> Hi >> >> I know that subset() is not intended for use in programming. However I >> am still curious to learn why, in non-inter

[R] non-interactive subset, `$`, and quotes

2014-06-26 Thread Benjamin Tyner
Hi I know that subset() is not intended for use in programming. However I am still curious to learn why, in non-interactive mode, if I take away the quotes around 'bar' Rscript -e "foo <- list(bar = iris); head(subset(foo$'bar', Species=='setosa'))" Sepal.Length Sepal.Width Petal.Length P

Re: [R] detecting the sourcing of site profile on Startup versus post-Startup

2014-04-24 Thread Benjamin Tyner
- > Sent from my phone. Please excuse my brevity. > > On April 23, 2014 6:11:09 PM PDT, Benjamin Tyner wrote: >> Thanks Duncan! Yes, I considered taking advantage of .First, but was >> concerned that the .First defined by the site profile could be masked &

Re: [R] detecting the sourcing of site profile on Startup versus post-Startup

2014-04-23 Thread Benjamin Tyner
gave was a simplifying assumption, sorry about that). On 04/23/2014 06:55 AM, Duncan Murdoch wrote: On 22/04/2014, 8:59 PM, Benjamin Tyner wrote: Greetings, Is there any way to programmatically detect whether a piece of code is being run within the initial (Startup) sourcing o

[R] detecting the sourcing of site profile on Startup versus post-Startup

2014-04-22 Thread Benjamin Tyner
Greetings, Is there any way to programmatically detect whether a piece of code is being run within the initial (Startup) sourcing of the site profile? For example, say I have a site profile, "/path/to/Rprofile.site". Is there any function "my_func" which would return different values for these tw

Re: [R] package environment versus namespace environment

2014-03-09 Thread Benjamin Tyner
earch path, then functions in the base namespace would no longer be able to see objects in utils? (I realize the answer may be vacuous, since library() does not honor pos=1). Regards, Ben On 03/09/2014 09:09 AM, Duncan Murdoch wrote: On 14-03-08 6:42 PM, Benjamin Tyner wr

Re: [R] package environment versus namespace environment

2014-03-08 Thread Benjamin Tyner
es back in alphabetical order.) Regards Ben On 03/07/2014 11:46 AM, Duncan Murdoch wrote: On 07/03/2014 10:16 AM, Benjamin Tyner wrote: Hello, I realize that a function in (for example, function "head1" below) is unable to see (without resorting to "::", anyw

[R] package environment versus namespace environment

2014-03-07 Thread Benjamin Tyner
Hello, I realize that a function in (for example, function "head1" below) is unable to see (without resorting to "::", anyway) objects in utils (for example, "head" below), since package:base is after package:utils on the search path. However, I'm wondering what is the machinery that allows a fu

[R] An unbalanced factorial experiment - function for GLM model with quasi-Poisson distribution

2014-02-20 Thread Benjamin Jarcuska
like ANOVA's sums of squares type III) should I used in R? Thank you for your answers! Sincerely, Benjamin Jarcuska Inst Forest Ecol, Slovak Acad Sci Zvolen, Slovakia [[alternative HTML version deleted]] __ R-help@r-project.org ma

[R] Two x axes - top and bottom

2014-02-08 Thread Benjamin Ward (ENV)
Hi, fellow R users, I've been asked to make a plot with two datasets each with a different x axis, and it's been suggested one be at the top and the other at the bottom of the graph. I normally use ggplot2, and I know how to plot multiple datasets by simply + a new geom with a different data op

[R] Change points in R

2014-01-30 Thread Benjamin Ward (ENV)
Hi R helpers, I have a set of data best shown in this below graph. Each coloured line represents a statistic calculated across pairs of DNA sequences. And for each coloured line, I would like to identify breakpoints - so identify the chunks where the values are high, for example, in the light

Re: [R] APA two-way bar-chart

2014-01-04 Thread Benjamin Telkamp
ary(fun.y = mean, geom = "bar", position = "dodge") + stat_summary(fun.data = mean_cl_normal, geom = "errorbar", position = position_dodge(width=0.90), width = 0.2) + labs (x = "PT Manipulatie", y = "Percentage Balworpen") you would make me so happy! th

Re: [R] APA two-way bar-chart

2014-01-04 Thread Benjamin Telkamp
Hi Jim, thanks for your response, i really appreciate it, I'm attaching a picture of a two way bar chart made by excel (i hate excel and i want to be desperately able to do these things with R) hope you could give me some command lines that will help me, Thanks Again. Benjamin Telkamp

[R] APA two-way bar-chart

2014-01-03 Thread Benjamin Telkamp [via R]
y one give me a hint or solution? See Below for my commands thank you so much, Benjamin Telkamp (Email: inzichtel...@gmail.com, The Netherlands) these are the commands used so far which just didn't give me what i wanted library(ggplot2) library(foreign) install.packages("Hmisc"

Re: [R] Extract values from vector and repeat by group

2013-11-18 Thread Benjamin Gillespie
November 2013 16:48 To: R help Cc: Berend Hasselman; Benjamin Gillespie Subject: Re: [R] Extract values from vector and repeat by group Hi, ?merge() sometimes change the order. For example: df1 <- df[-12,] df2 <- df1 merge(df1, df1[df1$time == 4, c("group", "var")], by.x

[R] Extract values from vector and repeat by group

2013-11-17 Thread Benjamin Gillespie
Hi all, I hope you can help. I have a data frame 'df': group=c(rep(1,8),rep(2,10),rep(3,11)) var=rnorm(29) time=c(seq(1,8),seq(1,10),seq(1,11)) df=data.frame(group,var,time) I would like to extract the value from 'var' for each 'group' at 'time'=4 and repeat these extracted values in a new vec

Re: [R] Small p from binomial probability function.

2013-10-16 Thread Benjamin Ward (ENV)
in so 0.05, 0.50 and 0.95, or even just 0.05 and 0.95? Then uniroot finds the root of this function and doing so find me the small p I need? Best, Ben. From: Rolf Turner [rolf.tur...@vodafone.co.nz] Sent: 11 October 2013 02:11 To: Benjamin Ward (EN

  1   2   3   4   >