Advertisement
Guest User

Untitled

a guest
Sep 26th, 2017
48
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.47 KB | None | 0 0
  1. data1 = RandomVariate[NormalDistribution[5, 1.2], 1000];
  2. skd1 = KernelMixtureDistribution[data1];
  3. custom[a_] := ProbabilityDistribution[PDF[skd1, x+a], {x, 0, 10}, Assumptions -> a >= 0 && a < 2.5];
  4.  
  5. data3 = RandomVariate[NormalDistribution[6, 1.2], 1000];
  6. FindDistributionParameters[data3, custom[a], ParameterEstimator -> {"MaximumLikelihood", Method -> {"FindMaximum"}}]
  7.  
  8. custom[a_] := ProbabilityDistribution[a PDF[skd1, x], {x, 0, 10}, Assumptions -> a >= 0 && a < 2.5];
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement