zweiseitige ks Testschleife, erhalten p.value

Question

Ich habe eine Spalte von Daten, aus denen ich randomisierte Unter Stichproben von 50% nehme. Ich führe einen zweiseitigen ks-Test, um die Verteilung von 50% der Daten gegen 100% der Daten zu vergleichen, um zu sehen, ob die Verteilung noch eine signifikante Anpassung ist.

Um meine Ziele zu erreichen, möchte ich dies als eine Schleife von sagen 1000 ausführen, um einen durchschnittlichen p-Wert von 1000 randomisierten Unterproben zu erhalten. Diese Codezeile gibt mir einen einzigen p-Wert für eine zufällige Teilmenge von 50% meiner Probe:

dat50=dat[sample(nrow(dat),replace=F,size=0.50*nrow(dat)),] 
ks.test(dat[,1],dat50[,1], alternative="two.sided") 

Ich brauche eine Codezeile, die diese 1000-mal ausgeführt wird den resultierenden (verschiedenen) p-Wert Speicher jeden Zeit in einer Spalte, die ich dann mitteln kann. Der Code, den ich versuche, zur Arbeit zu kommen sieht wie folgt aus:

x <- numeric(100) 
for (i in 1:100){ 
    x<- ks.test(dat[,7],dat50[,7], alternative="two.sided") 
    x<-x$p.value 
} 

jedoch diese nicht speichert mehr p-Wert

Auch dies versucht:

get.p.value <- function(df1, df2) { 
    x <- rf(5, df1=df1, df2=df2) 
    p.value <- ks.test(dat[,6],dat50[,6], alternative="two.sided")$p.value 
} 
replicate (2000, get.p.value(df1 = 5, df2 = 10)) 

Ich hoffe, das ist klar und ich würde jede Hilfe bei der Lösung dieses Problems sehr schätzen!

Q

Answer 1

In Ihrer for-Schleife Sie x in jeder Iteration überschreiben was bedeutet, dass Sie nur den p-Wert für die letzte Iteration sparen. Versuchen Sie stattdessen:

x <- numeric(100) 
for (i in 1:length(x)) 
    x[i] <- ks.test(dat[,17], dat[sample(nrow(dat), replace=F, size=0.5*nrow(dat)),7])$p.value 

Sie das gleiche Ergebnis mit bekommen replicate mit:

replicate(100, ks.test(dat[,7], dat[sample(nrow(dat), replace=F, size=0.5*nrow(dat)),7])$p.value)