web-dev-qa-db-ja.com

担当者のエラー：「回」引数が無効です

次のコードを10000回繰り返し実行しようとすると、次のエラーが発生します。rep（G1 [、2]、G1 [、3]）のエラー：無効な 'times'引数。そのため、そのエラーを修正するためにコードを変更する方法がわかりません。基本的に、年間の故障時間と修理時間8736時間の式を使用して、発電機のパフォーマンスの時系列を作成したいので、発電機が稼働していない時の時系列を時系列で表示します。開始条件は、発電機が最初の1時間稼働していることです。確かにこれをシミュレートするためのより洗練された解決策がありますが、私はそれを見つけることができません。コメントやヘルプをいただければ幸いです。

MTTF<-2940 # MEDIUM TIME TO FAIL(hours)
MTTR<-60 # MEDIUM TIME TO REPAIR (hours)
TTF<--MTTF*log(runif(100))# equation for Time to fail 
TTR<--MTTR*log(runif(100))# equation for Time to repair
mix<-rep(0,length(TTF)+length(TTR))
sw<-rep(0,length(TTF)+length(TTR))
for(i in 1:length(TTF)){
mix[2*i-1]<-TTF[i]
sw[2*i-1]<-1
mix[2*i]<-TTR[i]
}
cmix<-cumsum(mix)
ccmix<-cbind(cmix[1:which(cmix>8736)],sw[1:which(cmix>8736)])
ccmix[dim(ccmix)[1],1]<-8736


G1<-round(ccmix)
# transform binary values
G1[G1 == 1] <- 12 # is the capacity of the generator

G1 <- cbind(G1, c(G1[1,1], diff(G1[,1])))
a1 <- rep(G1[,2], G1[,3]) ## GENERATING 8736 Values

したがって、望ましい出力は、ONのとき12の8736値であり、OFFのとき0です。

rtime-seriesmontecarlo
6
kelamahim

何を確認するG1[,3]です。 times引数の負の値が原因でエラーが発生する可能性があります

10
Adela