1

I would like to simulate random variates from a transformed distribution of a joint distribution and a constant.

i.e.

joint=ProductDistribution[NormalDistribution[],BetaDistribution[1,2],BetaDistribution[2,2]];
transform=TransformedDistribution[3*joint....]

I am unsure how to complete the transformed distribution code to put in the variables.

Jim
  • 93
  • 7

2 Answers2

1
joint = ProductDistribution[NormalDistribution[], BetaDistribution[1, 2], BetaDistribution[2, 2]];
q = TransformedDistribution[{x, y, z} w, {{x, y, z} \[Distributed]  joint, 
                                          w \[Distributed] UniformDistribution[]}]
ListPlot3D@RandomVariate[joint, 100]

Mathematica graphics

Dr. belisarius
  • 115,881
  • 13
  • 203
  • 453
1
joint = ProductDistribution[NormalDistribution[], BetaDistribution[1, 2], BetaDistribution[2, 2]];
td = TransformedDistribution[3 {x, y, z} , Distributed[{x, y, z}, joint]];

Through@{Mean, Variance}@td
(* {{0,1,3/2}, {9,1/2,9/20}} *)

PDF[td, {x, y, z}]

enter image description here

ListPlot3D[RandomVariate[td, 50]]

enter image description here

See also: this answer to a closely related question by the same OP.

kglr
  • 394,356
  • 18
  • 477
  • 896