J'essaye ce code R pour faire sortir le nombre de gènes qui ne suivent pas la loi normale pour un échantillon
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
#tester la normalité des variables (gènes) en utilisant
#le test de kolmogorov
N <- 7129
n <- 38
R <- vector()
for (i in 1:N) {
  x <- as.matrix(data[i,])
  test<- ks.test( x, rnorm(n))  
  R <- append(R, test$p.value )
  }
sum(R<.05)
et j'ai en retour
> sum(R<.05)
[1] 7129 # càd toute la liste

Je m'attends plutôt à ~3000
Voyez-vous ce qui ne va pas ?
Merci bien