Perlゼミ

簡単な計算を行おう | Rstats入門(統計解析・科学技術計算)

簡単な計算を行ってみましょう。Rstatsでは、以下の演算子を使って、計算が行えます。

四則演算

足し算、引き算、掛け算、割り算の記号は、R言語と同じです。

累乗

累乗はPerlの演算規則に従って「**」を使って行います。R言語の「^」とは異なるので注意してください。

剰余

剰余はPerlの演算規則に従って「%」を使って行います。R言語の「%%」とは異なるので注意してください。

整数商

Rstatsでは整数商を計算するための演算子は用意されていません。R言語の「%/%」は利用できないので注意してください。整数商を求める場合は、tranc関数を使って行います。割り算を行った後に、整数部を取り出します。

r->tranc($x1 / $x2)

計算の例を示します。Rstatsでは、ベクトルを使った計算が基礎になるので、c_関数を使って作成します。c_関数はベクトルを作成するための関数です。ひとつの数値であっても、ひとつの要素を持つベクトルとして計算します。

use strict;
use warnings;
use Rstats;

my $x = (c_(1) + c_(2) - c_(3) * c_(4)) / c_(5) ** c_(6);
print $x;

出力結果です。計算結果が出力されます。

[1] -0.000576

c_()をいつでも使うのは多少面倒に感じるかもしれませんが、Perlの上にR言語を実装しているために必要です。Rstatsで計算をするときは、いつでもベクトルを使って計算をするということを覚えてください。

今後のサンプルでは簡便のための、上記の「use Rstats」の下の行から書いていきますので、必要に応じて補ってください。

ルートの計算、対数の計算を行う場合には、Rstatsにはたくさんの関数が用意されています。たとえば、ルートの計算を行うためには、「sqrt関数」を使用します。

my $x = r->sqrt(c_(2));
print $x;

出力結果。

[1] 1.4142135623731

以下にRstatsで用意されている数学関数を紹介します。

Rstats入門(統計解析・科学技術計算)