Question

现在我正在实现pmf模型，如下所示：

pmf_code = """

data {

int<lower=0> K; //number of factors
int<lower=0> N; //number of user
int<lower=0> M; //number of item
int<lower=0> D; //number of observation
int<lower=0> D_new; //number of pridictor 
int<lower=0, upper=N> ii[D]; //item 
int<lower=0, upper=M> jj[D]; //user
int<lower=0, upper=N> ii_new[D_new]; // item
int<lower=0, upper=N> jj_new[D_new]; // user
real<lower=0, upper=5> r[D]; //rating
real<lower=0, upper=5> r_new[D_new]; //pridict rating

}

parameters {
row_vector[K] i[M]; // item profile
row_vector[K] u[N]; // user profile
real<lower=0> alpha;
real<lower=0> alpha_i;
real<lower=0> alpha_u;

}
transformed parameters {
matrix[N,M] I; // indicator variable
I <- rep_matrix(0, N, M);
for (d in 1:D){
    I[ii[d]][jj[d]] <- 1;
}
}
model {
for (d in 1:D){
    r[d] ~ normal(u[jj[d]]' * i[ii[d]], 1/alpha);
}

for (n in 1: N){
   u[n] ~ normal(0,(1/alpha_u) * I);
}
for (m in 1:M){
   i[m] ~ normal(0,(1/alpha_i) * I);
}
}
generated_quantities{
for (d in 1:D_new){
   r_new[d] <- normal(u[jj_new[d]]' * i[ii_new[d]], 1/alpha);
}
}
"""

但在这行代码中出现No matches for: real ~ normal(matrix, real)错误：

for (d in 1:D){
   r[d] ~ normal(u[jj[d]]' * i[ii[d]], 1/alpha);
}

但jj[d]应为int，表示user的ID。 u[int]应该是row_vector k因素，i[ii[d]]也是sum。它们的产物应该是一个真正的价值，为什么斯坦说它是一个矩阵？

提前致谢！

编辑：

通过添加No matches for: row vector ~ normal(int, matrix) for (n in 1: N){ u[n] ~ normal(0,(1/alpha_u) * I); }产品解决了初始问题。但现在发生了另一个错误：

在pmf模型中，(1/alpha_u) * I是单位矩阵。所以{{1}}的乘积也是一个矩阵。但斯坦只接受矢量或实际值作为方差。我想知道如何将其转换为矢量或单个值。感谢。

Answer 1

Stan是一种静态类型语言。 vector乘以row_vector会产生matrix。在这种情况下，它是1x1矩阵。有关返回类型operator*，请参阅手册v2.9.0版本的第355页。

而不是相乘，而是使用dot_product()。

将两个1 * k向量相乘，但得到无匹配：pystan中的real~normal（矩阵，实数）错误

1 个答案: