Keyword: 分位数, 高速
概要
本サンプルは分位数を求めるC言語によるサンプルプログラムです。 本サンプルは以下に示される並べ替えのされていないベクトルについて高速のアルゴリズムを用いて特定の分位数を求め、出力します。
※本サンプルはnAG Cライブラリに含まれる関数 nag_double_quantiles() のExampleコードです。本サンプル及び関数の詳細情報は nag_double_quantiles のマニュアルページをご参照ください。
ご相談やお問い合わせはこちらまで
入力データ
（本関数の詳細はnag_double_quantiles のマニュアルページを参照）| このデータをダウンロード |
nag_double_quantiles (g01amc) Example Program Data 22 5 0.5 0.729 0.861 0.44 0.791 0.001 0.062 0.912 0.27 0.141 0.32 0.133 0.654 0.285 0.553 0.438 0.316 0.696 0.718 0.293 0.704 0.029 0.0 0.25 0.73 0.9 1.0
- １行目はタイトル行で読み飛ばされます。
- ２行目にベクトルの要素の数(n)を指定しています。
- ３行目に分位数の数(nq)を指定しています。
- ４〜５行目にベクトルの要素(rv)を指定しています。
- ６行目に５個の分位(q)を指定しています。
出力結果
（本関数の詳細はnag_double_quantiles のマニュアルページを参照）| この出力例をダウンロード |
nag_double_quantiles (g01amc) Example Program Results
Quantile Result
0.0000 0.0010
0.2500 0.2737
0.7300 0.6986
0.9000 0.7848
1.0000 0.9120
- ４から８行目に分位と分位数の値が出力されています。
ソースコード
（本関数の詳細はnag_double_quantiles のマニュアルページを参照）
※本サンプルソースコードはnAG数値計算ライブラリ（Windows, Linux, MAC等に対応）の関数を呼び出します。
サンプルのコンパイル及び実行方法
| このソースコードをダウンロード |
/* nag_double_quantiles (g01amc) Example Program.
*
* CLL6I261D/CLL6I261DL Version.
*
* Copyright 2017 Numerical Algorithms Group.
*
* Mark 26.1, 2017.
*/
#include <nag.h>
#include <stdio.h>
#include <nag_stdlib.h>
#include <nag_stddef.h>
#include <nagg01.h>
int main(void)
{
/* Scalars */
Integer exit_status = 0, i, n, nq;
/* Arrays */
double *vec = 0, *quants = 0, *quant_vec = 0;
/* Nag Types */
NagError fail;
INIT_FAIL(fail);
/* Skip heading in data file */
scanf("%*[^\n]");
printf("nag_double_quantiles (g01amc) Example Program Results\n");
scanf("%ld", &n);
scanf("%ld", &nq);
if (n >= 1 && nq >= 1) {
if (!(vec = nAG_ALLOC(n, double)) ||
!(quants = nAG_ALLOC(nq, double)) ||
!(quant_vec = nAG_ALLOC(nq, double)))
{
printf("Allocation failure\n");
exit_status = -1;
goto END;
}
}
else {
if (n < 1) {
printf("Invalid n.\n");
}
else {
printf("Invalid nq.\n");
}
exit_status = 1;
goto END;
}
for (i = 0; i < n; ++i)
scanf("%lf", &vec[i]);
for (i = 0; i < nq; ++i)
scanf("%lf", &quants[i]);
/* nag_double_quantiles (g01amc).
* Find quantiles of set of values of data type double
*/
nag_double_quantiles(n, vec, nq, quants, quant_vec, &fail);
if (fail.code != NE_NOERROR) {
printf("Error from nag_double_quantiles (g01amc).\n%s\n", fail.message);
exit_status = 1;
goto END;
}
printf(" Quantile Result\n\n");
for (i = 0; i < nq; ++i) {
printf(" %7.4f %7.4f\n", quants[i], quant_vec[i]);
}
END:
nAG_FREE(vec);
nAG_FREE(quants);
nAG_FREE(quant_vec);
return exit_status;
}