作者PhysicsBoy (翟南)
看板Statistics
标题[程式] R 软体出现引数不正确讯息
时间Sun Jul 24 02:35:33 2011
------------------------------------------------------------------------
[软体程式类别]: 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
1F：→ gsuper:太多变数了 难怪没人理你 问问题前要先简化问题 07/27 03:45
2F：→ gsuper:而且你没秀出 "A" 里面长什麽样子 07/27 03:46