The sample data has the following format:

data$u - sample u; data$v - sample v. 

1. Create variable joint_distribution- matrix of joint distribution P[u,v].
2. Create variable u_Marginal equal to the marginal distribution P[u].
3. Create variable v_Marginal equal to the marginal distribution P[v].
4. Create variable u_Conditional_v equal to the vector of conditional probabilities P[u|v=4].
5. Create variable v_Conditional_u equal to the vector of conditional probabilities P[v|u=3].

Collect all results in variable res:

res <-list(Joint_distribution=joint_distribution, u_Marginal = u_Marginal, v_Marginal = v_Marginal, u_Conditional_v = u_Conditional_v, v_Conditional_u = v_Conditional_u )

u "v"
2 3
2 3
2 2
1 2
1 2
1 2
1 3
3 2
2 3
2 3
2 3
1 3
2 2
2 2
2 3
1 4
1 2
1 4
2 4
2 2
2 2
2 3
2 3
1 3
1 3
2 1
2 3
3 2
1 2
1 3
1 4
1 2
1 4
3 3
2 3
2 3
1 2
3 3
2 3
1 2
1 4
1 2
1 3
2 2
1 3
3 2
3 3
1 2
2 2
1 4
2 3
1 3
2 2
2 3
1 4
1 1
1 3
3 2
1 4
2 4
1 2
1 4
2 3
3 3
1 1
2 2
1 4
1 2
1 4
1 4
2 3
1 4
1 3
2 3
3 2
3 2
1 3
2 4
2 3
1 3
2 3
1 2
2 2
3 3
2 4
1 3
2 2
2 2
1 4
3 1
1 4
1 3
1 4
1 2
1 1
1 3
2 4
3 3
2 4
1 1