[程式] R 軟體出現引數不正確訊息
------------------------------------------------------------------------
[軟體程式類別]: R
[程式問題]: 出現'times'引數不正確訊息
[軟體熟悉度]:低
[問題敘述]:小弟在處理一個MORAN PROCESS的問題，其中會用到一個矩陣
若該矩陣元素為純數字時則無問題，但改成變數時就出現'times'引數不正確
想將矩陣A裡面的元素改為N1 N2 N3變數取代
若改成變數執行時，到下列程式碼see myMoran那行後就出現該訊息了@@
請教大大們該如何解決
詳見小弟程式碼如下:
[程式範例]:
myMoran <- function(x,A,u,n) {
tmp <- matrix(nrow=n+1,ncol=length(x))
tmp[1,] <- x
for (i in 1:n) {
f <-
c(sum(A[1,]*x)-A[1,1],sum(A[2,]*x)-A[2,2],sum(A[3,]*x)-A[3,3])/(sum(x)-1)
p <- round((x*f/sum(x*f))*x)
k1 <- sample(rep(1:3,p),1)
if (runif(1) < u) k1 <- sample(c(1:3)[-k1],1)
k2 <- sample(rep(1:3,x),1)
x[k1] <- x[k1] + 1
x[k2] <- x[k2] - 1
tmp[i+1,] <- x
}
return(tmp)
}
N <- 30
N1 <- 14
N2 <- 3
N3 <- N - N1 - N2
x <- c(N1,N2,N3)
u <- 0.1
n <- 1000
see <- myMoran(x,A,u,n)
-----------------------------------------------------------------------------
--
離開你的視線~讓情緒完完整整的宣洩~
我想今夜就這樣吧~
反正孤獨也無所謂~
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 114.34.118.44
→
07/27 03:45, , 1F
07/27 03:45, 1F
→
07/27 03:46, , 2F
07/27 03:46, 2F