https://imgur.com/HcH103o
я хочу вычислить эту формулу
где альфа ij - элемент матрицы A, размерность которого (n, m)
, где р - вектор, которыйдлина m
где d - вектор, длина которого m
я хочу получить результат d
вот мой плохой код
n <- 5
m <- 7
A <- matrix(runif(m*n,0,1),n,m)
d <- vector("numeric",length = m)
p <- vector("numeric",length = m)
p <- runif(m,0,1)
for(i in 1:n){
pk <- NULL
for(k in 1:m){
pk[k]<- p[k]*A[4,2]
}
denominator<- sum(pk)
d1[i]<- A[i,j]*p[j]/denominator
}
n <- 5
m <- 7
A <- matrix(runif(m*n,0,1),n,m)
d <- vector("numeric",length = m)
p <- vector("numeric",length = m)
p <- runif(m,0,1)
den <- NULL
densum <- NULL
densumsum <- NULL
d <- NULL
for( j in 1:m){
for(i in 1:n ){
for(k in 1:m){
den[k]<- A[i,j]*p[k]
}
densum[i]<- A[i,j]*p[j]/sum(den)
}
d[j]<- sum(densum)
}