「統計学関連なんでもあり」の過去ログ---018

★ R ** 計算誤差 ★

　92　R ** 計算誤差　　青木繁伸　　2002/02/04 (月) 21:26
　　118　Re: R ** 計算誤差　　青木繁伸　　2002/02/07 (木) 11:54
　　95　Re: R ** 計算誤差　　伊達　　2002/02/05 (火) 12:17

92.　R ** 計算誤差　　青木繁伸　　2002/02/04 (月) 21:26

Macintosh 版の
R : Copyright 2001, The R Development Core Team
Version 1.4.0 (2001-12-19)
において，
> sum(scale(1:10)^3)
[1] -4.440892e-16
となってしまいますが，ほかのプラットホームでも同じでしょうか?
理論的には0で，よほど変なプログラムを書かない限りこんな値が出るはずないのですが。ま，-4.440892e-16 ほとんどゼロだからどうでもよかろうといえなくもないですが，気持ち悪いです。

　　　　　[このページのトップへ]

118.　Re: R ** 計算誤差　　青木繁伸　　2002/02/07 (木) 11:54

大きさの違う数値を足すときには注意しないといけないというのを忘れていました。
http://aoki2.si.gunma-u.ac.jp/Hanasi/Algo/letsc/1-body.html
の，3番目の項に書いてあるのにね(^_^)

my.sum <- function(x) { sum(x[order(abs(x))]) } というのを使うとこの問題は解消されます。　　（注:これでもだめです。実数値が近似値であることがそもそもの原因） #中林さんありがとうございました。

　　　　　[このページのトップへ]

95.　Re: R ** 計算誤差　　伊達　　2002/02/05 (火) 12:17

Windows版です。OSはWindows2000です。

R 1.4.0 - A Language and Environment
Copyright 2001,
The R Development Core Team

> sum(scale(1:10)^3)
[1] -4.440892e-16

となりました。

　　　　　[このページのトップへ]