I need to perform 10000 T tests
#I have two populations with different means
Popc1<-rnorm(100000,10,2)
Popc2<-rnorm(100000,8,2)
#I created two sets of samples - each set has 10000 samples, and I made a
matrix of 20 rows and 10000 columns to fit the data
sampc1<-matrix(,20,10000)
for(j in 1:10000){sampc1[1:20,j]<-sample(Popc1,20)}
sampc2<-matrix(,20,10000)
for(j in 1:10000){sampc2[1:20,j]<-sample(Popc2,20)}
#I wrote a blank matrix with 1 row and 10000 columns. I am planning to fill
this matrix with the t values that i calculate for each pair of samples
t<-matrix(,1,10000)
#What i need to do is to calculate the t value for each pair of samples: So
if i were to do it manually, for the first sample in each set, I would do
something like:
t.test(sampc1[,1],sampc2[,1])
#And i would continue doing this all the way until the 10000th pair of
samples
so t.test(sampc1[,2],sampc2[,2])
t.test(sampc1[,3],sampc2[,3]) and so on
#But i need to do this for all 10000 pairs of samples:
#I tried to write a loop
for(j in 1:10000) {t[,1:10000]<-t.test(sampc1[, j],sampc2[, j])$statistic}
#but this fills each column in the matrix with the same exact t value.
#Can someone tell me how to write the code to perform each pair of t tests?
I think i am missing something in my loop
Thanks
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
