DataMuseum.dk

top - download

    Length: 866 (0x362)
    Types: TextFile
    Names: »lsd.c«

└─⟦87ddcff64⟧ Bits:30001253 CPHDIST85 Tape, 1985 Autumn Conference Copenhagen
    └─ ⟦this⟧ »cph85dist/stat/src/lsd.c« 

TextFile

#include "unixstat.h"
PGM(lsd,Least Significant Differences and Confidence Intervals,0.0,1985)
/* Copyright (c) 1982 Gary Perlman (see Copyright file) */

double	critf ();
double
scheffe (p, df1, df2, mse, n) double p, mse;
	{
	return (sqrt (df1 * critf (p, df1, df2) * mse * 2 / n));
	}
double
lsd (p, df1, df2, mse, n) double p, mse;
	{
	return (sqrt (critf (p, df1, df2) * mse * 2 / n));
	}
main (argc, argv) char **argv;
	{
	double	prob;
	int 	df1;
	int 	df2;
	double	mse;
	int 	n;
	ARGV0;
	if (argc != 6)
		USAGE (prob df1 df2 mse n)
	prob = atof (argv[1]);
	df1 = atoi (argv[2]);
	df2 = atoi (argv[3]);
	mse = atof (argv[4]);
	n = atoi (argv[5]);
	puts ("Least Significant Difference and Scheffe's Method");
	printf ("LSD(%.2f) = %.3g\n", prob, lsd (prob, 1, df2, mse, n));
	printf ("S(%.2f)   = %.3g\n", prob, scheffe (prob, df1, df2, mse, n));
	exit (0);
	}