1

ロジスティック回帰で、もともと型 ts から変換された xts 時系列を使用していますが、

Error in `*.default`(x[good, , drop = FALSE], w) : non-conformable arrays

例:

success <- as.xts(ts(sample(0:10, 100, replace=T), start=1970, fr=12))
failure <- 10-success
x <- as.xts(ts(rnorm(100), start=1970, fr=12))
glm(cbind(success, failure) ~ x, family=binomial)

ただし、シリーズがクラスtsまたは単なるベクトルに変換されている場合は、エラーは発生しません。

glm(cbind(as.ts(success), as.ts(failure)) ~ as.ts(x), family=binomial)
glm(cbind(as.vector(success), as.vector(failure)) ~ as.vector(x), family=binomial)

元の xts シリーズで作業中にエラーを回避する方法はありますか?

4

1 に答える 1