[R言語]ポアソン分布乱数を使った2変数サンプル十万個(サイズ9)の相関係数

前提

R言語での計算です。

実現したいこと

ポアソン分布乱数を使って、xとyという２つの変数のサンプルサイズ9のサンプルを十万個生成し、サンプルの相関係数の平均と分散を計算したい。

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

警告メッセージ: cor(x1,y1)で、標準偏差が０です。

該当のソースコード

R

kN1<-(9) #サンプルサイズkNs<-(100000)　#サンプル数lamda<-(2)　#ポアソン分布パラメータcors<-numeric(kNs)for(i in seq(kNs)){+x1<-rpois(kN1,lamda)+y1<-rpois(kN1,lamda)+cors[i]<-cor(x1,y1)}mean(cors)　#相関係数の平均sd(cors)　#相関係数の分散

試したこと

サンプル数(kNs)が1000,10000,50000のケースは計算が可能であった。

補足情報（FW/ツールのバージョンなど）

オープンソースのフリーソフトウェアR