Displaying 2 results from an estimated 2 matches for "tukeyline0".
Did you mean:
tukeyline
2017 May 31
2
stats::line() does not produce correct Tukey line when n mod 6 is 2 or 3
...}
+ R_rsort(z,n);
+ yint = 0.5 * (z[il(n, 0.5)] + z[iu(n, 0.5)]);
}
for( i = 0 ; i < n ; i++ ) {
- w[i] = yint + slope*x[i];
- z[i] = y[i] - w[i];
+ w[i] = yint + slope*x[i];
+ z[i] = y[i] - w[i];
}
coef[0] = yint;
coef[1] = slope;
}
-void tukeyline0(double *x, double *y, double *z, double *w, int *n,
- double *coef)
+void tukeyline0(double *x, double *y, double *z, double *w, int *indx, int *n,
+ double *coef)
{
- line(x, y, z, w, *n, coef);
+ line(x, y, z, w, indx, *n, coef);
}
SEXP tukeyline(SEXP x, SEXP y, SEXP c...
2017 May 31
4
stats::line() does not produce correct Tukey line when n mod 6 is 2 or 3
Seriously, if a method gives a wrong result, it's wrong. line() does NOT
implement the algorithm of Tukey, even not after the patch. We're not
discussing Excel here, are we?
The method of Tukey is rather clear, and it is NOT using the default
quantile definition from the quantile function. Actually, it doesn't even
use quantiles to define the groups. It just says that the groups