ダメ出し：NA の存在するデータの多変量解析にご注意

Outlines of Factor Analysis (2); estimation / rotation でのことだけど

一番最後の

> 因子得点の基本的な特徴を再確認

> summary(result.fa$scores)

> ##       ML2               ML1        
> ##  Min.   :-2.2878   Min.   :-2.206  
> ##  1st Qu.:-0.6640   1st Qu.:-0.604  
> ##  Median :-0.1128   Median : 0.148  
> ##  Mean   :-0.0094   Mean   : 0.018  
> ##  3rd Qu.: 0.6678   3rd Qu.: 0.809  
> ##  Max.   : 2.2048   Max.   : 2.411  
> ##  NA's   :3         NA's   :3

因子得点の平均値が 0 になっていないのは，分析に NA を含むデータを除かなかったため。

fac_data <- na.omit(fac_data)

として，同じように分析を進めれば，以下のようにちゃんと平均値は0になる。fac_data では NA を含むデータが少ないので影響も小さいが，NA がたくさんあるデータだとどうなるかな?

これは，fa のバグなんだろうか?

> summary(result.fa$scores)
      ML2               ML1         
 Min.   :-2.2664   Min.   :-2.2769  
 1st Qu.:-0.6572   1st Qu.:-0.6321  
 Median :-0.1007   Median : 0.1331  
 Mean   : 0.0000   Mean   : 0.0000  
 3rd Qu.: 0.6697   3rd Qu.: 0.7946  
 Max.   : 2.1988   Max.   : 2.4392