Есть много способов, таких как:
p = fill(NaN, n, periods)
p[:, 1] .= 0.5
или
p = [j == 1 ? 0.5 : NaN for i in 1:n, j in 1:periods]
или аналогично:
p = [ifelse(j == 1, 0.5, NaN) for i in 1:n, j in 1:periods]
Все они имеют преимущество, заключающееся в том, что они выделяют толькоодна матрицаПервый из них немного быстрее, но требует двух операторов.