Can someone please help me with the coding (in R) and the answers for this question? Thank you.

Task 1 (Discrete Joint Distribution) Consider the joint probability mass function of a pair of random variables (X, Y), px,y(i, j) = P(X = i, Y = j) for i, j = 1, ...,6. In R, we take px,r(i, j) = pXY [i, j] where pXY