Lognormalverteilung
Hallo,
Hatte vielleicht jemand schon mal eine ähnliche Anforderung bzw. ein ähnliches Problem:
Ich möchte lognormalverteilte Zufallszahlen erzeugen, die dann einen bestimmten Mittelwert und eine bestimmte Standardabweichung haben sollen.
Wenn ich folgendermaßen direkt mit der "lognormal-Funktion" simuliere, bekomme ich einen Mittelwert von 1.6565948
und eine Standardabweichung von 2.1783043.
do i=1 to 10000;
x=rand('lognormal');
output test;
end;
run;
Kann ich diese Funktion irgendwie direkt umformen, um einen bestimmten Mittelwert und eine bestimmte Standardabweichung zu bekommen?
Ich denke, dass ich dazu eher über eine Normalverteilung simulieren muss. Leider funktioniert es nicht so einfach wie ich es im Folgenden probiert habe:
mu=10;
sd=33;
do i =1 to 10000;
x= exp(mu+sd*rannor(95));
output test2;
end;
run;
Denn so lege ich die Parameter für die Standardabweichung fest und für die Lognormalverteilung bekomme ich einen Mittelwert von 4.7686594E52 und eine Standardabweichung von 4.6720425E54.
Es bleibt also die Frage wie ich die Parameter für die Standardnormalverteilung wählen muss, damit ich eine Lognormalverteilung mit gewünschtem Mittelwert und Standardabweichung bekomme.
Vielen Dank im Voraus für eure Inputs!
- Anmelden oder Registrieren um Kommentare zu schreiben
Mit der PDF-Funktion können
Mit der PDF-Funktion können Sie ihr Problem glaube ich lösen
(siehe SAS-Online Doc):
Syntax
PDF('LOGNORMAL',x<,,>)
where
x
is a numeric random variable.
is a numeric location parameter.
Default: 0
is a numeric scale parameter.
Default: 1
Range: > 0
The PDF function for the lognormal distribution returns the probability density function of a lognormal distribution, with the location parameter and the scale parameter . The PDF function is evaluated at the value x. The equation follows: