Statistics function incorrectly computing median

		Miguel Garcia-Blanco – Statistics function incorrectly computing median Miguel Garcia-Blanco Statistics function incorrectly computing median - Octave 3.0.0 (i686-pc-msdosmsvc) The statistics function seems to be incorrectly computing the median: > x = 0:1; >... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Jan 6, 2008, at 1:30 PM, Miguel Garcia-Blanco wrote: > Octave 3.0.0 (i686-pc-msdosmsvc) > > The statistics... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Jan 6, 2008, at 2:14 PM, Ben Abbott wrote: > > On Jan 6, 2008, at 1:30 PM, Miguel Garcia-Blanco wrote: > >>... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		ororo ororo ----- Original Message ----- From: "Ben Abbott" <bpabbott@...> To: "Vercelli" <ororo@...> Sent: Sunday,... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Jan 6, 2008, at 10:40 PM, Vercelli wrote: > > ----- Original Message ----- From: "Ben Abbott"... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Miguel Garcia-Blanco Miguel Garcia-Blanco >>> What would the 1st quartile, median and 3rd quartile be for a >>> population of [0, 1, 2, 3, 4, 5]? >>> >>> My... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Jan 7, 2008, at 12:03 AM, Miguel Garcia-Blanco wrote: >>>> What would the 1st quartile, median and 3rd quartile... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Miguel Garcia-Blanco Miguel Garcia-Blanco > Please see the attached file and compare my results to R .... and I'd > like to see the results if they differ from... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Jan 7, 2008, at 7:17 AM, Miguel Garcia-Blanco wrote: >> Please see the attached file and compare my results to R... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Jan 6, 2008, at 8:36 PM, Ben Abbott wrote: > > On Jan 7, 2008, at 7:17 AM, Miguel Garcia-Blanco wrote: > >>>... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Miguel Garcia-Blanco Miguel Garcia-Blanco > I took another look and found a method that is more consistent with > the original algorithm. The original calculates... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		John W. Eaton John W. Eaton On 13-Jan-2008, Ben Abbott wrote: | I took another look and found a method that is more consistent with | the... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Jan 14, 2008, at 7:17 PM, John W. Eaton wrote: > On 13-Jan-2008, Ben Abbott wrote: > > | I took another look... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Jan 14, 2008, at 7:17 PM, John W. Eaton wrote: > On 13-Jan-2008, Ben Abbott wrote: > > | I took another look... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Miguel Garcia-Blanco Miguel Garcia-Blanco > I looked at the routines in R and Maxima and have concluded that they > assume the probability of each sample is a... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Jan 19, 2008, at 5:13 AM, Miguel Garcia-Blanco wrote: > > I looked at the routines in R and Maxima and have... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Miguel Garcia-Blanco Miguel Garcia-Blanco > It is clear, that a change is merited. I agree it makes sense to consider > "discrete", and "continuous"... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Jan 20, 2008, at 7:25 AM, Miguel Garcia-Blanco wrote: >> It is clear, that a change is merited. I agree it makes... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott Michael Godfrey ran some calculations using Matlab's quantile function. We can now compare Matlab's quantile (), R's... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Miguel Garcia-Blanco Miguel Garcia-Blanco > Michael Godfrey ran some calculations using Matlab's quantile function. > > We can now compare Matlab's quantile (),... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Muthiah Annamalai-3 Muthiah Annamalai-3 > >>> Michael Godfrey ran some calculations using Matlab's quantile >>> function. >>> >>> We can now compare... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Jan 21, 2008, at 2:19 AM, Miguel Garcia-Blanco wrote: >> Michael Godfrey ran some calculations using Matlab's... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Miguel Garcia-Blanco Miguel Garcia-Blanco > Are any of the licenses for any of those GPL compatible? > Acklam's code seems to be free: "This is a collection... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Jan 21, 2008, at 8:44 AM, Miguel Garcia-Blanco wrote: >> Are any of the licenses for any of those GPL... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Miguel Garcia-Blanco Miguel Garcia-Blanco > In the octave scripts, the routines take samples and sample probabilities > as inputs. > > Thus x = [1, 1, 2, 3] &... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Monday, January 21, 2008, at 09:36AM, "Miguel Garcia-Blanco" <miguel.01@...> wrote: >> In the octave... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Miguel Garcia-Blanco Miguel Garcia-Blanco > I was referring to Octave's statistics functions, specifically > discrete_inv.m > > Octave's present algorithms... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Jan 21, 2008, at 5:53 PM, Miguel Garcia-Blanco wrote: >> I was referring to Octave's statistics functions,... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Miguel Garcia-Blanco Miguel Garcia-Blanco > ... So what I plan to do is to leave discrete_???.m as it is, and have > empirical_???.m work as R's method #7. I'll... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Jan 22, 2008, at 5:50 AM, Miguel Garcia-Blanco wrote: >> ... So what I plan to do is to leave discrete_???.m as... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Miguel Garcia-Blanco Miguel Garcia-Blanco > Due to the need for consistency between the group of functions, > empirical_inv/pdf/cdf/rnd, I'm leaning your... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Jan 22, 2008, at 10:27 AM, Miguel Garcia-Blanco wrote: >> >> How is the following list for a short plan? >> >>... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Miguel Garcia-Blanco Miguel Garcia-Blanco > Ok. I decided to place Matlab compatibility as a high priority (when > applicable). I do not have access to Matlab's... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Jan 26, 2008, at 6:39 AM, Miguel Garcia-Blanco wrote: > My tests (using 40 samples with 6 observations... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Jan 26, 2008, at 6:39 AM, Miguel Garcia-Blanco wrote: > My tests (using 40 samples with 6 observations... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Miguel Garcia-Blanco Miguel Garcia-Blanco > I struggled quite a bit with #3. To be honest, I had little hope it > was working properly. The part below was added... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Jan 27, 2008, at 7:57 AM, Miguel Garcia-Blanco wrote: >> I struggled quite a bit with #3. To be honest, I had... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Miguel Garcia-Blanco Miguel Garcia-Blanco > Thanks Miguel, you script was very helpful. > > I've modified the algorithm to be consistent with the other cases,... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Monday, January 28, 2008, at 07:57AM, "Miguel Garcia-Blanco" <miguel.01@...> wrote: >> Thanks Miguel, you... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott – PATCH: Statistics function incorrectly computing median Ben Abbott PATCH: Statistics function incorrectly computing median - Here's the ChangeLog, patch, and new files (I've zipped them with their directory structure in place). 2008-01-23... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		John W. Eaton John W. Eaton On 28-Jan-2008, Ben Abbott wrote: | Here's the ChangeLog, patch, The patch failed to apply cleanly, probably... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Jan 29, 2008, at 2:14 AM, John W. Eaton wrote: > On 28-Jan-2008, Ben Abbott wrote: > > | Here's the ChangeLog,... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott – Re: Statistics function incorrectly computing median Ben Abbott Re: Statistics function incorrectly computing median - On Jan 29, 2008, at 2:14 AM, John W. Eaton wrote: > > Shouldn't prctile and quantile work on N-d arrays? The Matlab... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Jan 29, 2008, at 2:14 AM, John W. Eaton wrote: > > Shouldn't prctile and quantile work on N-d arrays? The Matlab... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott – Re: [CHANGESET]: Statistics function incorrectly computing median Ben Abbott Re: [CHANGESET]: Statistics function incorrectly computing median - changeset is attached. # HG changeset patch # User Ben Abbott <bpabbott@...> # Date 1204244909 18000 # Node... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		John W. Eaton John W. Eaton On 28-Feb-2008, Ben Abbott wrote: | changeset is attached. | +2008-02-28 Ben Abbott <bpabbott@...> | + |... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Mar 5, 2008, at 4:50 PM, John W. Eaton wrote: > On 28-Feb-2008, Ben Abbott wrote: > > | changeset is... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Mar 5, 2008, at 4:50 PM, John W. Eaton wrote: > On 28-Feb-2008, Ben Abbott wrote: > > | changeset is... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Jaroslav Hajek-2 Jaroslav Hajek-2 On Thu, Mar 6, 2008 at 3:44 AM, Ben Abbott <bpabbott@...> wrote: > > On Mar 5, 2008, at 4:50 PM, John W. Eaton... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Mar 6, 2008, at 2:46 AM, Jaroslav Hajek wrote: > On Thu, Mar 6, 2008 at 3:44 AM, Ben Abbott <bpabbott@...>... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Jaroslav Hajek-2 Jaroslav Hajek-2 On Thu, Mar 6, 2008 at 1:28 PM, Ben Abbott <bpabbott@...> wrote: > > > On Mar 6, 2008, at 2:46 AM, Jaroslav... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Mar 6, 2008, at 2:46 AM, Jaroslav Hajek wrote: > On Thu, Mar 6, 2008 at 3:44 AM, Ben Abbott <bpabbott@...>... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		David Bateman – RE: [CHANGESET]: Statistics function incorrectly computing median David Bateman RE: [CHANGESET]: Statistics function incorrectly computing median - Why can't for q = 1:n inv(k(q)) = interp1 (cdf, x, p(k(q)), "linear", "extrap"); ... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Jaroslav Hajek-2 Jaroslav Hajek-2 On Thu, Mar 6, 2008 at 1:28 PM, Ben Abbott <bpabbott@...> wrote: > > > On Mar 6, 2008, at 2:46 AM, Jaroslav... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Thursday, March 06, 2008, at 09:16AM, "Bateman David-ADB014" <David.Bateman@...> wrote: >Why can't... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Thursday, March 06, 2008, at 09:25AM, "Jaroslav Hajek" <highegg@...> wrote: >On Thu, Mar 6, 2008 at 1:28 PM,... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Jaroslav Hajek-2 Jaroslav Hajek-2 On Thu, Mar 6, 2008 at 6:02 PM, Ben Abbott <bpabbott@...> wrote: > > On Thursday, March 06, 2008, at 09:25AM,... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Jaroslav Hajek-2 Jaroslav Hajek-2 Hi Ben, in case you're interested, I've created a loop-free version of __quantile__ (attached as __quantile_vec__).... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Mar 9, 2008, at 10:43 AM, Jaroslav Hajek wrote: > Hi Ben, > > in case you're interested, I've created a... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Jaroslav Hajek-2 Jaroslav Hajek-2 On Sun, Mar 9, 2008 at 5:33 PM, Ben Abbott <bpabbott@...> wrote: > > On Mar 9, 2008, at 10:43 AM, Jaroslav Hajek... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Mar 9, 2008, at 1:15 PM, Jaroslav Hajek wrote: > On Sun, Mar 9, 2008 at 5:33 PM, Ben Abbott <bpabbott@...>... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Jaroslav Hajek-2 Jaroslav Hajek-2 On Sun, Mar 9, 2008 at 6:23 PM, Ben Abbott <bpabbott@...> wrote: > > On Mar 9, 2008, at 1:15 PM, Jaroslav Hajek... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Mar 9, 2008, at 1:48 PM, Jaroslav Hajek wrote: > On Sun, Mar 9, 2008 at 6:23 PM, Ben Abbott <bpabbott@...>... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott A vectorized changeset is attached. Ben # HG changeset patch # User Ben Abbott <bpabbott@...> # Date... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott I thought I'd remove roundb() from __quantile__.m and found an error in patched statistics.m. I'll send a new... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		Ben Abbott Ben Abbott On Mar 25, 2008, at 8:45 PM, Ben Abbott wrote: > I thought I'd remove roundb() from __quantile__.m and found an error... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink

		John W. Eaton John W. Eaton On 25-Mar-2008, Ben Abbott wrote: | | On Mar 25, 2008, at 8:45 PM, Ben Abbott wrote: | > I thought I'd remove... Loading...   Reply More Rate this Message: Reply to author Delete Print Show in classic view Show in thread view Permalink