Hi,
On Mar 29, 2007, at 8:55 AM, Erich Neuwirth wrote:
> It is accessible at http://dx.doi.org/10.1007/s00180-007-0023-6
> If you cannot access it, write to me, I will send you a copy.
>
> Currently, we are working on a cross-platform alternative to using
> COM to connect the spread-sheet to R. Th
Hi,
I would like to potentially become involved with this project. I am
still learning about R-project internals but I have been digging
around a lot in the internal R-code. That said, I am initimately
familiar with OOo since I was involved as a volunteer for that
project previously. I
Hi,
From an svn update run just a few minutes ago, make check-all fails
running code in 'R-intro.R' ... OK
comparing 'R-intro.Rout' to '../../r-devel/r-devel/R/tests/R-
intro.Rout.save' ...608c608
< alternative hypothesis: true location shift is not equal to 0
---
> alternative hypothesis: tru
Hi,
iff is shorthand typically used in mathematical proofs. You have
probably found the one group (users of R) that have the highest
probability of knowing the meaning of "iff". That said, I agree that
it should not be used in official documentation. It should be
replaced with "if and o
Hi,
One final question ... is there an automated build system anyplace
for package developers and R -developers to make sure their code
works on all platforms. I was able to check with Linux, and MacOSX
(both Intel and PPC) but I do not own a Windows box and was unable to
test/debug the b
Hi,
> If you unsuccessfully attempted to upload a file once but got
> connected, it will prevent you from trying it again as the file
> already exists.
Ah! That must be it. I never got a successful completion message,
just a hang with a message about going into PASSIVE mode. So some or
Hi,
Thanks anyway but I can connect as well and was even able to cd into
incoming but it would not allow me upload any files. That is what
the debug log showed below.
Can you upload things into incoming? Perhaps there is simply
something funny about my ftp client? I will try from my MacO
ate file.
ftp>
Any ideas?
Kevin
On Aug 3, 2006, at 6:46 PM, Kevin B. Hendricks wrote:
> Hi,
>
> Thanks to lots of help and hints, I now have a working validated
> external version of all of the igroup functions. The source Package
> is called Rigroup_0.80.0.tar.gz and it pass
Hi,
Thanks to lots of help and hints, I now have a working validated
external version of all of the igroup functions. The source Package
is called Rigroup_0.80.0.tar.gz and it passes R CMD check of my Linux
box and will install and pass all tests on my Mac OSX(Intel and PPC)
machines and
Hi,
My last word on this topic until I get a working external R package ...
The igroup code has now been validated both with and without NAs and
with and without removing them. Thanks to Bill, Tom, Thomas, and
everyone for your helpful comments and hints.
The results for my validation run a
Hi Tom,
> Whether or not the R core developers want to merge these functions
> in base
> R, they would make a great little package on CRAN. That way others
> could
> easily use them, and for yourself, the package automatically gets
> updated
> with new versions of R. It sounds like you're do
Hi Thomas,
Here is a comparison of performance times from my own igroupSums
versus using split and rowsum:
> x <- rnorm(2e6)
> i <- rep(1:1e6,2)
>
> unix.time(suma <- unlist(lapply(split(x,i),sum)))
[1] 8.188 0.076 8.263 0.000 0.000
>
> names(suma)<- NULL
>
> unix.time(sumb <- igroupSum
Hi Tom,
> Now, try sorting and using a loop:
>
>> idx <- order(i)
>> xs <- x[idx]
>> is <- i[idx]
>> res <- array(NA, 1e6)
>> idx <- which(diff(is) > 0)
>> startidx <- c(1, idx+1)
>> endidx <- c(idx, length(xs))
>> f1 <- function(x, startidx, endidx, FUN = sum) {
> + for (j in 1:length(res)) {
Hi Bill,
After playing with this some more and adding an implementation to
handle NAs in the data vector, I have run into the problem of what to
return when the only data values for a particular bin (or level) in
the data vector were NAs and the user selected na.rm=T
1. Should it return 0 f
Hi Bill,
So you wrote one routine that can calculate any single of a variety
of stats and allows weights, is that right? Can it return a data
frame of any subset of requested stats as well (that is what I was
thinking of doing anyway).
I think someone can easily calculate all of those thin
Hi Bill,
>>>sum : igroupSums
Okay, after thinking about this ...
# assumes i is the small integer factor with n levels
# v is some long vector
# no sorting required
igroupSums <- function(v,i) {
sums <- rep(0,max(i))
for (j in 1:length(v)) {
sums[[i[[j <- sums[[i[[j + v
Hi Bill,
> Splus8.0 has something like what you are talking about
> that provides a fast way to compute
> sapply(split(xVector, integerGroupCode), summaryFunction)
> for some common summary functions. The 'integerGroupCode'
> is typically the codes from a factor, but you could compute
> it in
Hi,
> There was a performance comparison of several moving average
> approaches here:
> http://tolstoy.newcastle.edu.au/R/help/04/10/5161.html
>
Thanks for that link. It is not quite the same thing but is very
similar.
The idea is to somehow make functions that work well over small sub-
sequ
Hi,
I was using my installed R which is 2.3.1 for the first tests. I
moved to the r-devel tree (I svn up and rebuild everyday) for my "by"
tests to see if it would work any better. I neglected to retest
"merge" with the devel version.
So it appears "merge" is already fixed and I just need
Hi Developers,
I am looking for another new project to help me get more up to speed
on R and to learn something outside of R internals. One recent R
issue I have run into is finding a fast implementations of the
equivalent to the following SAS code:
/* MDPC is an integer sort key made fro
Hi,
Is there anybody interested in this patch to fix bug 8141 - the
rewriting substituteList in coerce.c?
If so, I have updated it to today's tree. Testers welcome as well.
If not, please let me know.
Thanks,
Kevin
Attached is both the patch and a gzipped version of the patch that
is t
Hi,
Attached is both the patch and a gzipped version of the patch that is
the fix for bug 8141 - rewriting substituteList in coerce.c to use a
loop instead of recursion to walk the list.
The new version passes all of my tests (make check-all, etc) and I
have used it with no negative impac
Hi,
Just in case anyone is interested in bug 8141. With some further
testing, I found and fixed a small bug in my code and simplified it
a bit.
Here is the latest version of the coerce.c routine substituteList
that fixes the C stack overflow problem as documented in Bug 8141.
Comments w
Hi,
I have a proposed fix for Bug 8141 that passes make check-all on my
machine and that will actually NOT overflow the C stack even for the
larger problems than the test case given in 8141.
The basic idea is to not use recursion to walk the list elements and
instead use a loop building up
Hi,
Just in case this is an issue. I freely give all of my patches and
changes to the R source code to whomever holds or keeps the copyright/
rights for this project. If it matters I have not signed any FSF
agreements but I have done so with Sun for OpenOffice.org and would
be willing to
Hi,
Since no one answered my questions about the "theoretically correct
values of named for that LANGSXP object, I made an assumption that
the problem was in the creation of the call2 object to begin with.
Yes, the bug was there.
Inline is the patch to fix bug 7924.
It seems that in subset.
Hi,
On Jun 5, 2006, at 8:02 AM, Peter Dalgaard wrote:
> Not quite: more like freshly-made-not-assigned,
> assigned-but-only-once, assigned-maybe-more-than-once.
So for my particular case ...
> call2 <- Quote(f(arg[[1]]))[c(1,2,2,2)]
> 0: 0x9e7d18 LANGSXP Object with length 1, named 1
>
Hi,
Okay I threw together a quick dump_object routine and found something
that I don't think is correct when call2 is created.
> call2 <- Quote(f(arg[[1]]))[c(1,2,2,2)]
> get("call2")
I use the do_get break to find the SEXP value I want
Breakpoint 1, do_get (call=0xc2d530, op=0x52bd30, args
Hi,
I am working my way through R-exts.pdf and have found R_PV but it
does not do a full dump of an object.
So I was wondering since I can call C code from with gdb (with call)
if there was a routine someplace that somebody has that would dump
all of the objects and fields and any important
Hi,
Again, sorry for the length of this post. Once I get my new office I
will get a website set up on my work machine and will simply post a
link to the log since I doubt many people are truly interested in
these logs.
To further analyze what is happening, I added my own routine in
main.
Hi,
Please forgive me but I am just learning how to debug in this
environment.
Info on Bug 7924, as far as I understand the comment, if someone knew
how calling identical(call1,call2) actually helped set NAMED it might
lead them to a solution.
Here is my attempt to track down when how this
Hi Thomas,
> A key fact here is that NAMED(object) is 0 if the object is not
> (part of) an R variable, 1 if it is (part of) exactly one R
> variable, 2 if it is part of more than one R variable.
> The point is that NAMED=0 or 1 objects can be safely modified, but
> NAMED=2 have to be copi
Hi All,
Well I finally have found the time to get svn working and I have
successfully built my own tuned atlas (multi-threaded version) libs
and have both the r-devel and r-patched trees building daily on my box.
The problem is I still do not have a good idea of the layout and
design of R,
Hi,
I have built my own tuned atlas libs enabling multiple processors and
pthread support. Then based on the note on page 23 of the "R
Installation and Administration" manual, I added --disable-R-
profiling to the configuration options
Unfortunately, the code in R/src/main/eval.c [updated f
Hi,
I just joined the list and wanted to introduce myself. I am a
professor of operations management at a Canadian University and an
empirical researcher (using mainly archival data from finance and
accounting primarily with econometric methods). I have finally made
the leap to remove al
35 matches
Mail list logo
