Hallo,
ich bin im Rahmen einer Monte Carlo Simulation auf ein Problem gestoßen. Nachfolgend beschreibe ich es. Hoffentlich könnt ihr mir folgen!
zunächst habe ich bei R Studio eine Datenmatrix mit 100.000 Zeilen und 4 Spalten konstruiert. Jede Zeile besteht aus einer abhängigen und 3 unabhängigen Variablen. Dann habe ich eine Regression der Datenmatrix durchgeführt. Die daraus resultierenden Beta Werte meiner unabhängigen Variablen nenne ich nachfolgend "wahre" Beta Werte.
Dann habe ich 1000 Stichproben aus meiner Datenmatrix a jeweils 25 Zeilen gezogen.
Jetzt möchte ich jede einzelne Stichprobe regressieren und die jeweiligen Beta Werte auf Wahrheit überprüfen. Dafür möchte ich einen t Test benutzen mit der Nullhypothese H0: Beta Wert dieser Regression = wahrer Beta Wert. Jede Stichprobe der insgesamt 1000 soll regressiert und auf Wahrheit getestet werden. Letztlich interessiert mich die prozentuale Ablehnungsquote der t Tests in den 1000 Stichproben.
Meine Idee war daher eine Schleife mit einem Zähler zu konstruieren, der jede einzelne Ablehnung der Nullhypothese zählt. Letztlich sollte mir die Schleife dann die gesamte Anzahl der abgelehnten t Tests "ausspucken". Leider bin ich nach langen probieren nicht in der Lage die Schleife zu konstruieren. Hier ein Auszug aus meinen Quellcode:
# generate 1000 random samples without replacement for n= 25
s25 <- mydata[sample(nrow(mydata), 25, replace=F), ]
r25 <- replicate(1000, s25, simplify = F)
Also ganz konkret möchte ich "r25" in die Schleife packen und jedes einzelne 25-teilige Sample regressieren (einschließlich t Test und einen Zähler falls es abgelehnt worden ist). Ist das überhaupt möglich? Ich bin für jeden Beitrag dankbar!
Mit besten Grüßen
Moritz