0

Rで実行中のウィンドウのスピアマン相関を計算したいと思います。これまでのところ、パッケージrunningから使用してきましたがgtools、ピアソン相関しか取得できませんでした。パラメータを変更しようとしましたfunが、エラー以外は何も生成できませんでした。

私の入力はdata.frameです：

row.names   Small   Large
1   1   97.80341    88.71192
2   2   97.46807    87.96206
3   3   97.18862    88.13904
4   4   97.76615    87.67329
5   5   97.09081    87.52425
6   6   97.16067    87.85493
7   7   97.73820    88.43712

など。これは、私が操作しようとした基本的な実行コマンドです。

corr_table <- running (mytable$Large, mytable$Small,fun=cor, width=20, by=10, allow.fewer=FALSE, pad=FALSE, align="left")

2 番目の質問は、各「ウィンドウ」の pvalue を追加するにはどうすればよいですか?

ありがとう！

4

1 に答える 1

0

OK は、両方の質問に対する解決策を見つけました。

corr_table <- running (mytable$Large, mytable$Small,fun=cor.test, method="spearman", width=60, by=30, allow.fewer=FALSE, pad=FALSE, align="left")
newlist <- t(corr_table)
newlist <- as.matrix (newlist)
y <- newlist[,4]
x <- seq(from=1, to=92) #change the "to" according to number of values in Y
pvalue <- newlist[,3]
plot(x,y, xlim=c(1,92),ylim=c(-1,1), col=ifelse(pvalue < 0.05, "red", "black")) #change the "xlim" according to number of values in Y
于 2015-06-22T11:42:41.107 に答える