Displaying 20 results from an estimated 10000 matches similar to: "options(errorfn=traceback)"
2013 Feb 07
4
Hard Stop?
is it possible to throw a stop() that is so hard that it will escape
even tryCatch?
/iaw
----
Ivo Welch (ivo.welch at gmail.com)
2010 Aug 22
2
on abort error, always show call stack?
Dear R Wizards---is it possible to get R to show its current call
stack (sys.calls()) upon an error abort? I don't use ESS for
execution, and it is often not obvious how to locate how I triggered
an error in an R internal function. Seeing the call stack would make
this easier. (right now, I sprinkle "cat" statements everywhere, just
to locate the line where the error appears.) Of
2013 Feb 09
2
character strings with embedded commands: perl "/gee" ?
dear R experts---I am trying to replicate a perl feature. I want to be
able to embed R commands inside a character string, and have the
string be printed with the command executed. my perl equivalent is
my $a=10;
my $teststring = "the expression, $a+1, is ::$a+1::, but add one
more for ::$a+2::\n";
$teststring =~ s/::(.*?)::/$1/gee;
print $teststring;
of course, R does not use
2010 Nov 24
1
segfault interest?
in a long program, I ran into
?*** caught segfault ***
address 0xdc3f9b48, cause 'memory not mapped'
Traceback:
?1: rep.int(seq_len(nx), rep.int(rep.fac, nx))
?2: rep.int(rep.int(seq_len(nx), rep.int(rep.fac, nx)), orep)
?3: expand.grid(seq_len(nx), seq_len(ny))
?4: merge.data.frame(d, ss)
?5: merge(d, ss)
?6: valid.range(opt)
?7: eval.with.vis(expr, envir, enclos)
?8: eval.with.vis(ei,
2004 Aug 21
4
loadhistory() in .Rprofile ?
dear wizards: my .Rprofile has just one command for testing,
loadhistory("~/.Rhistory")
but this gives me an error on R startup:
Error: couldn't find function "loadhistory"
Invoking loadhistory() as the first interactive command works fine;
incidentally, I believe loadhistory() in the .Rprofile worked in
earlier or other platform R releases, too.
Is the .Rprofile
2013 Feb 26
1
Light Libraries
Dear R-Help group:
I have been tinkering with how I want my personal standard library
functions to look like. They are not designed to be professional and
heavyweight, but lightweight. There are probably dozens of little bugs,
because I don't know or have not properly taken care of a variety of
internal R code issues. still, I like how this ended up, and there is no
learning curve, so I
2013 Feb 06
5
First R Package --- Advice?
Dear R experts---
after many years, I am planning to give in and write my first R
package. I want to combine my collection of collected useful utility
routines.
as my guide, I am planning to use Friedrich Leisch's "Creating R
Packages: A Tutorial" from Sep 2009. Is there a newer or better
tutorial? this one is 4 years old.
I also plan on one change---given that the
2012 Mar 30
4
list assignment syntax?
Dear R wizards: is there a clean way to assign to elements in a list?
what I would like to do, in pseudo R+perl notation is
f <- function(a,b) list(a+b,a-b)
(c,d) <- f(1,2)
and have c be assigned 1+2 and d be assigned 1-2. right now, I use the clunky
x <- f(1,2)
c <- x[[1]]
d <- x[[2]]
rm(x)
which seems awful. is there a nicer syntax?
regards, /iaw
----
Ivo Welch
2011 Oct 10
5
multicore by(), like mclapply?
dear r experts---Is there a multicore equivalent of by(), just like
mclapply() is the multicore equivalent of lapply()?
if not, is there a fast way to convert a data.table into a list based
on a column that lapply and mclapply can consume?
advice appreciated...as always.
regards,
/iaw
----
Ivo Welch (ivo.welch at gmail.com)
2019 Apr 05
2
Deep Replicable Bug With AMD Threadripper MultiCore
The following program is whittled down from a much larger program that
always works on Intel, and always works on AMD's threadripper with
lapply but not mclappy. With mclapply on AMD, all processes go into
"suspend" mode and the program then hangs. This bug is replicable on an
AMD Ryzen Threadripper 2950X 16-Core Processor (128GB RAM), running
latest ubuntu 18.04. The R version
2011 Jul 02
5
%dopar% parallel processing experiment
dear R experts---
I am experimenting with multicore processing, so far with pretty
disappointing results. Here is my simple example:
A <- 100000
randvalues <- abs(rnorm(A))
minfn <- function( x, i ) { log(abs(x))+x^3+i/A+randvalues[i] } ?## an
arbitrary function
ARGV <- commandArgs(trailingOnly=TRUE)
if (ARGV[1] == "do-onecore") {
?library(foreach)
?discard <-
2010 Jan 08
4
fast lm se?
dear R experts---I am using the coef() function to pick off the coefficients
from an lm() object. alas, I also need the standard errors and I need them
fast. I know I can do a "summary()" on the object and pick them off this
way, but this computes other stuff I do not need. Or, I can compute (X'
X)^(-1) s^2 myself. Has someone written a fast se() function?
incidentally, I think
2013 Apr 04
6
categorized complete list of R commands?
every time I read the R release notes for the next release, I see many
functions that I had forgotten about and many functions that I never knew
existed to begin with. (who knew there were bibtex facilities in R?
obviously, everyone except me.)
I wonder whether there is a complete list of all R commands (incl the
standard packages) somewhere, preferably each with its one-liner AND
2010 Jun 11
3
lm without error
this is not an important question, but I wonder why lm returns an
error, and whether this can be shut off. it would seem to me that
returning NA's would make more sense in some cases---after all, the
problem is clearly that coefficients cannot be computed.
I know that I can trap the lm.fit() error---although I have always
found this to be quite inconvenient---and this is easy if I have only
2010 Aug 30
4
different interface to by (tapply)?
dear R experts:
has someone written a function that returns the results of by() as a
data frame? ??of course, this can work only if the output of the
function that is an argument to by() is a numerical vector.
presumably, what is now names(byobject) would become a column in the
data frame, and the by object's list elements would become columns.
it's a little bit like flattening the by()
2012 May 31
2
print.data.frame to string?
dear R experts---is there a function that prints a data frame to a string?
cat() cannot handle lists, so I cannot write cat("your data frame is:\n",
df, "\n").
regards, /iaw
----
Ivo Welch (ivo.welch@gmail.com)
[[alternative HTML version deleted]]
2012 May 09
2
big quasi-fixed effects OLS model
dear R experts---now I have a case where I want to estimate very large
regression models with many fixed effects---not just the mean type, but
cross-fixed effects---years, months, locations, firms. Many millions of
observations, a few thousand variables (most of these variables are
interaction fixed effects). could someone please point me to packages, if
any, that would help me estimate such
2011 May 15
4
"Low Pain" Unicode Characters in pdf graph?
Dear R-experts---is there a relatively low-pain way to get unicode
characters into a plot to a pdf device?
pdf(file="cardsymbols.pdf")
plot( 0, xlim=c(0,5), ylim=c(0,5), type="n")
text(1,1, "♠")
text(2,2, "♥")
text(3,3, "♦")
text(4,4, "♣")
dev.off()
(these are the characters that I need the most
2004 Jul 07
3
fast NA elimination ?
dear R wizards: an operation I execute often is the deletion of all
observations (in a matrix or data set) that have at least one NA. (I
now need this operation for kde2d, because its internal quantile call
complains; could this be considered a buglet?) usually, my data sets
are small enough for speed not to matter, and there I do not care
whether my method is pretty inefficient (ok, I
2013 Jan 03
6
Bounty on Error Checking
Dear R developers---I just spent half a day debugging an R program,
which had two bugs---I selected the wrongly named variable, which
turns out to have been a scalar, which then happily multiplied as if
it was a matrix; and another wrongly named variable from a data frame,
that triggered no error when used as a[["name"]] or a$name . there
should be an option to turn on that throws an