Hi,
I have a matrix with 4 rows (different groups), and 10,000 columns. Each cell contains a p-value.
In addition, I have a vector named alpha, with 1000 values from 0 to 1.
I need to compare each alpha value, to the 10,000 p-values of the first matrix raw.
I need to count how many p-values are smaller than the alpha.
This is my code:
alpha<-seq(from=0,to=1,by=0.001)
for (i in 1:length(alpha)) {
treshold<-alpha[i]
for (j in 1:N) {
my.true<-p.values[1,j]<treshold
results<-matrix(my.true, nrow = 1000, ncol = N, byrow = TRUE)[i,j]
}
}
Mi idea was to get TRUE or FALSE for each value, and then count how many TRUEs I got for each alpha. Finally divide these values by 10,000 in order to get the percentage of the p-values smaller than alpha.
My code doesn't work. It seems like it runs for the first alpha and doesn't continue to the following one.
Please assist 