関連情報
ホーム > 製品 > NAG数値計算ライブラリ > サンプルソースコード集 > 二変量Frank copulaから疑似乱数行列を生成 (C言語/C++)

二変量Frank copulaから疑似乱数行列を生成

C言語によるサンプルソースコード
使用関数名:nag_rand_bivariate_copula_frank (g05rfc)

Keyword: 二変量, Frank copula, 疑似乱数

概要

本サンプルは二変量Frank copulaから疑似乱数行列の生成を行うC言語によるサンプルプログラムです。 本サンプルは以下の式で表されるFrank copulaのC-12.0の分布をもつ13の2変量疑似乱数を生成し出力します。

frankのデータ 

※本サンプルはNAG Cライブラリに含まれる関数 nag_rand_bivariate_copula_frank() のExampleコードです。本サンプル及び関数の詳細情報は nag_rand_bivariate_copula_frank のマニュアルページをご参照ください。
ご相談やお問い合わせはこちらまで

出力結果

(本関数の詳細はnag_rand_bivariate_copula_frank のマニュアルページを参照)

この出力例をダウンロード
nag_rand_bivariate_copula_frank (g05rfc) Example Program Results

Uniform variates with copula joint distrbution
  0.636373    0.141124
  0.106480    0.896697
  0.746035    0.184315
  0.798270    0.125439
  0.104592    0.998172
  0.492451    0.690103
  0.384260    0.625047
  0.787070    0.165409
  0.498187    0.529766
  0.671703    0.290232
  0.050537    0.955405
  0.258030    0.818964
  0.623806    0.301387

  • 3〜16行目に生成された疑似乱数が出力されています。

ソースコード

(本関数の詳細はnag_rand_bivariate_copula_frank のマニュアルページを参照)

※本サンプルソースコードはNAG数値計算ライブラリ(Windows, Linux, MAC等に対応)の関数を呼び出します。
サンプルのコンパイル及び実行方法


このソースコードをダウンロード
/* nag_rand_bivariate_copula_frank (g05rfc) Example Program.
 *
 * CLL6I261D/CLL6I261DL Version.
 *
 * Copyright 2017 Numerical Algorithms Group.
 *
 * Mark 26.1, 2017.
 */
/* Pre-processor includes */
#include <stdio.h>
#include <math.h>
#include <nag.h>
#include <nag_stdlib.h>
#include <nagg05.h>

#define X(I, J) x[order == Nag_ColMajor?((J-1)*pdx + I-1):((I-1)*pdx + J-1)]

int main(void)
{
  /* Integer scalar and array declarations */
  Integer exit_status = 0;
  Integer i, lstate, pdx, sdx;
  Integer *state = 0;

  /* Double scalar and array declarations */
  double *x = 0;

  /* NAG structures */
  NagError fail;

  /* Use row major order */
  Nag_OrderType order = Nag_RowMajor;

  /* Set the number of variates */
  Integer n = 13;

  /* Choose the base generator */
  Nag_BaseRNG genid = Nag_Basic;
  Integer subid = 0;

  /* Set the seed */
  Integer seed[] = { 1762543 };
  Integer lseed = 1;

  /* Set the theta parameter value */
  double theta = -12.0e0;

  /* Initialize the error structure */
  INIT_FAIL(fail);

  printf("nag_rand_bivariate_copula_frank (g05rfc) "
         "Example Program Results\n\n");

  /* Get the length of the state array */
  lstate = -1;
  nag_rand_init_repeatable(genid, subid, seed, lseed, state, &lstate, &fail);
  if (fail.code != NE_NOERROR) {
    printf("Error from nag_rand_init_repeatable (g05kfc).\n%s\n",
           fail.message);
    exit_status = 1;
    goto END;
  }

  /* Set matrix size and principal dimension according to storage order */
  pdx = (order == Nag_ColMajor) ? n : 2;
  sdx = (order == Nag_ColMajor) ? 2 : n;

  /* Allocate arrays */
  if (!(x = NAG_ALLOC((pdx * sdx), double)) ||
      !(state = NAG_ALLOC(lstate, Integer)))
  {
    printf("Allocation failure\n");
    exit_status = -1;
    goto END;
  }

  /* Initialize the generator to a repeatable sequence */
  nag_rand_init_repeatable(genid, subid, seed, lseed, state, &lstate, &fail);
  if (fail.code != NE_NOERROR) {
    printf("Error from nag_rand_init_repeatable (g05kfc).\n%s\n",
           fail.message);
    exit_status = 1;
    goto END;
  }

  /* Generate variates */
  nag_rand_bivariate_copula_frank(order, state, theta, n, x, pdx, sdx, &fail);
  if (fail.code != NE_NOERROR) {
    printf("Error from nag_rand_bivariate_copula_frank (g05rfc).\n"
           "%s\n", fail.message);
    exit_status = 1;
    goto END;
  }

  /* Display the results */
  printf("Uniform variates with copula joint distrbution\n");
  for (i = 1; i <= n; i++) {
    printf(" %9.6f   %9.6f\n", X(i, 1), X(i, 2));
  }

END:
  NAG_FREE(x);
  NAG_FREE(state);

  return exit_status;
}


Results matter. Trust NAG.

Privacy Policy | Trademarks