Use vec(any(...)) instead of the sum:
julia> vec(any(P[2][:,2] .== permutedims(R[2]+1:R[2]+D[2]), dims = 2))
11-element BitVector:
1
1
1
1
1
1
1
1
1
0
0
julia> findall(vec(any(P[2][:,2] .== permutedims(R[2]+1:R[2]+D[2]), dims = 2)))
9-element Vector{Int64}:
1
2
3
4
5
6
7
8
9