Ich habe eine Zufallsvariable X
und eine Transformation f
und ich würde gerne die Wahrscheinlichkeitsverteilungsfunktion von , zumindest ungefähr kennen. In Mathematica gibt es TransformedDistribution, aber in R. konnte ich nichts Ähnliches finden. Wie ich schon sagte, wäre auch eine Art Näherungslösung in Ordnung.Hat R etwas Ähnliches wie TransformedDistribution in Mathematica?
3
A
Antwort
3
Sie können das distr
Paket überprüfen. Zum Beispiel, sagen, dass y = x^2+2x+1
, wo x
normalerweise mit einem Mittelwert 2 und Standardabweichung verteilt 5. Sie können:
require(distr)
x<-Norm(2,5)
y<-x^2+2*x+1
#[email protected] gives random samples. We make an histogram.
hist([email protected](10000))
#[email protected] and [email protected] are the density and the cumulative functions
[email protected](80)
#[1] 0.002452403
[email protected](80)
#[1] 0.8891796
gut aussieht. Aber ich bekomme y @ p (-1) = 0,08723655 das ist komisch, da y niemals negativ sein sollte ... – frank
Ich ging in das hinein und es sieht so aus, dass der Ausdruck 'x^2 + 2 * x + 1' interpretiert wird als die Summe der * unabhängigen * Zufallsvariablen, dh der 'x^2'-Teil und der' 2 * x'-Teil werden getrennt voneinander genommen. Wenn Sie "y <- (x + 1)^2" definieren, erhalten Sie das richtige Verhalten. Schätze, dass dies die Aufmerksamkeit des Betreuers erfordert. – nicola
Super, danke für die ganze Arbeit. Ungeachtet des Fehlers war genau das, wonach ich suchte. – frank