From ivo at iugrina.com  Fri Dec 14 23:03:37 2007
From: ivo at iugrina.com (Ivo Ugrina)
Date: Fri, 14 Dec 2007 23:03:37 +0100
Subject: [Rd] Improvement of SignRank functions
Message-ID: <4762FDB9.1060603@iugrina.com>

I took some time and liberty and tried to improve
existing implementation of SignRank functions
in R. (dsignrank, ...)

As I have seen they've been based on csignrank.
So I modified csignrank and, I believe,
improved calculation time and memory efficiency.

The idea is basically the same. I use the same recursion
as original author used with one slight modification.
I am generating Wilcoxon SignRank density from the
beginning (for n=2,3,...) with the help of recursion formula.

What is changed?
There is no need for SINGRANK_MAX in src/nmath/nmath.h anymore.
Only functions:
static void w_free()
void signrank_free()
static void w_init_maybe(int n)
static double csignrank(int k, int n)
in src/nmath/signrank.c have been altered.
There was no change to dsignrank, psignrank, ...

I've tried to make as little changes as possible.
So, to compile with new functions only src/nmath/signrank.c
needs to be changed with the one given in attachment.

I hope it really is an improvement.

With respect,
-- 
Ivo Ugrina
ICQ: 47508335 | www.iugrina.com
-------------------------------
baza matematickih pojmova
http://baza.iugrina.com
---------------------------
anime, manga, Japan fanzin
http://yoshi.iugrina.com
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: signrank.c
Url: https://stat.ethz.ch/pipermail/r-devel/attachments/20071214/ab3480ad/attachment.c 


[Rd] Improvement of SignRank functions
Ivo Ugrina