必要な複製の数を計算する関数をrで作成しました。私の問題は、デルタとシグマの値を行列に入れてから関数に入れる必要があることです。シグマ/デルタにおけるRmarkdownでの行列の問題
required.replicates <- function(delta, sigma, z.alpha = 1.959964,
z.beta=0.8416212)
{
r.exact <- (2*(z.alpha + z.beta)^2)*(sigma/delta)^2
result <- list(delta, sigma, z.alpha, z.beta)
names(result) <- c("delta", "sigma", "z.alpha", "z.beta")
r.integer <- ceiling(r.exact)
replicate.output <- list(result, r.exact, r.integer)
names(replicate.output) <- c("result","r.exact","r.integer")
return(replicate.output)
}
sigma.vec <- c(2,4,6,8,10,12)
delta.vec <- c(1,2,5,8,10)
sigma.mat <- matrix(rep(sigma.vec,times=length.s), nrow =length.s,
byrow=TRUE)
delta.mat <- matrix(rep(delta.vec,times=length.s), ncol = length.s,
byrow=FALSE)
required.replicates(delta= delta.mat, sigma=sigma.mat)
エラー:非適合配列
私は、エラーメッセージが出てい理由はわかりません。
あなたは何をしたいですか? 'sigma.mat'は'(10 * 6) 'の行列ですが、' delta.mat'は '(5 * 10)'なので、 'sigma.mat/delta.mat'が生成するものは分かりません。 – herbaman
申し訳ありません。私がする必要があるのは、sigma.vecとdelta.vecのすべての値をrequired.replicatesの式に入れることです。 –