R言語：コレスポンデンス分析時の欠損値の扱いについて

Question

発生している問題・エラーメッセージ

Rでコレスポンデンス分析を実施したいのですが、欠損値の処理がうまくいっていないようでエラーとなってしまいます。欠損値の処理の仕方か、もしくは元のファイルの修正が必要なのか、はたまたcorresp関数は欠損値を処理できないのでしょうか。。。初心者なものでどうぞ宜しくお願い致します。

該当のソースコード

x <- read.table("sales.txt", header=TRUE, sep=",",na.strings="*")
x

SIZE PRICE DISTANSE SALES
1 M 1500 2.3 1
2 L 1800 1.2 3
3 S 1300 2.4 2
4 M 1500 1.1 2
5 S 1300 0.8 2
6 L 1800 1.8 4
7 M 1500 1.3 2
8 M 1500 0.1 1
9 M 1500 0.2 1
10 S 1300 0.7 1
11 L 1800 NA 2
12 M 1500 2.3 4
13 M 1500 1.5 3
14 M 1500 1.3 3
15 S 1300 NA 2
16 L 1800 0.3 1
17 S 1300 0.4 2
18 L 1800 1.3 1

x.co <- corresp(x, nf=2)

x%%1 でエラー: 二項演算子の引数が数値ではありません

試したこと

na.rm = TRUEを加えてみましたが、同様のエラーでした。

x.co <- corresp(x, nf=2,na.rm = TRUE)

x%%1 でエラー: 二項演算子の引数が数値ではありません

Answer 1

試しているのはMASSライブラリのcorresp()でしょうか。だとすると、欠損値の有無ではなく、データの構造が基本的に違うのではないでしょうか。

参考:
http://d.hatena.ne.jp/bob3/20120219/p1

Comments

[fuusuke]

>KojiDoi様

[fuusuke]

回答いただきありがとうございました。ご指摘の通りMASSライブラリのcorresp()です。データ構造が違うということで、別のテストデータで試したところ、異なるエラーメッセージとなりました。

[fuusuke]

> x <- read.table("trip.txt", header=TRUE, sep=",", row.names=1,na.strings="*")
> x
温泉 テーマパーク リゾート
家族で 45 35 11
いっぱい遊ぶ 9 42 35
ゆっくり 28 23 25
豪華に 9 20 40
近場で 25 35 3
さわぐ 13 41 10
脱日常 20 25 10
リフレッシュ 33 12 12
のびのび 20 18 NA
> x.co <- corresp(x, nf=2)
if (any(x < 0 | x%%1 > 10 * sqrt(.Machine$double.eps))) warning("negative or non-integer entries in table") でエラー:
TRUE/FALSE が必要なところが欠損値です
> x.co <- corresp(x, nf=2, na.rm = TRUE)
if (any(x < 0 | x%%1 > 10 * sqrt(.Machine$double.eps))) warning("negative or non-integer entries in table") でエラー:
TRUE/FALSE が必要なところが欠損値です
上記テストデータから「のびのび ,20,18 ,NA」の行を削除すると欲しいグラフが出力されるのですが、corresp関数では欠損値は扱うことができないということでしょうか？caパッケージがインストールされていないこともあり、可能であればcorrespで実現させたいと考えています。ご回答いただけると助かります。よろしくお願いいたします。

[KojiDoi]

NAがあると確かにエラーで計算続行できないようですね。その行をあきらめるかデータを補完するかするしかないでしょう。たとえばこの例ではNAは0に置き換えて差し支えないように見えますね。
x[is.na(x)] <- 0

[fuusuke]

丁寧に回答いただきありがとうございました！！