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

スチューデントt-copulaから疑似乱数行列を生成

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

Keyword: スチューデントt-copula, 疑似乱数

概要

本サンプルはスチューデントt-copulaから疑似乱数行列の生成を行うC言語によるサンプルプログラムです。 本サンプルは以下に示される共分散行列をもつスチューデントt-copulaから10個の疑似乱数を生成し出力します。

スチューデントt-copulaのデータ 

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

出力結果

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

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

     0.6445    0.0527    0.4082    0.8876 
     0.0701    0.1988    0.8471    0.3521 
     0.7988    0.6664    0.2194    0.5541 
     0.8202    0.0492    0.7059    0.9341 
     0.1786    0.5594    0.7810    0.2836 
     0.4920    0.2677    0.3427    0.5169 
     0.4139    0.2978    0.8762    0.7145 
     0.7437    0.9714    0.8931    0.2487 
     0.4971    0.9687    0.8142    0.1965 
     0.6464    0.5304    0.5817    0.4565 

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

ソースコード

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

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


このソースコードをダウンロード
/* nag_rand_copula_students_t (g05rcc) 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*pdx + I):(I*pdx + J)]

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

  /* NAG structures */
  NagError fail;
  Nag_ModeRNG mode;

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

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

  /* Set the number of variables and variates */
  Integer m = 4;
  Integer n = 10;

  /* Input the covariance matrix */
  double c[] = { 1.69e0, 0.39e0, -1.86e0, 0.07e0,
    0.39e0, 98.01e0, -7.07e0, -0.71e0,
    -1.86e0, -7.07e0, 11.56e0, 0.03e0,
    0.07e0, -0.71e0, 0.03e0, 0.01e0
  };
  Integer pdc = 4;

  /* Set the degrees of freedom */
  Integer df = 10;

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

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

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

  printf("nag_rand_copula_students_t (g05rcc) 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;
  }

  pdx = (order == Nag_ColMajor) ? n : m;
  x_size = (order == Nag_ColMajor) ? pdx * m : pdx * n;

  /* Calculate the size of the reference vector */
  lr = m * m + m + 2;

  /* Allocate arrays */
  if (!(r = NAG_ALLOC(lr, double)) ||
      !(x = NAG_ALLOC(x_size, 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;
  }

  /* Set up reference vector and generate variates */
  mode = Nag_InitializeAndGenerate;
  nag_rand_copula_students_t(order, mode, n, df, m, c, pdc, r, lr, state,
                             x, pdx, &fail);
  if (fail.code != NE_NOERROR) {
    printf("Error from nag_rand_copula_students_t (g05rcc).\n%s\n",
           fail.message);
    exit_status = 1;
    goto END;
  }

  /* Display the variates */
  for (i = 0; i < n; i++) {
    printf("  ");
    for (j = 0; j < m; j++)
      printf("%9.4f%s", X(i, j), (j + 1) % 10 ? " " : "\n");
    if (m % 10)
      printf("\n");
  }

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

  return exit_status;
}


Results matter. Trust NAG.

Privacy Policy | Trademarks