Paul BIGNIER [Tue, 2 Apr 2013 08:07:33 +0000 (10:07 +0200)]
For future changes in them (commit to follow)

Change-Id: I2d90a59995ef38577682ec53a9e379657032e864

@@ -1,32 +1,32 @@

// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) 2000 - INRIA - Carlos Klimann
-//
+//
// This file must be used under the terms of the CeCILL.
// This source file is licensed as described in the file COPYING, which
// you should have received as part of this distribution.  The terms
// are also available at
// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
-//
+//

function [f,p]=ftest(samples)
-//
-//Given a  number nc of samples each  of them composed of nr
-//observations this fonction computes in f the Fischer ratio
-//(it is the  ratio between  nr times  the variance   of the
-//means of  samples  and the  mean of the  variances of each
-//sample).  Additionnally it gives (in p) the p-value of the
-//computed  Fischer ratio.  samples   is a matrix  with type
-//nrxnc.
-//
-//References:  Wonacott, T.H. & Wonacott, R.J.; Introductory
-//Statistics, J.Wiley & Sons, 1990.
-//
-  [nr nc]=size(samples)
-  xbar=mean(samples,'r')
-  xbarbar=mean(xbar)
-  sx2=sum((xbar-xbarbar).^2)/(nc-1)
-  sp2=(sum(sum((samples-ones(nr,1)*xbar).^2),'r'))/(nc*(nr-1))
-  f=(nr*sx2)/sp2
-  p=1-cdff("PQ",f,(nc-1),(nc*(nr-1)))
+    //
+    //Given a  number nc of samples each  of them composed of nr
+    //observations this fonction computes in f the Fischer ratio
+    //(it is the  ratio between  nr times  the variance   of the
+    //means of  samples  and the  mean of the  variances of each
+    //sample).  Additionnally it gives (in p) the p-value of the
+    //computed  Fischer ratio.  samples   is a matrix  with type
+    //nrxnc.
+    //
+    //References:  Wonacott, T.H. & Wonacott, R.J.; Introductory
+    //Statistics, J.Wiley & Sons, 1990.
+    //
+    [nr nc]=size(samples)
+    xbar=mean(samples,"r")
+    xbarbar=mean(xbar)
+    sx2=sum((xbar-xbarbar).^2)/(nc-1)
+    sp2=(sum(sum((samples-ones(nr,1)*xbar).^2),"r"))/(nc*(nr-1))
+    f=(nr*sx2)/sp2
+    p=1-cdff("PQ",f,(nc-1),(nc*(nr-1)))
endfunction
index 1244513..d6ef333 100644 (file)
@@ -1,55 +1,55 @@

// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) 2000 - INRIA - Carlos Klimann
-//
+//
// This file must be used under the terms of the CeCILL.
// This source file is licensed as described in the file COPYING, which
// you should have received as part of this distribution.  The terms
// are also available at
// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
-//
+//

function [f,p]=ftuneq(varargin)
-//
-//This function computes the F ratio for samples of unequal
-//size.
-//
-//"The most  efficient design is  to make all  samples the
-//same  size n.   However when  this is  nor  feasible, it
-//still  is possible  to modify  the  ANOVA calculations."
-//Note  that  the  definition  of  xbarbar  is  no  longer
-//mean(xbar), but  rather a weighted  average with weights
-//ni.  Additionnally  it gives (in  p) the p-value  of the
-//computed Fischer ratio.
-//
-//Given a number  a of samples each of  them composed of n_i
-//(i from 1  to a) observations this fonction  computes in f
-//the Fischer  ratio (it is  the ratio between nr  times the
-//variance  of the  means of  samples  and the  mean of  the
-//variances of each sample).
-//
-//References:  Wonacott, T.H. & Wonacott, R.J.; Introductory
-//Statistics, J.Wiley & Sons, 1990.
-//
-//
-  [lhs,rhs]=argn(0)
-  data=[]
-  total=0
-  sse=0
-  for i=1:rhs
-    len(i)=size(varargin(i),'*'),
-    mat=matrix(varargin(i),len(i),1),
-    partial=sum(mat),
-    xbar(i)=partial/len(i),
-    total=total+partial,
-    data=[data; mat]
-    sse=sse+sum((mat-xbar(i)).^2)
-  end
-  xbarbar=total/sum(len)
-  ssa=len'*((xbar-xbarbar).^2)
-  msa=ssa/(rhs-1)
-  deglib=sum(len-1)
-  mse=sse/deglib
-  f=msa/mse
-  p=1-cdff("PQ",f,(rhs-1),deglib)
+    //
+    //This function computes the F ratio for samples of unequal
+    //size.
+    //
+    //"The most  efficient design is  to make all  samples the
+    //same  size n.   However when  this is  nor  feasible, it
+    //still  is possible  to modify  the  ANOVA calculations."
+    //Note  that  the  definition  of  xbarbar  is  no  longer
+    //mean(xbar), but  rather a weighted  average with weights
+    //ni.  Additionnally  it gives (in  p) the p-value  of the
+    //computed Fischer ratio.
+    //
+    //Given a number  a of samples each of  them composed of n_i
+    //(i from 1  to a) observations this fonction  computes in f
+    //the Fischer  ratio (it is  the ratio between nr  times the
+    //variance  of the  means of  samples  and the  mean of  the
+    //variances of each sample).
+    //
+    //References:  Wonacott, T.H. & Wonacott, R.J.; Introductory
+    //Statistics, J.Wiley & Sons, 1990.
+    //
+    //
+    [lhs,rhs]=argn(0)
+    data=[]
+    total=0
+    sse=0
+    for i=1:rhs
+        len(i)=size(varargin(i),"*"),
+        mat=matrix(varargin(i),len(i),1),
+        partial=sum(mat),
+        xbar(i)=partial/len(i),
+        total=total+partial,
+        data=[data; mat]
+        sse=sse+sum((mat-xbar(i)).^2)
+    end
+    xbarbar=total/sum(len)
+    ssa=len'*((xbar-xbarbar).^2)
+    msa=ssa/(rhs-1)
+    deglib=sum(len-1)
+    mse=sse/deglib
+    f=msa/mse
+    p=1-cdff("PQ",f,(rhs-1),deglib)
endfunction