• ベストアンサー

統計ソフトRで4倍精度の計算はどうするのでしょうか?

初心者ですが、ご教示お願いします。 Fortran などのような4倍精度の計算はRではどうするのでしょうか? 宜しくお願い致します。

  • iykm
  • お礼率67% (98/145)

質問者が選んだベストアンサー

  • ベストアンサー
  • backs
  • ベストアンサー率50% (410/818)
回答No.1

Rのことについてでしたら、群馬大学の青木先生の掲示板(http://hpmboard3.nifty.com/cgi-bin/thread.cgi?user_id=NCD01156&msg_per_page=10&def=10)でご質問されるとよろしいかと思います。

iykm
質問者

お礼

ご教示ありがとうございます。 早速、質問いたしました。

関連するQ&A

  • 単精度 倍精度

    単精度 倍精度という言葉をたまに聴きます。 camに使うときや、「体積計算などするときには倍精度でないと計算が出来ない」とか「単精度から倍精度に変換したから数値がおかしい」などと聞くことがあるのですが、そもそも「単精度 倍精度」どう意味なのか、具体的にどういう数値なのかご存知の方がいらっしゃったら教えてください。

  • 単精度計算と倍精度計算

    F(x)=(1-x・x)1/2乗において、xの値を0.9、0.99、…0.99999と変化させた場合の単精度計算と倍精度計算の結果を求めるプログラムを教えてください。

  • FFTの単精度と倍精度の計算時間

    FFTのプログラミングで計算時間を短縮するため、三角関数の計算にチェビシェフの三項漸化式を用いたところ、計算時間はかなり短縮されました。ただ精度が心配なので、組み込み関数の三角関数の倍精度結果と、漸化式の単精度および倍精度の結果を比較したところ、2000点を越えるあたりから単精度ではかなりの誤差が生じましたが、倍精度では問題ありませんでした。そこで質問ですが、単精度と倍精度の違いとして、「精度とメモリが倍」は分かるのですが、計算時間はどうなるのでしょうか。 参考:チェビシェフの三項漸化式 cos{(j+1)*h} = 2*cos(h) * cos(j*h) - cos{(j-1)*h} sin {(j+1)*h} = 2*cos(h) * sin(j*h) - sin{(j-1)*h}

  • 計算精度(限りなくゼロに近く)

    現在、fortranで大きな数値計算をしています。 そこで、計算回数が多く、その計算過程において大きい値と小さい値とを何回も計算していると誤差が次第に大きくなっていきます。 A=BとなるようにA(又はB)を求めるときの判別としてabs(A-B)<10**(x) としています。 このxに-300ぐらいの値をいれ判別したいのですがそれは可能ですか? もちろん計算には倍精度を使っていますがabs(A-B)が10**(-17)付近になると値が一定になってしまいます。これがfortranの限界なのでしょうか?

  • 単精度と倍精度

    「10進表示で0.3をIEEE745規格の 単精度(32ビット)と倍精度(64ビット)で表示しなさい。ただし、表示数値<0.3となるように。」 という問題が課題が与えられたのですが、わかりますでしょうか? 自分で調べてやったら単精度が 010000000 0100100110011001001100 となったのですが、出したときに「違う!」と言われました。 どうかよろしくお願いします。

  • Fortranの倍精度実数について

    こんにちは。 現在、Fortran 90でプログラムを作成しています。 その中で、整数の倍精度実数への型変換についての疑問がわきましたので、質問させていただきます。 以下の2つのプログラムで、計算がより速く、より精度よくできるのはAとBどちらなのでしょうか。 実際は下記のプログラムが、ループの中に入っているので少しでも計算時間を短くしたいのです。 よろしくお願いいたします。 program A(毎回dxを足す) ----------------------------- real(8) :: pi, dx real(8), dimension(1000) :: x integer i, j pi=atan(1d0)*4d0 dx=pi/5000d0 x=0d0 j=1 do i=1, 1000 x(i)=x(j)+dx j=i write(*,*)x(i) end do ----------------------------- program B(毎回dx*ループ回数を計算する) ----------------------------- real(8) :: pi, dx real(8), dimension(1000) :: x integer i, j pi=atan(1d0)*4d0 dx=pi/5000d0 x=0d0 do i=1, 1000 x(i)=dx*dble(i) write(*,*)x(i) end do -----------------------------

  • 倍の計算

    先程はありがとうございます。 30を37.5にしたら何倍かの計算ですが、 何倍か計算する場合、なんで37.5割る30と計算するのか教えて下さい。 私は算数が苦手で大きな数割る小さな数と解釈していますが、間違った解釈だとは分かっています。 こういうのは何年生で習うのですか? 勉強をやり直したいのですが。

  • Visual Studioのバージョンによる計算精度

    Visual StudioのC++で開発を行っております。 現在はちょと古いですがVisual Studio .NET 2003を使用しており、今度Visual Stuido 2005にバージョンあげて再ビルドしようと思ってます。今後はVisual Studio 2008も考えております。 ただ計算精度として、各バージョンで違いがでてしまうのではないかと少し不安です。基本的には倍精度での一致すれば問題はないと思ってます。

  • Octaveで単精度で計算するには

    数値解析ソフトOctaveで単精度で数値計算を行うには、どうしたらいいでしょうか。

  • 統計ソフトRについての質問です。

    統計ソフトRについての質問です。 四分位数の計算方法(コマンドではない)を知りたいのですが、どこを参照すればよいのでしょうか? Rでは四分位数の計算方法が複数あると聞いたのですが、どのように計算しているか、ヘルプのどこら辺に載っていますか?