Hi James,

Based on your suggestion, I rewrote the code as:

```
transformed parameters {
vector[N] mu;
vector[N] lambda_tweedie;
vector[N] alpha_tweedie;
vector[N] beta_tweedie;
real p_star = log((p-1)/(2-p));
for (k in 1:s) {
for (j in 1:n_s){
mu[(k*n_s-n_s+j)] = exp(alpha[k] + X[(k*s-s+j), 1]*beta);
lambda_tweedie[(k*n_s-n_s+j)] = 1/phi*mu[(k*n_s-n_s+j)]^(2-p)/(2-p);
alpha_tweedie[(k*n_s-n_s+j)] = (2-p)/(p-1);
beta_tweedie[(k*n_s-n_s+j)] = 1/phi*mu[(k*n_s-n_s+j)]^(1-p)/(p-1);
}
}
}
model {
// other prior definition
p_star ~ normal(0, 10);
alpha ~ multi_normal_prec(zeros, delta*(T - rho*W_orignal));
for (k in 1:s) {
for (j in 1:n_s){
//mu[(k*n_s-n_s+j)] = exp(alpha[k] + X[(k*s-s+j), 1]*beta);
//lambda_tweedie[(k*n_s-n_s+j)] = 1/phi*mu[(k*n_s-n_s+j)]^(2-p)/(2-p);
//alpha_tweedie[(k*n_s-n_s+j)] = (2-p)/(p-1);
//beta_tweedie[(k*n_s-n_s+j)] = 1/phi*mu[(k*n_s-n_s+j)]^(1-p)/(p-1);
}
}
}
}
```

Here, alpha, beta, p, phi, and rho are all unknown parameters I want to figure out. Eventually, I received the following error messages:

```
SYNTAX ERROR, MESSAGE(S) FROM PARSER:
Base type mismatch in assignment; variable name = mu, type = real; right-hand side type=vector
error in 'model4d2c54af6689_model' at line 46, column 31
-------------------------------------------------
44: for (k in 1:s) {
45: for (j in 1:n_s){
46: mu[(k*n_s-n_s+j)] = exp(alpha[k] + X[(k*s-s+j), 1]*beta);
^
47: lambda_tweedie[(k*n_s-n_s+j)] = 1/phi*mu[(k*n_s-n_s+j)]^(2-p)/(2-p);
-------------------------------------------------
PARSER EXPECTED: <expression assignable to left-hand side>
Error in stanc(file = file, model_code = model_code, model_name = model_name, :
failed to parse Stan model 'model' due to the above error.
In addition: Warning message:
In readLines(file, warn = TRUE) :
```

What should I do next? Thank you for your time.

With appreciation,

Tairan