similar to: What algorithm is R using to calculate mean?

Hello, I am trying to write FORTRAN code to do the same as some R code I have. I get (small) differences when using the sum function in R. I know there are numerical routines to improve precision, but I have not been able to figure out what algorithm R is using. Does anyone know this? Or where can I find the code for the sum function? Regards, Rampal Etienne

The algorithm does make a differece. You can use Kahan's summation algorithm (https://en.wikipedia.org/wiki/Kahan_summation_algorithm) to reduce the error compared to the naive summation algorithm. E.g., in R code: naiveSum <- function(x) { s <- 0.0 for(xi in x) s <- s + xi s } kahanSum <- function (x) { s <- 0.0 c <- 0.0 # running compensation for lost

Dear Tomas, Where do I find these files? Do they contain the code for the sum function? What do you mean exactly with your point on long doubles? Where can I find documentation on this? Cheers, Rampal On Mon, Feb 18, 2019, 15:38 Tomas Kalibera <tomas.kalibera at gmail.com wrote: > See do_summary() in summary.c, rsum() for doubles. R uses long double > type as accumulator on systems

Specifically: https://svn.r-project.org/R/trunk/src/main/summary.c And if you don't want to deal with Subversion, you can look at the read-only github mirror: https://github.com/wch/r-source/blob/e5b21d0397c607883ff25cca379687b86933d730/src/main/summary.c#L115-L131 On Thu, Feb 21, 2019 at 11:57 AM David Winsemius <dwinsemius at comcast.net> wrote: > > > On 2/20/19 2:55 PM,

Hi all,? The survdiff() from survival package has an argument "rho" that implements Fleming-Harrington weighted long rank test.? But according to several sources including "survminer" package (https://cran.r-project.org/web/packages/survminer/vignettes/Specifiying_weights_in_log-rank_comparisons.html), Fleming-Harrington weighted log-rank test should have 2 parameters

Dear list, I spent about two hours searching on the message archive, with no avail. I have a list of people that have to pass an on-line test, but only a fraction of them do it. Moreover, as they input their names, the resulting string do not always match the names I have in my database. I would like to do two things: 1. Match any strings that are 90% the same Example: name1 <- "Harry

Hi, I am working on vxml to record voice. I have trouble with getting url of recorded voice. I want to save and I am using java to get record parameter from url and it returns a string which is audio/basic:len(123123):p0x5a6e6241, but I want to get file object or base64 string with parameter or to relate returning string with path in asterisk server, are there any way to do this? --

> On Feb 13, 2018, at 4:02 PM, array chip via R-help <r-help at r-project.org> wrote: > > Hi all, > > The survdiff() from survival package has an argument "rho" that implements Fleming-Harrington weighted long rank test. > > But according to several sources including "survminer" package

> On Feb 14, 2018, at 5:26 PM, David Winsemius <dwinsemius at comcast.net> wrote: > >> >> On Feb 13, 2018, at 4:02 PM, array chip via R-help <r-help at r-project.org> wrote: >> >> Hi all, >> >> The survdiff() from survival package has an argument "rho" that implements Fleming-Harrington weighted long rank test. >>

All, I'm trying to send an SMS directly from asterisk but it doesn't seem to be working. The SMS() function does create an outgoing file but doesn't deliver the SMS. Can anyone help me to understand how SMS() works. Thanks. extensions.conf example: same => n,SMS(hello,a,17654307001,"hello nick") - nick

https://bugs.freedesktop.org/show_bug.cgi?id=34371 Summary: [natty] Video corruption on kernel 2.6.38-1-generic (and on -3) and nVidia Corporation GT216 [GeForce GT 230M] (rev a2) Product: xorg Version: 7.6 Platform: x86-64 (AMD64) OS/Version: Linux (All) Status: NEW Keywords:

We have around 30+ users using one computer (Win2kTS) So I guess that means one smbd process is handling all 30 users' requests? You can watch the smbd process for this machine max out the CPU at times Is there anyway of getting multiple smbd processes per client/computer that connects Thanks, Paul -------------- next part -------------- A non-text attachment was scrubbed... Name:

I need a hosting provider just like Dreamhost.net's shared hosting service, but with a CentOS shell instead of Debian. Any pointers? Matt

Previously I had problems with data(foo) and Gass and I had problems with help(). His diagnostic was "POSIX regular expressions not available" and mine was "invalid regular expression". Martin Maechler has no problems with Solaris 2.5 and suggested a regular expression library was needed, e.g., regex-x.xx. Solaris 2.4 has regex.o in /usr/ccs/libgen.a and regex.h in

R-0.62.1 sparc-sun-solaris2.4: demo(graphics) stops at: > data("iris") Error in sub(pattern, replacement, x, ignore.case, extended) : invalid regular expression or these direct entries give: > data(iris) Error in sub(pattern, replacement, x, ignore.case, extended) : invalid regular expression > > help(iris) Error in gsub(pattern, replacement, x, ignore.case, extended) :

I'm open to flames if this has already been posted. http://www.sparkfun.com/commerce/product_info.php?products_id=8300 As of April 17th, SF listed a new Ogg Vorbis decoder chip that is seriously listing for $10 USD. I'm excited. I think I'll buy a couple and integrate them into an IR streaming project I've been tinkering with. -Chris Harrington

I'm moving my Ogg Vorbis collection off of my Linux server and onto my laptop. I plan on using iTunes to play my collection, but, that's a whole 'nother can of worms. I'm having trouble moving my complete collection over because dbPowerAmp (an application I used to love) made some dumb decisions about naming my files. For example, in track/album/artist names that contained a

Need your help! I'm very confised and tired, tried a lot of variants but still at the start, even lost all understanding how the permissions and masks works. Please write me a config for the following case: There is a directory: '/home/shared/door', owned by 'michael:office' shared via samba as '[door]' How can I configure the FS and samba to let 'michael'

Can someone point me to an FAQ, or tell me how to access shared PC drives from Solaris 8? I've heard of the smbclient thing - does this allow seamless access - ie can I run this, then say cd /sharedpc/c/somefolder? Thanks for any help Paul

Hi, I am currently doing a research/development project with Rails as part of my Computer Science degree (my final year project). I have been looking for those who have had previous experience with Rails in industry. Using Rails to create applications for a business, or using Rails at work to create applications. Attached is a txt document with a few questions, designed to get the information I